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ABSTRACT 

Gimballess  inertial  reference  equipment  might  provide  six  possible 
advantages  over  conventional  stable  member  systems.     These  are  savings 
in  power,    size,   weight,    and  expense,   flexibility  of  packaging,    and  ease 
of  maintenance.     The  theoretical  problems  involved  in  gimballess  equip- 
ment are  the  question  of  what  inertial  instruments  measure,    and  how 
best  to  use  the  information  they  generate.     The  engineering  problems 
concern  the  difficulty  of  torquing  gyroscopes  with  accuracies  approaching 
one  part  in  a  hundred  million,    and  the  design  of  inertial  instruments  :o 
withstand  the  high  angular  rates  the}r  will  experience  in  this  application. 
The  inertial  package  to  be  considered  contains  three  delta-modulated  inte- 
grating gyros  and  three  delta-modulated  integrating  accelerometers. 

Delta  modulation  is  a  pulse  modulation  system  which  delivers 
output  pulses,    representing  increments  of  the  input  or  one  of  its 
integrals  or  derivatives,    in  synchronism  with  a  clock.     The  analysis 
of  delta-modulated  loops  with  either  a  two-  or  three-level  relay  is 
a  straightforward  matter  using  z-transforms  and  modified  step-by-step 
transient  methods.     Application  of  delta  modulation  to  integrating  gyros 
and  accelerometers  leads  to  instruments  whose  output  pulses  represent 
increments  of  rotation  and  velocity. 

Transformation  of  a  vector  from  one  Cartesian  coordinate  frame 
to  another  may  be  represented  in  three  essentially  different  ways: 
through  matrix-tensor  operations,    four-parameter  techniques,    and 
vector  algebra.     Comparison  of  these  to  determine  which  is  most 
suitable  for  machine  computation,   both  from  the  standpoint  of  generating 
the  transformation  parameters  and  of  accomplishing  the  transformation, 
shows  that  a  matrix  transformation,    with  straightforward  generation  of 
direction  cosines,    is  best.   The  most  attractive  choice  of  machine  for 
the  computation  is  the  digital  differential  analyzer,    since  the  angular  and 
linear  motion  information  is  in  incremental  form.   High-order  integration 
rules  with  larger  angular  increments  cannot  be  used  to  reduce  computer 
speed  while  maintaining  accuracy  because  of  an  irreducible  uncertainty 
in  the  digital  angular  information., 


The  transformation  computer  may  be  analyzed  using  frequency  domain 
techniques.     With  these  techniques,    a  signal  flow  graph  representing  the 
computer  can  be  developed,  which  gives  great  insight  into  the  dominant 
mode  of  operation  of  the  computer,   and  which  provides  an  estimate  of 
truncation  error.     Quantization  may  be  treated  statistically. 

This  investigation  has  considered  one  way  of  instrumenting  gimballess 
inertial  reference  equipment  in  detail,   and  has  provided  a  survey  of  some 
of  the  possibilities  in  this  field.     It  has  shown  the  theoretical  feasibility 
of  such  equipment,   and  has  also  indicated  that  the  major  development 
effort  required  to  permit  practical  realization  of  such  equipment  with 
a  specified  accuracy  must  be  concerned  with  inertial  instruments  designed 
especially  for  this  application. 
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Associate  Professor  of 
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CHAPTER  1 


INTRODUCTION 


LA    Background 

The  problem  of  arriving  at  a  destination  has  been  important  to  man 
from  the  beginning.     At  first,   getting  there  was  merely  a  matter  of 
following  a  line  of  sight  to  the  destination.     Later,   landmarks  whose 
position  with  respect  to  the  destination  were  known  were  used  as 
guideposts  in  making  journeys  out  of  sight  of  the  destination.     The  stars, 
which  can  be  used  out  of  sight  of  land,    represent  an  important  class  of 
these  reference  points.     Celestial  bodies  can  be  observed  only  on  clear 
nights,   and  of  course   other  landmarks,  by  and  large,   must  be  visible 
to  be  useful.     With  the  invention  of  radio,   and  later  radar,   man-made 
radiation  could  be  used  to  distinguish  landmarks  during  periods  of 
reduced  visibility.     In  modern  military  and  transport  applications  many 
situations  exist  in  which  terrestrial  and  celestial  reference;    points 
are  either  unavailable  or  unreliable  as  a  basis  for  guidance.     In  such 
cases  there  is  a  need  for  navigation  systems  which  do  not  depend  on 
any  external  information.     Inertial  navigation  systems,   which  deduce 
the  position  of  a  vehicle  from  on-board  measurements,    are  systems 
of  this  type.     The  measured  acceleration  is  integrated  twice  with  respect 
to  time  to  give  displacement  for  an  initial  position. 

The  majority  of  inertial  navigation  systems  today  use  two  or  three 
single-axis  accelerometers  mounted  on  a  stable  platform  to  measure 
acceleration.     The  platform  is  either  kept  non-rotating  with  respect 
to  inertial  space  or  precessed  so  that  the  input  axes  of  the  accelerometers 
are  kept  coincident  with  some  slowly  rotating  set  of  coordinate  axes, 
such  as  local  geographic  coordinates.     The  stabilized  platform  is 


articulated  with  respect  to  its  support  by  a  system  of  gimbals.     The 
gyros  mounted  on  the  platform  detect  the  angular  motion  of  the  platform, 
and  the  gimbals  are  torqued  in  response  to  gyro  signals  to  drive  the 
platform  back  to  its  null  position. 

Another  configuration  for  inertial  navigation  systems  has  the  inertial 
instruments  fixed  to  the  vehicle  and  hence  is  gimballess.     Single  axis 
accelerometers  measure  vehicle  acceleration  and  resolve  it  into  vehicle 
coordinates.     The  guidance  equations  may  be  solved  in  body  coordinates, 
in  which  case  gyro  information  about  vehicle  rotations  is  used  directly 
in  the  guidance  equations;  or  they  may  be  solved  in  an  inertial  coordinate 
frame,    in  which  case  the  gyro  information  is  used  to  generate  a  coordinate 
transformation  which  operates  on  the  accelerometer  data  to  give  acceleration 
resolved  into  inertial  axes. 

There  are  six  possible  advantages  available  in  gimballess  guidance 
systems.     With  the  removal  of  gimbals,    weight,    space,   and  power 
required  are  reduced.     Further,   the  not  inconsiderable  expense  of 
gimbals  is  eliminated.     Since  the  inertial  package  in  a  gimballed  system 
must  assume  practically  a  spherical  shape,    it  is  necessary  to  find  a 
spherical  space  in  the  vehicle  in  which  to  put  the  guidance  system. 
Without  gimbals  to  dictate  packaging,   much  more  flexibility  in  packaging 
is  available.     Finally,  assembly  and  maintainence  of  the  inertial  package 
is  easier  without  the  encumbrance  of  gimbals. 
l.B    The  Problem  of  Gimballess  Inertial  Reference  Equipment 
1.  B.  1    Theoretical  Aspects 

The  theoretical  problems  of  gimballess  inertial  navigation  are 
two:  the  measurement  of  acceleration  in  body  fixed  coordinates  and 
the  determination  of  the  orientation  of  the  body  fixed  coordiante  frame 
with  respect  to  the  inertially  fixed  coordinate  frame  from  gyro  information. 

The  three  single  axis  accelerometers  measure  inertial  acceleration 
and  resolve  it  into  components  along  body-fixed  axes.     Coriolis  effects 
are  not  inherently  present  in  the  measured  acceleration  components; 
they  arise  only  in  the  intepretation  of  these  components  in  non-inertial 
reference  frames.     If  such  frames  are  used  in  the  guidance  equations, 
the  Coriolis  terms  must  be  computed  separately  using  angular  information 
from  the  gyros.     If  the  guidance  equations  are  set  in  inertial  coordinates, 
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the  three  measured  components  need  only  be  transformed  into  the  inertial 
coordinate  system. 

The  generation  of  a  vector  transformation  from  gyro  information  is 
a  straightforward  process.     The  differential  equation  governing  the 
transformation  parameters  may  be  derived  and  instrumented  in  any 
convenient  way. 
1.  B.  2    Engineering  Aspects 

The  engineering  problems  involved  in  gimballess  applications  of 
inertial  guidance  are  many,   but  the  two  most  important  involve  the 
environment  to  which  the  inertial  instruments  will  be  subjected  and  the 
speed  of  computation  required  to  keep  the  computed  orientation  of  the 
body -fixed  coordinate  frame  close  to  its  actual  orientation. 

While  inertial  instruments  mounted  on  a  stable  platform  are  isolated 
from  the  angular  motions  of  the  vehicle  carrying  them,   they  do  experience 
all  the  linear  motion  and  most  of  the  vibrations  of  the  vehicle.     Body- 
mounted  instruments,    in  contrast,   experience  both  the  angular  and  linear 
motions,  but  they  may  be  effectively  shock  mounted  to  minimize  the 
vibrations  present.     Although  much  effort  has  been  expended  in  developing 
instruments  whose  performance  does  not  deteriorate  under  acceleration, 
practically  none  has  been  devoted  to  designs  which  might  experience  large 
angular  velocities.     For  this  reason  no  instruments  are  available  at 
present  which  can  be  used  in  a  gimballess  application. 

In  gimballess  inertial  reference  equipment,   the  gyros  are  used  to 

measure  angular  orientation  of  the  body-fixed  coordinate  frame  with 

respect  to  some  reference.     This  may  be  done  either  by  using  a  pick-off 

arrangement  with  free  gyros,   or  by  torquing  single-degree-of-freedom 

gyros  to  make  them  follow  body  rotations.     If  free  gyros  are  used,   the 

pick  -off s  must  be  accurate  to  about  five  thousandths  of  a  degree.     The 

required  accuracy  of  torquing  for  single-degree-of-freedom  gyros  may 

be  shown  by  the  following  considerations.     If  the  gyro  is  to  have  a  drift 

_  8 
rate  of  one  meru  (about  7  .10       radians  per  second),   the  uncertainty  in 

applied  torque  must  not  exceed  this  level.     The  maximum  angular  velocity 

which  a  typical  vehicle  might  experience  is  about  one  radian  per  second, 

so  that  the  maximum  torque  available  must  be  able  to  precess  the  gyro 


21 


at  that  rate.     Comparison  of  these  two  numbers  shows  that  for  gimballess 
applications,    and  the  above  specifications,   the  torquing  must  be  accurate 
to  about  seven  parts  in  a  hundred  million.     For  higher  acceptable  drift 
rates  or  lower  maximum  velocities,;  this  accuracy  requirement  is  less 
stringent. 

The  orientation  of  the  body-fixed  coordinate  frame  must  be  known 
to  a  high  degree  of  accuracy  (about  five  thousandths  of  a  degree);  thus 
at  reasonable  angular  velocities  very  high  computer  speeds  are  required 
to  keep  up  with  the  vehicle  rotation.     If  a  single  adder  is  to  be  used, 
add  times  on  the  order  of  five  microseconds  per  word  are  required. 
This  is  near  the  limit  of  available  performance. 
1.  C    Scope  of  the  Investigation 

A  gimballed  inertial  reference  system  delivers  acceleration  information 
in  stabilized  inertial  coordinates.     This  investigation  will  consider  an 
inertial  reference  system  with  single  axis  integrating  accelerometers 
and  single-degree-of-freedom  integrating  gyros  fixed  to  a  vehicle.     All 
the  instruments  will  be  pulse  torqued.     The  gyro  information  is  to  be 
used  to  generate  a  coordinate  transformation  which  will  operate  on 
accelerometer  information  to  give  inertially  stabilized  data.     Thus  it 
is  proposed  that  the  gimballed  inertial  package  be  replaced  by  a  body- 
mounted  inertial  package  and  a  transformation  computer. 

This  investigation  assumes  that  the  guidance  equations  will  be 
solved  in  inertial  coordinates;  no  examination  is  made  of  the  possibility 
of  doing  all  or  part  of  the  computation  in  body-fixed  coordinates.     Further, 
no  consideration  is  taken  of  the  feedback  through  the  guidance  computer 
and  vehicle  autopilot. 

There  are  no  studies  available  which  compare  the  attractiveness 
of  various  coordinate  frames  available  for  computation  in  gimballess 
guidance.     For  this  reason,   the  present  investigation  of  the  performance 
and  design  of  a  system  which  provides  inertially  stabilized  information 
from  a  gimballess  inertial  package  is  a  useful  one.     Further,   since  the 
action  of  the  inertial  package  and  transformation  computer  is  much 
faster  than  the  feedback  present,   the  equipment  may  profitably  be  examined 
as  an  open  loop  device. 
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Since  the  inertia!  package  consists  of  three  delta-modulated  single- 
axis  integrating  accelerometers  and  three  delta-modulated  single-degree  - 
of -freedom  integrating  gyros,   the  theory  of  delta  modulation  and  its 
application  to  inertial  instruments  will  be  considered  in  some  detail. 
Methods  for  determining  the  zero-input  oehavior  of  deita-moduiated  loops 
will  be  reviewed,  and  the  response  of  such  ±oops  to  impulses  and  steps 
will  be  found.     Criteria  for  the  choice  of  loop  parameters  will  be  set 
forth. 

The  available  methods  of  transforming  a  vector  from  one  Cartesian 
coordinate  frame  to  another  will  be  examined.     There  are  essentially 
three  ways  of  viewing  the  transformation:  the  matrix  of  direction  cosines, 
the  four  parameter  representations,   and  the  vector  representation.     The 
properties  of  each  will  be  considered  and  the  differential  equations 
governing  the  transformation  parameters  will  be  derived.     Finally,   the 
transformations  will  be  compared  to  determine  their  suitability  for 
machine  computation,    considering  the  form  of  instrument  output. 

The  transformation  computer  will  be  designed  to  instrument  the 
transformation  found  to  be  best  for  machine  computation.     The  possibility 
of  accumulating  a  numDer  of  increments  from  the  gyros  and  instrumenting 
more  complicated  algorithms  to  permit    the  use  of  a  computer  with  a 
slower  clock  rate    and  reduction  of  error  will  be  examined.     With  the 
design  of  the  computer  fixed,   a  detailed  analysis  will  be  carried  out 
to  determine  the  dominant  mode  of  the  computer  and  estimates  of 
truncation  and  quantization  errors. 

Finally  a  few  typical  values  will  be  used  to  illustrate  the  method 
of  determining  the  basic  system  parameters  and  at  the  same  time  to 
provide  a  feel  for  the  order  of  magnitude  of  the  numbers  involved. 
.  D    Previous  Results 

Delta  modulation  was  originated  as  a  transmission  system  (Refs. 
2,    5,    6,    14)  with  little  attention  being  paid  to  the  feedback  aspect  of 
the  transmission  loop.     Later,    with  interest  developing  in  contro_ 
systems  using  reiays,    control  theory  was  applied  to  deita-moduiated 
loops  (Refs.   1,    3,    4,   6-13,   15-19)  under  the  names  bang-bang,   on-off, 
and  contactor  control  systems.     The  zero-input  behavior  of  these  loops 
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received  most  of  the  attention,  with  phase  plane  methods  playing  an 
important  part  (Refs.   8-11).     In  1950  a  method  for  determining  the  exact 
zero  input  behavior  of  a  loop  of  any  order  by  means  of  difference  equations 
was  presented  by  Torng  and  Meserve  (Ref.   19)  and  the  z -transform 
equivalent  of  this  method  was  suggested  by  Bergen  in  his  discussion  of 
their  paper  (Ref.    1).     Schwartz  elaborated  Bergen's  method  and  applied 
it  to  the  determination  of  the  loop's  response  to  step  inputs  (Ref.    15). 
Beyond  this,   no  work  has  been  done  previously  in  determining  the  response 
of  the  loop  to  step  inputs . 

The  application  of  delta  modulation  to  inertial  instruments  first 
took  place  in  1957,    and  development  work  has  continued  so  that  these 
instruments  are  now  a  very  attractive  choice  for  any  inertial  system 
which  requires  instrument  outputs  in  digital  form  (Refs.    4,    7,   8,    17,    18). 

Coordinate  transformations  have  received  only  cursory  examination 
in  connection  with  gimballess  inertial  reference  equipment.     Although 
some  reports  consider  the  possibility  of  using  transformations  other 
than  direction  cosines,   they  do  not  examine  them  in  any  detail.     Of 
course,   the  literature  on  coordinate  transformations   is  available  in 
almost  any  text  or  reference  on  classical  mechanics  or  mathematical 
physics  (Refs.   30-39). 

The  analysis  of  numerical  processes  has  a  fairly  long  history, 
and  most  techniques  are  easily  applied  to  digital  computers  (Refs. 
41,    44,   46,    48,    51-53).     However  the  application  of  frequency  domain 
methods  to  the  analysis  of  computers  is  quite  rare  (Refs.   49-51,    55,    61), 
and  thus  most  analysis  to  date  is  quite  cumbersome.     No  attempt  has 
yet  been  made  to  apply  frequency  domain  analysis  to  digital  computers 
in  which  the  vairable  of  integration  has  not  been  time. 

A.  number  of  private  companies  have  worked  in  the  field  of 
gimballess  guidance,  both  on  their  own  initiative  and  under  contract. 
Their  studies  have  mostly  provided  a  brief  survey  of  the  possibilities 
in  this  field,   but  no  detailed  analyses.     Very  little  of  this  work  is 
available  in  the  open  literature  (Ref.   66). 
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CHAPTER  II 


DELTA -MODULATED  LOOPS 


2.  A    Introduction 

Delta  modulation  is  a  puise  modulation  system  patented  in  France 
in  1946  (Ref.   6)  and  first  described  in  the  literature  by  de    Jager    (Ref.    5) 
and  others  (Ref.   14)  in  1952  as  a  system  for  speech  transmission.     In 
1957,   a  similar  system  was  developed  giving  digital  outputs  directly 
from  instruments  and  interest  in  such  instruments  has  remained  active. 
(Refs.    4,    7,    8,    12,    17,    18) 
2.B    Description  of  the  Loop 

The  input  is  translated  into  a  sequence  of  pulses,   each  having  the 
same  value,  but  with  either  polarity  being  possible.     The  pulses  are 
obtained  by  means  of  the  negative -feedback  loop  shown  in  Fig.    2-1.     One 
or  the  other  of  the  iinear  plant  functions  may  be  unity.     When  the  loop  is 
used  for  transmission,   the  linear  plant  in  the  forward  path  is  usually 
unity,   while  in  instrumentation  applications  the  feedback  plant  usually 
contains  a  zero-order  hold.     In  any  case,   an  integration  occurs  somewhere 
in  the  linear  plant  to  assure  zero  average  error  (Ref.   3,   p.   53).     The 
relay  may  have  either  two  levels  or  three  levels.     The  mathematical 
model  for  the  sampler  is  an  impulse  modulator  synchronized  with  a  clock. 

Each  impulse  appearing  at  the  output  represents  a  unit  change  in  the 
input  or  one  of  its  integrals  or  derivatives.     This  representation  of  a 
signal  in  terms  of  increments  leads  to  the  name  "delta  modulation".     If 
the  integration  appears  in  the  forward  path,   each  output  pulse  represents 
a  change  in  the  integral  of  the  input.     If  the  integration  appears  in  the 
feedback  path,   each  pulse  represents  a  change  in  the  input  itself.     Two 
integrations  in  the  feedback  path  cause  each  puise  to  represent  a  change 


(a) 


Two  Level  Relay 


(b) 

Three  Level  Relay 
Fig.   2-1    Delta-Modulated  Loops 
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in  the  slope  of  the  input. 

2.C    Analysis  for  Zero  Input 

With  no  input,   a  loop  with  a  two-level  relay  will  experience  a  self- 
excited  oscillation  or  limit  cycle,    and  if  the  zero  level  is  narrow  enough, 
a  loop  with  a  three-level  relay  may  also  have  a  limit  cycle.     From  physical 
considerations,   it  is  clear  that  the  periods  of  these  limit  cycles  are 
integral  multiples  of  the  sampling. period.     Further,   if  the  linear  plant 
contains  an  integration,   the  number  of  positive  pulses  in  a  cycle  will 
equal  the  number  of  negative  pulses  (Ref.   3,   p.   53). 
2.C.1    Moding 

Each  limit  cycle  exhibited  by  delta-modulated  loops  produces  a 
characteristic  sequence  of  output  pulses.     These  sequences  are  referred 
to  as  the  modes  of  the  system  and  are  very  important  in  the  analysis  of 
input-output  relationships.     The  modes  of  a  system  are  designated  by 
the  output  sequence  they  produce.     For  example,    if  the  limit  cycle  in  a 
loop  with  a  two-level  relay  produces  a  single  positive  pulse  followed 
by  a  single  negative  pulse,   the  mode  designation  is  1-1.     Similarly,    a 
mode  producing  two  positive  pulses  followed  by  two  negative  pulses  is 
designated  2-2.     If  the  loop  contains  a  three -level  relay  and  has  a  limit 
cycle  producing  two  positive  pulses  followed  by  two  negative  pulses  and 
one  zero  pulse,   the  mode  is  termed  2-0-2-1.     The  zero  for  the  second 
digit  indicates  that  no  zero  pulse  occurs  between  the  positive  and  negative 
pulses.     Note  that  the  mode  designations  always  begin  with  the  positive 
pulses. 

The  concept  of  moding  is  applicable  only  to  systems  with  no  input, 
since  the  limit  cycle  is  disturbed  by  inputs.     The  description  of  outputs 
in  response  to  inputs  is  discussed  in  Section  2.  D. 
2.  C.  2    Describing  Function  Analysis 

In  his  famous  paper(Ref .    58)  Kochenburger  presented  a  method  of 
steady-state  sinusoidal  analysis  applicable  to  non-linear  systems  in 
which  non-linear  elements  are  represented  as  linear  gains.     These 
gains  are  determined  as  functions  of  the  amplitude  and  frequency  of 
the  inputs  to  the  non-linear  element.     Chow  (Ref.   3)  and  Russell  (Ref.   13) 
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first  used  the  describing  function  to  determine  limit  cycles  in  delta- 
modulated  loops.     Their  interest  was  in  adjusting  the  width  of  the  zero 
level  in  a  three-level  relay  to  eliminate  limit  cycles,  but  their  analysis 
is  applicable  to  any  sampled-data  loop  containing  any  non-linearity. 

Their  method  of  analysis  is  the  same  as  the  standard  describing 
function  analysis  applied  to  continuous  time  systems  with  the  addition 
of  a  variable  delay  due  to  sampling.    This  variable  delay,   6,   arises  because 
the  output  pulses  can  occur  only  at  sampling  instants,   even  though  the 
relay  had  switched  earlier.     This  delay,   which  ranges  from  zero  to  a 
full  sampling  period,    T   ,    is  easily  understood  from  an  examination  of 
Fig.   2-2.     The  marks  on  the  time  axis  indicate  clock  pulses.     In  each 


Fig.   2-2  Time  Delay  due  to  Sampling 

case,   the  negative  output  pulse  occurs  with  the  second  clock  pulse  shown 

and  6  ranges  from  almost  ah  entire  sampling  period  in  Fig.   2 -2(a)  to 

practically  zero  in  Fig.   2 -2(c).     The  maximum  phase  shift  due  to  sampling, 

<j>   ,   is  the  product  of  frequency  and  sampling  period,   and  is  given  by 
s 

T 

</>s  =  u/Ts  =  2tt     ^§  =  2iF~  (2-1) 

s 

where  T  is  the  period  corresponding  to  the  input  angular  frequency  w. 

Standard  describing  function  analysis  (Ref.   63,   Ch.   10,   pp.    559-580) 

is  carried  out  by  graphical  solution  of  the  equation 

1  +  G(jw)N(A,w)  =  0  (2-2) 

G(joj), the  transfer  function  for  the  linear  part  of  the  system,   is  plotted 

1 
on  a  Nyquist  diagram  with  -   sjs T — :,   the  negative  of  the  reciprocal  of 

the  describing  function  for  the  non-linear  element.     Intersection  of  the 
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two  curves  indicates  a  root  of  Eq.    (2-2)  and  predicts  a  limit  cycle.     In 
the  analysis  of  delta-modulated  loops  the  phase  shift  due  to  sampling 
must  be  added  to  the  plot;  this  is  most  conveniently  done  by  adding  <j> 

o 

as  a  phase  lag  at  frequencies  on  the  G(jco)  locus  for  which  a  limit  cycle 
is  possible,   that  is,    at  submultipies  of  the  sampling  frequency.     In  this 
application,   a  log  polar  plot  is  most  convenient  for  making  the  Nyquist 
plot.     It  is  termed  a  "moding  plot".  , 

In  Fig.   2-3,    a  moding  plot  for  G(jco)  =     .    , .  -,.    is  shown.     The 

negative  reciprocal  describing  function  for  a  two-level  relay  covering 
the  whole  negative  real  axis  is  shown.  As  can  be  seen  from  this  plot, 
limit  cycles  are  predicted  for  angular  frequencies  of  u  /2,  wa/4,  and 
co  /6.    At  co  /8  the  sampling  phase  is  at  most  ir  /4,   which  is  not  sufficient 

S  S 

to  reach  the  negative  real  axis.     It  is  clear  that  limit  cycies  will  not 
exist  at  lower  frequencies;  thus  the  modes  predicted  are  1-1,    2-2,   and 
3-3.     For  a  three-level  relay,   the  describing  function  reaches  a  maximum 
when  the  amplitude  of  the  relay  input  is  n/T  times  the  zero  level  half-width, 
A   .     Thus  the  reciprocal  of  the  describing  function  has  a  minimum,   and 
the  zero  level  half- width  may  be  adjusted  so  that  no  limit  cycle  is 
predicted.     An  example  of  this  calculation  is  shown  in  Appendix  B. 
2.C.3    Detailed  Analysis 

The  describing  function  is  only  an  approximate  tool,   albeit  a  powerful 
one.     For  this  reason  attempts  were  made  to  achieve  exact  answers. 
The  step-by-step  transient  method  was  discarded  as  being  too  cumbersome. 
Phase-plane  methods  received  more  attention,    and  although  limited  to 
second  order  systems  became  quite  elegant  (Refs.   9-11).     Torng  and 
Meserve  (Ref.   19)  presented  a  method  for  determining  the  exact  output 
of  the  linear  plant  of  a  delta-modulated  loop  at  sampling  instants  using 
discrete  orthogonal  functions.     Bergen  (Ref.   1),   in  his  discussion  of 
their  paper,   outlined  a  similar  method  using  z -transform  techniques. 
His  method  allows  the  possibility  of  determining  no-transient  initial 
conditions  for  the  limit  cycle.     By  manipulation,   the  state  of  the  iinear 
plant  at  each  sampling  instant  may  also  be  determined. 

Each  method  assumes  a  mode  at  the  relay  and  applies  it  as  an  input 
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Fig.   2-3    Moding  Plot 
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to  the  linear  plant.     The  output  of  the  linear  plant  resulting  from  this 
input  is  examined  at  the  sampling  instants  to  see  if  the  signal  has  the 
proper  polarity  to  sustain  the  postulated  mode.     Thus  if  the  1-1  mode 
be  possible,   the  linear  plant  output  at  sampling  instants  must  be  a 
repeating  sequence  of  two  numbers:  the  first  positive,   the  second  negative. 
Torng  and  Meserve  do  this  using  discrete  orthogonal  functions  in  the  time 
domain.     Bergen  operates  in  the  frequency  domain  with  z-transforms,    and 
introduces  initial  conditions  in  this  way. 

The  modes  which  must  be  investigated  are  limited,    and  may  be 
determined  by  a  describing  function  analysis  as  carried  out  above  (Sec  » 
2.  C.  2).     Torng  and  Meserve  list  the  modes  which  must  be  investigated 
for  limit  cycles  with  periods  of  two,   four,   six,    and  eight  sampling  periods 
(Ref.   19,    Table  II,   p.    303). 

For  a  loop  with  a  three-level  relay,    choice  of  the  zero-level  width 
so  that  a  describing  function  analysis  predicts  no  limit  cycle  may  be 
made.     However,   for  fast  systems  (T  fi  >  1),   this  choice  of  zero-level 
width  may  not  eliminate  low  energy  limit  cycles.     To  eliminate  these, 
a  transient  analysis  is  the  easiest.     The  linear  plant  output  is  assumed 
to  be  just  at  the  edge  of  the  zero  level  of  the  relay.     Thus  one  feedback 
pulse  drives  the  linear  plant  back  toward  null.     For  no  limit  cycle  to 
occur,   the  limit  of  the  linear  plant  motion  must  not  reach  the  outer 
edge  of  the  zero  level,   even  for  very  long  times.     The  result  of  this 
analysis  for  any  linear  plant  containing  an  integration  is  that  the  zero 
level  width  should  be  chosen  greater  than  the  product  of  the  sampling 
period  and  the  steady  state  velocity  of  the  plant  under  the  action  of  the 
relay  alone. 
2.C.4    Compensation 

As  Schwartz  (Ref.    15)  and  Sliwkowski  (Ref.   16)  point  out,   in  limit 
cycling  loops     inputs  may  not  be  sensed  until  one  full  period  of  the 
prevailing  mode  has  elapsed.     Another  view  of  this  problem  is  that  in 
limit  cycling,   the  loop  is  storing  information;  the  amplitude  of  the  limit 
cycle,   indicated  by  the  moding,    is  a  measure  of  the  information  storage. 
This  information  storage  should  be  kept  to  a  minimum  if  the  system 
is   to  be  used  to  best  advantage.     In  a  loop  with  a  three-level  relay  the 
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width  of  the  zero  level  is  a  measure  of  stored  information,   and  for  slow 

systems  (T   / 1  <  .1)  may  become  quite  large.     For  this  reason  it  may  be 
s 

desirable  to  compensate  the  system  so  that  all  but  the  1-1  mode  are 

eliminated,   or  so  that  the  zero  level  of  a  three -level  relay  is  reduced 

to  a  reasonable  size.     There  are  three  methods  of  compensation  which 

suggest  themselves:  dc  analog,    ac  analog,   and  discrete  or  digital. 

Compensation  is  best  viewed  in  terms  of  the  moding  plot  (Fig.   2-3). 

To  eliminate  the  2-2  mode,   more  than  90     phase  margin  must  be  provided 

for  G(jco)  at  w  /4.     Similarly  to  eliminate  the  3-3  mode  60     of  phase 
s 

margin  must  be  provided  at  oo  /6  and  so  on.     With  these  conditions  met, 

s 

the  extended  G(joo)  locus  will  not  intersect  the  describing  function  except 
at  the  1-1  mode.     For  a  three -level  relay,   the  zero  level  width  may  be 
adjusted  so  that  the  reciprocal  describing  function  reaches  down  to  just 
above  the  1-1  level. 

Dc  analog  compensation  methods  are  covered  in  detail  in  most 
elementary    automatic   control  textbooks  (Ref.    57,    Vol.   I,   Ch.   12;  Vol  II, 
Ch.   4;  Ref.   64,  Ch.    6;  for  example).     Adjusting  the  moding  plot  to  achieve 
the  objectives  listed  above  is  a  straightforward  task,     Ac  analog  compensation 
is  also  treated  in  many  texts  (Ref.    63  Ch.   6;  Ref.   64  Ch.   6;  for  example). 
The  desirability  of  this  form  of  compensation  stems  from  the  fact  that 
the  relay  is  often  a  phase -sensitive  switch  operating  on  information  in 
a  supressed-carrier  wave;  thus,   if  dc  compensation  were  used,   the 
signal  would  have  to  be  demodulated.     However,   the  signal  frequencies 
used  in  practice  are  very  close  to  the  carrier  frequencies,   and  de- 
modulation is  very  difficult  under  these  conditions.     For  this  reason,    ac 
compensation  is  quite  attractive.     It  should  be  noted  that  the  use  of  low- 
pass  to  band-pass  transformations  do  not  provide  sufficiently  accurate 
compensation  for  this  application. 

Discrete  or  digital  compensation  has  received  little  attention  in  the 
literature.     The  advantage  of  this  form  of  compensation  is  that  there  is 
no  problem  of  drift  as  compared  to  analog  compensation,    and  non-linear 
or  logical  compensation  is  possible.     Most  books  on  sampled -data  systems 
include  sections  on  the  techniques  of  linear  discrete  compensation  and 
Tou  (Ref.   62,   pp.   444  ff)  lists  three  methods  for  realizing  linear  discrete 
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transfer  functions:  computers,   delay  lines,   and  pulsed  RC  circuits. 

The  point  in  the  delta-modulated  loop  at  which  the  signal  is  discrete 
is  after  the  sampler;  it  is  here  that  the  discrete  compensator  must  be 
inserted.     If  a  pulsed  RC  circuit  is  used,   a  second  sampler  must  follow 
it. 

Linear  discrete  compensation  is  best  viewed  in  the  G(z)  plane.     The 
locus  of  the  linear  plant  has  extensions  similar  to  those  in  the  moding 
plot,  because  of  sampling  delays.     The  compensation  must  be  chosen  so 
that  the  extended  G  (s)  locus  does  not  intersect  the  negative  reciprocal 
describing  function  locus. 

Logical  discrete  compensation,   using  branching  operations  which 
influence  the  relay  on  the  basis  of  past  output  pulses,  has  been  tried  only 
once  (Ref.   20)  with  little  success.     It  would  seem  that  there  are  great 
possibilities  for  this  type  of  compensation. 
2.D    Input -Output  Relationships 

As  explained  above,   each  output  pulse  represents  a  unit  change  in 
the  input  or  one  of  its  integrals  or  derivatives.     Thus  with  a  finite  input, 
the  algebraic  sum  of  the  output  pulses  should  average  to  the  input  or  one 
of  its  integrals  or  derivatives  or,   in  the  case  of  zero  input,  to  zero. 
The  system  overloads  when  the  rate  of  change  of  the  input  (or  one  of 
its  integrals  or  derivatives)  exceeds  that  value  which  requires  all  the 
output  pulses  to  be  of  one  sign  (Ref.   2,  p.   145).     Therefore  a  convenient 
way  of  expressing  the  input  is  as  a  fraction  of  the  maximum  permissible 
signal. 

Since  the  position  and  number  of  integrators  in  the  loop  change  input - 
output  relationships,   a  general  analysis  is  very  cumbersome.     For  this 
reason,  the  analysis  discussed  here  considers  a  loop  with  one  integrator 
in  it.     The  linear  plant  in  the  feedback  path  is  moved  to  the  forward  path, 
with  the  exception  of  zero-order  hold,  if  one  is  present,   which  is  left 
in  the  feedback  path.     This  requires  the  reciprocal  of  the  moved  plant 
to  be  placed  in  front  of  the  loop.     Fig.   2-4  shows  the  rearranged  loop. 
The  dotted  lines  around  the  hold  indicate  that  it  may  or  may  not  be 
present.   GMs)  is  the  part  of  the  linear  plant  moved  out  of  the  feedback 
path.     The  analysis  which  follows  ignores  the  1/G'  (s)  in  front  of  the 
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loop  since  it  tends  to  confuse  the  issue. 
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Fig.   2-4  Delta-Modulated  Loop,   Rearranged 

2.D.1    Transient  Response 

The  following  analyses  are  exact  determinations  of  the  inpulse  and 
step  responses  of  delta  modulated  loops.     With  the  exception  of  Schwartz's 
work  (Ref.   15)  on  step  inputs  for  loops  with  two-level  relays,   they  have 
not  been  previously  described. 
2.D.  l.a    Impulse  Response 

The  height  of  an  impulse  (infinity)  exceeds  the  capacity  of  the  system; 
the  best  output  makes  up  the  integral  of  the  impulse  a  step  at  a  time  to 
within  the  quantization  error  of  the  loop.     Clearly  if  the  value  of  the 
impulse  is  too  low,   no  output  will  appear,   except  for  the  limit  cycle. 

The  analysis  of  the  loop's  impulse  response  is  best  done  by  a  step- 
by-step  transient  method.     The  applied  impulse  drives  the  linear  plant 
some  distance  from  null;  calculation  of  this  excursion  in  the  presence 
of  the  feedback  signal  is  a  simple  task.     The  linear  plant  output  at 
sampling  instants  is  examined  for  the  time  at  which  the  feedback  signal 
changes  sign.     If  the  loop  has  a  two-level  relay,   the  state  of  the  linear 
plant  is  compared  with  the  set  of  initial  conditions  corresponding  to 
each  mode  to  determine  into  which  mode  the  system  will  settle.     Usually 
this  will  be  the  highest  mode.     The  output  of  the  linear  plant  plant  in  a 
loop  containing  a  three-level  relay  is  examined  to  see  if  it  will  cause 
the  relay  to  have  an  output  of  opposite  sign  from  the  input.     If  not,   the 
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computation  is  complete;  if  so,   it  must  be  carried  further. 
2.D.1.D    Step  Response 

A  step  input  of  magnitude  D,   where  D  is  the  drive  level  of  the  relay, 
will  saturate  the  system.     For  this  reason  the  magnitudes  of  step  inputs 
are  given  as  a  fraction,    r,  of  D. 

For  a  periodic  sequence  to  appear  at  the  output,   the  input  must  be 
expressed  as  a  rational  fraction  of  the  saturating  input.     This  is  not  a 
serious  limitation,   since  any  number  less  than  one  may  be  well  approximated 
by  a  rational  number  with  a  small  denominator 

Suppose  the  input  is  a  rational  fraction,   say  U/U,   of  the  saturating 

signal.     Then  the  output  will  be  a  sequence  of  NW  pulses.     For  a  loop 

with  a  two-level  relay,   there  will  be  N     =  (W+U)/2  positive  pulses  and 

N     =  N(W-U)/2  negative  pulses.     This  may  be  shown  by  a  short  computation. 

The  relationship  between  U,    W,   N  ,   and  N    is 

N+-N~    =    U 

N++N"         W  (2-3) 

Rearrangement  of  this  ratio  yields 

N+    =    W+U 

N  W-U  (2-4) 

A  loop  with  a  three-level  relay  may  have  N     zero  pulses.     In  this  case, 

the  required  relation  between  U,   W,   N  ,   N   ,   and  N    is 

n+-n"         u 


N++N"+N°        W  (2-5) 

If  U  is  positive,   and  there  are  no  negative  pulses,  then 

N+  U 


N°  W-U  (2-6) 

With  ratios,  of  positive,  negative,   and  zero  pulses  determined,   it  only 
remains  to  determine  the  organization  of  these  pulses  within  the  sequence. 
This  may  be  done  by  determining  bounds  on  the  number  of  consecutive 
pulses  of  one  sign  under  various  conditions.     Different  bounds  are  of 
interest  depending  on  whether  the  loop  has  a  two-level  or  a  three -level 
relay. 
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For  a  loop  with  a  two -level  relay,   the  bounds  that  are  important  are 
the  maximum  number  of  pulses  of  one  sign  that  are  possible  and  the 
minimum  number  of  pulses  of  one  sign  that  are  possible  following  a 
given  number  of  pulses  of  the  opposite  sign. 

The  upper  bound  is  found  by  observing  the  greatest  possible  excursion 
from  null  and  counting  the  number  of  clock  periods  required  to  drive  the 
linear  plant  back  across  null.     The  maximum  excursion  from  null  occurs 
when  the  linear  plant  output  travels  away  from  null  at  maximum  rate  for 
a  full  sampling  period. 

The  lower  bound  is  found  by  assuming  minimum  initial  conditions. 
The  minimum  displacement  is  clearly  zero  plus,   and  the  minimum 
derivatives  result  from  maximum  velocity  at  the  previous  switching 
time  a  given  number  of  sampling  periods  earlier.     Under  these  initial 
conditions,   the  number  of  sampling  periods  required  for  the  system  to 
cross  null  is  found. 

The  bounds  to  be  determined  in  the  case  of  a  loop  with  a  three -level 
relay  are  the  maximum  number  of  pulses  of  the  same  sign  as  the  input 
and  the  maximum  number  of  opposite  sign. 

Determination  of  the  maximum  number  of  pulses  of  the  same  sign 
as  the  input  proceeds  from  the  assumption  that  the  linear  plant  has 
travelled  at  maximum  rate  under  the  influence  of  the  input  alone  for  a 
full  sampling  period.     At  this  point,   positive  pulses  begin  to  drive  the 
plant  back  to  null.     The  pulses  required  to  bring  the  linear  plant  output 
within  the  zero  level  are  counted. 

The  possibility  of  pulses  of  a  polarity  opposite  to  the  input  may  be 
examined  by  a  similar  analysis.  It  is  assumed  that  the  linear  plant  is 
driven  across  the  zero  level  under  the  influence  of  the  feedback  signal 
and  input  for  a  full  sampling  period  starting  from  just  outside  the  zero 
level.  The  feedback  signal  switches  off,  and  the  system  output  is 
examined  at  successive  sampling  times  to  see  whether  it  gets  outside 
the  zero  level  on  the  other  side. 

Examples  of  the  above  computations  are  shown  in  Appendices  A  and  B. 

With  these  bounds  set,   the  sequences  having  the  required  average 
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may  be  written  down.     These  sequences  are  designated  in  a  way  similar 
to  that  used  for  modes.     For  example,   a  sequence  from  a  loop  with  a 
two -level  relay  consisting  of  three  positive  pulses  followed  by  a  negative 
pulse  is  designated  3,   1;  one  having  four  positive  pulses,   one  negative 
pulse,  two  positive  pulses,   and  another  negative  pulse  is  designated 
4,  1,   2,   1.     A  loop  with  a  three -level  relay  having  an  output  sequence  of 
a  positive  pulse  followed  by  a  zero  pulse  would  have  that  sequence  de- 
signated 1,  1,   0,   0.     Nothe  that  each  of  these  sequences  has  an  average 
of  one -half.     The  permissible  sequences,  which  are  few  because  of  the 
elimination  process  just  described,  may  be  tested  by  Bergen's  method 
to  see  if  they  can  actually  occur.     The  sequence,     with  the  input,  is  assumed 
to  drive  the  linear  plant,   and  the  signs  of  the  linear  plant  output  are 
examined  at  sampling  times  to  see  if  the  postulated  sequence  is  sustained. 

With  inputs  smaller  than  about  one -tenth  of  the  saturating  signal, 
the  period  of  the  output  sequence  becomes  rather  too  long  for  the  analysis 
above  to  handle  conveniently,   and  a  modified  step -by -step  transient 
procedure  is  preferable. 

In  a  loop  with  a  two -level  relay  and  a  small  input,  the  prevailing  mode 
will  continue  for  a  number  of  cycles.     The  feedback  signal  has  zero  average 
over  a  full  period  of  the  mode,   and  its  effect  may  be  ignored.     This  is 
perhaps  best  visualized  in  phase  space  with  the  mode  trajectory  moving 
slowly  through  the  space  under  the  influence  of  the  input.     When  the 
position  coordinate  of  one  of  the  sampling  points  of  the  trajectory  moves 
to  a  new  region  of  phase  space  (from  positive  to  negative,  for  example), 
the  prevailing  mode  will  change.     The  mode  changes  either  by  dropping 
a  pulse  or  by  adding  a  pulse.     In  any  case,  the  total  state  of  the  linear 
plant  after  the  mode  changes  may  be  found  and  compared  with  the  initial 
states  of  the  other  modes.     Thus  the  new  mode  of  the  system  is  found. 
The  trajectory  of  the  new  mode  is  assumed  to  move  under  the  influence  of 
the  input,  with  a  new  initial  position  and  an  initial  velocity  which  is  the 
velocity  of  the  plant  relative  to  the  new  prevailing  mode.     This  process 
is  continued  until  the  entire  sequence  is  determined. 

For  the  loop  with  a  three -level  relay,  the  appropriate  calculation 
is  a  step-by-step  transient  response.     Assume  an  initial  condition,   and 
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follow  the  plant  through  until  the  sequence  is  determined.     In  almost 
all  cases  only  one  positive  output  pulse  will  occur  between  zero  pulses. 
2.  D.  1.  c    Effects  of  Imperfect  Integration 

Thus  far,   it  has  been  assumed  that  the  integrator  has  been  perfect. 
In  some  cases,   the  integrator  is  in  reality  a  first  order  lag  with  a  long 
time  constant.     This  may  be  viewed  as  a  small  elastic  restraint.     With 
an  elastic  restraint  present,  there  is  a  limit  below  which  inputs  cannot 
be  detected.     This  is  also  best  seen  in  terms  of  the  mode  trajectory  in 
phase  space.     If  the  steady-state  response  of  the  linear  plant  to  a  small 
input  is  insufficient  to  move  the  trajectory  far  enough  to  cause  the  mode 
to  change,   no  input  will  be  detected.     Of  course,   with  a  three -level  relay, 
this  critical  value  of  linear  plant  output  is  the  edge  of  the  zero  level. 
2.D.l.d    Ramps  and  Higher  Order  Inputs 

Ramps  and  higher  order  inputs  (e.g.   parabolic)  cannot  be  followed 
by  the  loop  over  infinite  time  since  the  amplitude  of  such  a  signal  increases 
without  limit.     A  sawtooth  wave  may  be  analyzed  in  terms  of  the  frequency 
response  described  below.     For  ramps  which  rise  to  a  certain  value  and 
then  remain  constant,   a  step-by-step  transient  computation  seems  to  be 
the  most  convenient.     The  bounds  on  number  of  pulses  are  useful  in  this 
computation. 
2.D.  2    Frequency  Response 

Since  the  frequency  response  method  of  servomechanism  design  is 
familiar  to  most  engineers,   the  following  modification  of  the  describing 
function  is  presented. 

A  sinusoidal  signal  is  assumed  at  the  relay.     The  standard  describing 
function  is  used  to  represent  the  relay,   assuming  enough  low  pass  filtering 
in  the  loop  to  make  the  describing  function  valid.     The  fundamental  of  the 
relay  is  passed  through  the  sampler  and  hold,   giving  the  feedback  signal. 
With  the  input  to  the  relay,   which  is  the  output  of  the  linear  plant,  known, 
the  error  signal  can  be  determined.     The  input  to  the  loop  is  the  sum  of 
the  error  and  feedback  signals.     With  the  relay  output  taken  as  the  loop 
output,   the  describing  function  for  the  entire  loop  is  determined  as  the 
ratio  of  relay  output  to  loop  input. 

Three  comments  must  be  made  about  the  above  analysis.     First, 
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the  output  of  the  loop  has  been  taken  as  the  relay  output,   with  no  limit 
cycle  present.     Thus,   this  analysis,   for  frequencies  approaching  the 
limit  cycle  frequency,    is  open  to  question.     Second,   there  is  a  minimun 
amplitude  of  signal  which  can  be  detected, corresponding  either  to  the 
limit  cycle  amplitude  or  the  zero-level  width  of  the  relay.     Third,   the 
bandwidth  of  the  loop  is  limited  by  the  sampling  frequency,  and  signals 
above  this  frequency  will  not  be  passed;  instead  they  appear  as  steps. 
(Ref.  19). 
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CHAPTER  III 
DELTA- MODULATED  INSTRUMENTS 


3.  A    General  Comments 

The  delta-modulated  loops  discussed  in  Chapter  II  are  to  be  applied 
to  two  inertial  instruments,  the  gyroscope  and  the  accelerometer. 
The  gyro  to  be  used  is  the  single -degree-of-freedom  floated  integrating 
gyro  described  by  Draper,    Wrigley  and  Grohe  (Ref.   21).     The 
accelerometer  is  a  single-axis  floated  integrating  pendulum  type 
described  by  Hutchings  and  others  (Refs.   4,   7,   8,   12,   17). 
3.A.1    Brief  Description  of  the  Instruments 

The  integrating  gyro,   when  behaving  in  an  ideal  manner,  is 
described  by  the  following  differential  equation: 

JAoA  +  C^OA  =  HspWIA  +  Mtg  (3-la) 

where  J  is  the  moment  of  inertia  of  the  gyro  float  about  its  output 

axis; 

C  is  the  damping  coefficient  of  the  fluid  damper; 

H       is  the  spin  angular  momentum  of  the  gyro; 
sp 

WTA  is  the  angular  velocity  of  the  gyro  case  about  its  input  axis, 

with  respect  to  inertial  space; 
M.     is  the  torque  applied  to  the  instrument  float  by  the  torque 

generator; 
and  A        is  the  angle  through  which  the  float  has  turned  with 

respect  to  the  case  avout  the  instrument  output  axis. 
A  dot  over  a  quantity  (e.g.  A^.)  denotes  time  differentiation. 
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Eq.   (3 -la)  may  be  rewritten  as 

H  M 

7AOA  +  AOA=    ^WIA+  -T 


(3-lb) 


where  7  =  J/C  is  the  characteristic  time  or  time  constant  of  the  gyro. 
The  differential  equation  describing  the  accelerometer  is 

jiOA  +  C^OA  =  PaIA  +  Mtg  (3-2a) 

where  the  symbols  above  are  similarly  defined;  where  P  is  the 
pendulosity  of  the  accelerometer;  and  aT»     is  the  non- gravitational 
acceleration  of  the  instrument  with  respect  to  inertial  space  along 
its  input  axis. 
Eq.   (3 -2a)  may  be  rewritten  as 

P  Mt 

7AOA  +  AOA  =  C"  aIA  +  "TT  (3-2b) 

where  again  7  =  J/C. 

The  torque  generators  are  fed  from  a  constant  current  source,   with 
the  current  switched  between  two  opposing  windings  to  produce  torque 
of  different  signs,   or  switched  to  a  dummy  load  or  off  to  produce  zero 
torque.     The  required  sign  of  the  torque  is  determined  at  sampling 
times  and  is  either  opposite  to  the  sign  of  the  float  angle  or  zero, 
depending  on  the  logic  of  the  loop.     Therefore, the  torque  generator 
acts  as  a  zero-order  hold  in  the  feedback  path. 
3.  A.  2    Form  of  the  Loops 

The  delta-modulated  instruments  may  be  represented  as  shown 
in  Fig.   3-1.     These  loops  are  identical  in  form  to  each  other  and  to 
the  loops  analyzed  in  Appendices  A  and  B.     The  output  pulses  represent 
changes  in  the  integral  of  the  input,   so  that  the  gyro's  output  pulses 
indicate  changes  in  the  angular  orientation  of  the  gyro  case  about  its 
input  axis,   and  the  accelerometer's  output  pulses  indicated  changes 
in  velocity  along  the  accelerometer's  input  axis. 
3.B    Delta-Modulated  Gyros 
3.B.1    Detailed  Description 

A  detailed  physical  description  of  the  single -degree -of -freedom 
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floated  integrating  gyro  is  presented  in  Reference  21.     Current 
practice  is  to  center  the  gyro  float  in  the  case  with  a  magnetic 
suspension  system  similar  to  that  described  in  Reference  26,   as 
well  as  to  float  it.     Nominal  values  of  the  characteristics  of  some 
gyros  are  given  in  Reference  25.     Fig.    3-2  shows  a  schematic 
drawing  of  the  instrument. 

The  differential  equation  governing  the  gyro  may  be  derived 
from  Euler's  equations  of  motionfor  arigid  body  about  principal 
axies.   It  is  assumed  that  the  bearings  of  the  gyro  are  infinitely  stiff, 
so  that  only  motions  about  the  output  axis  need  be  considered.     The 
bearing  stiffness  assumption  is  quite  well  satisfied  for  gyros  even 
under  high  angular  velocities.     The  set  (IA,   OA,   SRA)  forms  a  right- 
handed  triad.     The  equation  is 


Mapp  =  "Hsp(WIA  cos  AOA  -  W3RA    sin  AOA>  +  J(AOA+  Wca> 
+(IIA  -  W^IA2  "  WSRA2  )  Sin  AOA    cos  AOA 

+(IIA"  ISA)(c0s2a0A  -  sin2  AOA)WIAWSRA  (3-3a) 

M;        =  M.      -  CA~A  (3-3b) 

app  tg  OA 

where  M         is  the  torque  applied  to  the  instrument  float  about 

its  output  axis; 
WnR  a  is  the  angular  velocity  of  the  gyro  float  about  its  spin 

reference  axis,   with  respect  to  inertial  space; 
W       is  the  angular  velocity  of  the  case  about  its  output  axis,   with 

respect  to  inertial  space; 
I,,   is  the  moment  of  inertia  of  the  gyro  float  about  its  input  axis; 

L.   is  the  moment  of  inertia  of  the  gyro  float  about  its  spin  axis. 

These  equations  may  be  rewritten  as 

lAOA  +  AOA  -  TP  WSRA    sin  AOA  ■  -C*  WIA  cos  AOA  + 
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Mtg ,„2 


~7Wca  +    ~&+  (ISA  -  W(W  IA  "  W  SRA)  sin  AOA  cos  AOA/C 

+  (ISA  '  IIA)  (cos2  AOA  -  sin'  AOA)WIAWSRA/C 

(3-3c) 

The  last  two  terms  on  the  right  side  of  Eq.    (3 -3c)  are  anisoinertia 

effects  and  are  negligible  in  practical  gyros.     They  will  not  be 

considered  further  here. 

Since  A~A  stays  well  below  a  degree  in  practice,   the  standard 

small  angle  approximations  may  be  made  in  Eq.    (3 -3c): 

H  H  .  M 

7AOA  +  AOA  -  -TP- WSRAAOA  =  "I?  WIA  "  7  Wca  +  "C* 

(3-3d) 

For  WCT3A   =  W       =0,   Eq.    (3 -3d)  is  identical  with  Eq.    (3 -lb). 

3.B.2    Rotational  Coupling 

Eq.   (3-3d)  shows  that  rotations  of  the  gyros  about  all  axes  affect 
its  motion.     The  coupling  effects  of  these  motions  will  be  considered 
in  detail. 
3.B.2.a    Non-Reversing  Rotations 

Examination  of  Eq.    (3 -3d)  shows  that  an  angular  velocity  about 
the  spin  reference  axis  appears  as  an  elastic  restraint  term  in  the 
gyro  equation.     This  elastic  restraint  causes  imperfect  integration  in 
the  loop.     As  noted  in  Chapter  II,    Section  2.D.I.  c,   imperfect  integration 
in  the  loop  produces  a  dead  zone,  and  inputs  within  this  dead  zone 
cannot  be  detected.     This  is  a  significant  phenomenon  in  gimballess 
inertial  reference  equipment,   since  large  angular  velocities  may 
occur  about  any  axis. 

Angular  acceleration  of  the  case  about  its  output  axis  produces 
a  torque  on  the  instrument  float  of  the  inertia,  reaction  category.     As 
can  be  seen  from  Eq.    (3 -3d),   this  torque  is  proportional  to  the  gyro 
time  constant  and  acts  as  a  spurious  input  to  the  gyro. 

These  torques  may  be  compensated  out  either  directly  in  the 
gyro  or  in  the  computer,   if  it  is  found  that  they  seriously  affect  the 
performance  of  the  gyro,   in  a  manner  similar  to  that  used  for 
compensation  of  compliance  effects  in  present  gyros. 
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3.B.2.b    Vibratory  Rotations 

Vibratory  rotations  of  the  gyro  about  one  axis  at  a  time  do  not 
affect  the  performance  of  the  instrument.     When  this  type  of  motion 
occurs  about  two  axes  at  a  time,   however,   it  can  produce  a  drift. 

Consider  sinusoidal  motion  of  the  case  about  the  spin  reference 
axis  and  the  input  axis.     With  sinusoidal  motion,  the  gyro  equation  of 
motion,  ignoring  torque  generator  torque,   is 

TAOA  +  AOA  -  -TT  WSRA  sin  <2*fSRAt)AOA  =^T  WIA  sin  »"**  +  *> 

(3-4) 
A  solution  for  this  equation  using  series  techniques,   which  coverges 
for  WgRA  Hg    /C  <  1,  has  been  given  by  Haff  and  Meltzer  (Ref.   27,   Ch.    5), 

A  further  discussion  of  the  series  solution  and  its  convergence  is  given 
in  Reference  29. 

Haff  and  Meltzer  show  that  if  nf„R.    =  mL.,   a  constant  drift  is 

predicted  by  the  series  solution  which  they  propose.     A  similar  effect 
occurs  if  a  vibration  about  the  output  axis  occurs  with  the  vibration  about 
the  spin  reference  axis.     Further,   if  fqRA  is  harmonically  related  to 

the  clock  frequency,   drift  may  also  be  present,   as  can  be  seen  by  writing 
the  Fourier  series  expansion  for  the  torque  generator  torque. 

Angular  vibrations  of  the  gyro  about  any  two  axes  produce  a  conical 
motion.     Computation  of  the  orientation  of  the  gyro  may  be  made  in  such 
a  way  as  to  introduce  a  constant  angular  velocity  about  the  cone  axis  in 
the  computed  orientation,   when  none  actually  exists.     Haff  and  Meltzer 
present  a  general  treatment  of  this  effect  for  a  stabilized  platform 
(Ref.   27,   Ch.   3,   Sec.   3.4).     It  is  clear  from  their  discussion  how  the 
computation  implied  by  a  stable  platform  leads  to  this  error.     This 
coning  error  does  not  exist  in  gimballess  systems. 
3.B.3    Linear  Motions 

Many  writers  (Ref.    22,  for  example)  have  treated  the  effects  of 
linear  motions  on  gyros  and  they  will  not  be  considered  here. 
3.B.4    Choice  of  Loop  Parameters 

While  the  parameters  of  the  loop  must  be  chosen  so  that  production 
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gyros  and  available  electronic  components  may  be  used,   a  certain 
amount  of  latitude  exists  in  choosing  some  parameters. 

The  spin  angular  momentum  of  the  gyro  is  chosen  on  the  basis 
of  the  allowable  drift  rate.     The  higher  the  drift  acceptable,  the 
lower  the  angular  momentum  may  be. 

The  torque  level,   M,   required  in  the  torquer  depends  on  the 
maximum  angular  velocity  that  the  system  must  follow,   W         .     The 

product  of  angular  momentum  and  W  gives  the  lower  bound  on  M. 

For  gyros  with  drift  characteristics  acceptable  in  inertial  guidance 
applications,  this  bound  may  be  quite  high,  and  torquers  other  than 
microsyns  may  be  required. 

The  clock  frequency  may  be  determined  from  the  maximum 
angular  velocity  expected  and  the  increment  size  in  the  loop,   A0,  . 
The  quotient  W         /A9  gives  the  lower  bound  on  clock  frequency. 

If  an  ac  signal  generator  is  used,  the  carrier  must  be  a  multiple  of 
the  clock  frequency.     This  requirement  follows  from  the  desirability  of 
peak  detection  in  the  switch.     If  the  carrier  frequency  is  not  a.multiple 
of  the  clock  frequency,  the  phasing  of  the  clock  and  the  carrier  will 
not  remain  constant,   and  the  clock  pulses  cannot  be  adjusted  to  appear 
at  the  peaks  of  the  carrier  sinusoid. 

Since  it  is  desirable  to  have  the  uncompensated  mode  of  the  loop 
as  low  as  possible,   or  alternatively  to  have  the  relay  zero -level  width 
as  narrow  as  possible,   the  gyro  time  constant,   J/C,   should  be  as 
low  as  possible.     Since  the  polar  moment  of  inertia  of  the  float  is 
usually  fixed  by  the  choice  of  angular  momentum,  the  time  constant 
depends  only  on  the  damping  coefficient,   C,   which  can  be  varied  by 
varying  the  operating  temperature  of  the  instrument  and  the  clearances 
between  float  and  case.     The  higher  the  damping,  the  lower  the  time 
constant,  but  the  smaller  the  angle  through  which  the  float  travels. 
There  is  a  certain  value,  A     .    ,  below  which  the  polarity  of  the  signal 
cannot  be  determined,   so  that  an  indefinite  increase  in  damping  is  not 
practicable.     A  short  computation,   suggested  by  J.    E.   Miller  of  the 
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M.  I.  T.  Instrumentation  Laboratory,   gives  an  expression  for  the 
optimum  damping  in  a  loop  with  a  two -level  relay,   considering  only 
the  dynamics  of  the  gyro,  which  usually  dominate. 
The  differential  equation  governing  the  gyro  is 

jAOA+  CAOA  *  HspWIA  +  Mtg  (3-5a) 

Integration  of  this  expression  over  an  integral  number  of  sampling 
periodsg  and  assumption  of  zero  initial  conditions  for  A~.   gives 

JAQA  +  CAQA=  Hgp(0IA  -  9.nd)  (3_5b) 

where  0,..  is  the  angle  through  which  the  gyro  case  has  turned  about 

its  input  axis;  and 
9.     ,  is  the  indicated  value  of  #Ta»  obtained  by  counting  output 

pulses. 
The  quantity  (0TA  -  9.    _,)  may  be  recognized  as  information  stored 

in  the  float  (see  Ch.   II,   Sec.   2.C.  4).     When  the  torque  switches  sign, 

the  value  of  A~A  is  A  and  the  value  of  A^A  is  M(l-a)/(l+a  )C. 

OA  min  OA  '  * 

.Is. 

Here  a  =  e      1    and  the  loop  is  assumed  to  be  moding  n-n.     Thus  at 
switching, 

j   =     JM(l-gn)  +  C 
S~VW>)      HspAmin  (3_5c) 

By  differentiating  Eq.    (3 -5c)  with  respect  to  C,   and  setting  the 
derivative  equal  to  zero,   an  extreme  value  of  I     may  be  found,   which 

can  be  shown  to  be  a  minimum.     The  dependence  of  a  and  n  on  C  may 
be  neglected. 

dla.      JM(l-a")+Amin   =0 
dC       HpdW>      Hsp 

or 
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:  =  I  "aI1)  J 


(l+an)A     .  (3-6) 

min  x        ' 

To  show  that  this  value  of  C  indeed  provides  a  minimum  in  I  , 

d2I  :  ; 

s  is  computed: 

dC2 

d%  =  2         JM(l-/),n 

dC2  H     C3(l+an) 

sp 

which  indicates  a  minimum. 

This  optimum  value  of  C  has  been  chosen  to  minimize  stored 

information.     The  corresponding  time  constant  is 

J         l(l+an)JA     . 

mm 

C      V     (l-»n)M  (3-7) 

For  loops  with  a  three -level  relay,  the  damping  should  be  chosen 
so  that 

Ac  =  Amin  (3-8) 

The  expressions  for  A     are  derived  in  Appendix  B,   Sec.   B.  1. 

It  is  possible  to  show  a  relationship  between  the  elastic  restraint 
present  in  a  gyro  and  the  threshold  of  the  instrument.     The  differential 
equation  of  the  gyro  is 

JAOA  +  CAOA  +  ^OA  ■  HspWIA +  Mtg  (3-9)  | 

where  K  is  the  elastic  restraint  acting  on  the  float  from  all  causes, 
including  electromagnetic  reaction  torque  and  the  effect  of  angular 
velocity  about  the  spin  reference  axis.     If  the  inertia-reaction  torque 
is  neglected  (it  is  usually  much  smaller  than  the  damping  torque)  and 
the  average  angle  of  the  float,  A~.,   is  considered,   limit  cycles  may  be 
ignored  (See  Ch.   II,   Sec.    2.D.I.  b)    Under  these  assumptions,   the  solution 
of  Eq.  (3-9)  is 

—  %WIA  "f* 

AOA  =     K  (1"e         }  (3-10) 
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In  time  T,  the  average  angle  of  the  float  reaches  a  value  AAa,   which 

represents  one  increment  in  the  loop's  output.     This  value  is  given  by 

A        =^T 
AA0       C     s 

Writing  Eq.    (3-10)  for  t  =  T  and  solving  for  T  yields 

sp    IA  (3-11) 

T  becomes  very  large  as(A  .  ^/II  W.)K  approaches  unity.     For  this 

term  greater  than  or  equal  to  one,   T  is  not  defined.     Physically,  this 
corresponds  to  the  case  where  elastic  restraint  completely  masks 
the  input  (See  Ch.   II,   Sec.   2.D.I.  c).     Using  only  the  first  and  second 
order  terms  in  the  series  expansion  of  the  logarithm,  for  small 
K*  gives 

CAA0  1        AA9 

T*  lv^(1+ Jlv^ K) 


(3-12) 


Since  MT    =  H     A9,  this  may  be  written 
s         sp  J 


WIA(2'WIA'^")  (3-13) 

Thus  the  percentage  error  in  time  of  occurrence  of  an  output  pulse 
is 

1       A0        K_ 

T"2'   WIA   '   C  (3-14) 

By  manipulating  Eq.   (3-14),  many  interesting  relations  between  the 

time  error  and  the  loop  parameters  can  be  seen.     The  error  in  time 

decreases  with  elastic  restraint  and  maximum  input  angular  velocity. 

It  is  inversely  proportiqnal  to  the  actual  input  angular  velocity  and 

the  square  of  the  damping  coefficient.     Since  the  optimum  damping 

coefficient  depends  inversely  on  A     .   ,   decreasing  A_.     decreases 
r  J  min  °     min 
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3.C    Delta -Modulated  Accelerometers 
3.C.1    Detailed  Description 

A  detailed  physical  description  of  a  pulse -restrained  accelerometer 
is  given  in  Reference  8.     Present  practice  is  to  center  the  float  in 
the  case  by  a  magnetic  suspension  system  described  in  Reference  25. 
Reference  4  gives  some  parameter  values  for  a  typical  accelerometer. 
Fig.   3-3  shows  a  schematic  drawing  of  the  instrument. 
The  accelerometer's  differential  equation  is 

JAOA+  CAQA  +  PapRA  sin  AQA  =  Pa^  cos  AQA  -  JW^ 

+  (1^  -  IpA)WIAWpA  +  Mtg  (3_15a) 

where  PA  and  PRA  refer  to  the  pendulum  axis  and  pendulum  reference 
axis  respectively.     (IA,   PRA,  OA)  is  a  right-handed  triad. 

The  last  term  on  the  right  hand  side  of  Eq.   (3 -15a)  is  an  anisoinertia 
term  which  is  negligible  in  practical  accelerometers  and  it  will  not 
be  considered  further  here.     Making  the  small  angle  approximation 
for  A^A  and  dividing  by  C  gives 

P  P  .  M 

7AOA  +  AOA  "  C"    aPRAAOA  =  c"    aIA    "  TWca  +  ~TT 

(3-15b) 

When  apRA  =  WQ     =0,   Eq.   (3-15b)  becomes  Eq.   (3-2b).     As  in  the 

case  of  the  gyro,  infinite  bearing  stiffness  has  been  assumed,   and 

torques  considered  only  about  the  output  axis. 

3.  C.  2    Rotational  Coupling  Effects 

As  can  be  seen  from  Eq.    (3 -15b),   only  rotations  of  the  instrument 
about  its  output  axis  affect  it,  because  of  the  assumed  stiffness  of 
the  bearings.     As  in  the  case  of  the  gyro,   angular  accelerations  of 
the  case  about  the  output  axis  look  like  inputs. 

If  the  case  motion  about  its  output  axis  is  vibratory  and  a  vibratory 
linear  motion  along  the  pendulum  reference  axis  occurs,   an  equation 
of  the  form  of  Eq.    (3-4)  governs  the  accelerometer  motion.     As  noted 
above,  if  the  frequencies  of  the  two  motions  are  harmonically  related, 
a  drift  in  the  loop's  output  may  appear. 
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3.C.3    Linear  Coupling  Effects 

Eq.   (3 -15b)  shows  that  an  acceleration  along  the  pendulum  reference 
axis  appears  as  an  elastic  restraint.     Imperfect  integration  results 
and  a  dead  zone  appears  in  the  instrument.     This  is  important  in  any 
application  of  this  accelerometer  since  it  may  experience  acceleration 
in  any  direction. 

If  a  vibration  occurs  along  the  pendulum  reference  axis  in 
conjunction  with  a.  vibratory  input,   a  drift  may  result  as  described 
above.     A  similar  situation  may  arise  if  vibrations  along  the  pendulum 
reference  axis  are  related  to  the  clock  frequency. 
3.  C.  4    Choice  of  Loop  Parameters 

Parameters  for  delta-modulated  accelerometers  may  be  chosen 
in  a  way  similar  to  that  used  when  choosing  gyro  parameters.     Those 
which  may  be  chosen  are  pendulosity,   torque  level,   clock  frequency, 
and  damping. 

The  pendulosity  of  the  accelerometer  is  chosen  on  the  basis  of 
the  level  of  torque  uncertainties  expected.     The  lower  the  torque 
uncertainty,   the  lower  the  pendulosity  may  be.     This  choice  involves 
the  same  sort  of  compromise  as  the  choice  of  angular  momentum  for 
a  gyro. 

The  required  torque  level  is  the  product  of  the  pendulosity  and 
the  highest  level  of  input  acceleration  expected.     This  value  is  practically 
always  lower  than  the  corresponding  value  for  a  gyro. 

The  clock  rate  is  determined  from  the  maximum  acceleration 
expected  and  the  velocity  increments  used.     Their  quotient  gives  the 
lower  bound  on  clock  frequency. 

The  choice  of  damping  coefficient  for  the  accelerometer  is  made 
just  as  in  the  case  of  a  gyro.  For  a  limit  cycling  accelerometer,  the 
optimum  damping  is 

1 


C  =i 


'(l-an)JM 


V  (1^n»Amin  (3-16) 


which  gives  a  time  constant 
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'(l+an)JA 

m: 


(l-<*n)M  (3-17) 

For  an  accelerometer  with  a  three -level  relay,   the  damping  coefficient 
should  be  chosen  so  that 

Ac  =  Amin  (3-18) 

The  remarks  about  elastic  restraint  made  with  respect  to  the 
gyro  apply  to  the  accelerometer  as  well.     The  equation  corresponding 
to  Eq.    (3-14)  is 

=  J_  _AV    K_ 
T      2     aIA     C  (3-19) 

3.D    Non -Ideal  Components 


Up  to  this  point,   very  little  has  been  said  about  the  non -ideal 
behavior  of  any  components.     The  following  sections  consider  some  of 
the  more  serious  hardware  problems. 
3.D.1    Gyros 

In  most  applications,  floated,   single -degree -of -freedom  integrating 
gyros  do  not  experience  significant  angular  velocities  about  any  axis. 
In  any  environment  in  which  they  do,  their  performance  will  deteriorate, 
simply  because  they  were  not  designed  to  take  such  punishment.     A  high 
angular  velocity  environment  affects  present-day  gyros  in  two  ways. 
First,   they  are  not  structurally  strong  enough  to  withstand  the  bending 
and  torsion  they  experience.     When  torque  is  applied  by  the  torque 
generator,   the  float  experiences  a  torsional  and  bending  stress  because 
of  the  high  spin  angular  momentum  it  has.     In  gimballess  applications, 
the  high  torques  involved  may  literally  tear  the  instrument  apart. 
A  second  effect  of  high  angular  velocities  is  high  bearing  loading.     The 
bearings  must  be  able  to  supply  reaction  forces  strong  enough  to  counter 
a  torque  equal  to  the  torque  generator  torque.     This  cannot  be  done  for 
long     periods  of  time,  but  in  a  floated  instrument  the  floatation  fluid 
provides  temporary    additional  bearing  support,   so  that  for  short  times 
the  float  will  stay  reasonably  well  centered.     While  at  the  present  time 
floatation  provides  necessary  additional  support  for  the  instrument, 


55 


magnetic  supports  are  becoming  capable  of  handling  the  entire  load 

alone. 

3.D.2    Torque  Generator  and  Current  Source 

In  the  analysis  above,   it  has  been  assumed  that  the  time  constant 

of  the  torque  generator  was  negligibly  small.     This  assumption  is 

violated  as  the  torque  from  a  given  instrument  is  increased.     Torque 

rises  as  the  square  of  the  input  current.     To  increase  this  current, 

either  the  voltage  must  be  increased  or  the  impedance  of  the  torquer 

decreased.     Increasing  the  voltage  is  not  an  attractive  solution,   since, 

in  general  transistors  drive  the  torquer,   and  their  voltage  limit  is 

quite  low.     Lowering  the  impedance  requires  lowering  the  resistance, 

which  in  turn  raises  the  time  constant  of  the  torquer,    L/R.     As  the 

time  constant  of  the  torque  motor  increases,  the  effect  of  the  lag 

makes  itself  felt  in  the  moding  plot.     The  value  of  the  torquer  time 

constant  above  which  the  lag  must  be  considered  is  about  T.     ="5"T  . 

tg  s 

More  torque  can,   of  course,  be  obtained  by  increasing  the  size  of  the 

torquer,   but  it  is  desirable  to  have  small  instruments. 

Microsyn  torquers  (Ref.  28)  have  been  successfully  used  in  pulse 
torquing  applications.     However,   they  are  not  high  torque  instruments, 
and  their  dependence  on  current  squared  makes  them  non -linear  when 
used  with  a  three -level  relay. 

Flapper-type  torquers  (Ref.    23)  provide  about  ten  times  as  much 
torque  as  a  corresponding  size  microsyn.     Their  main  disadvantage 
is  the  presence  of  a  fairly  large  negative  elastic  restraint,    which 
makes  the  torquer  positionally  unstable.     This  might  be  eliminated  by 
providing  a  positive  mechanical  or  electrical  spring,   or  using  ac 
voltage  for  torquing  and  tuning  out  the  restraint  term.     The  ac  torquer 
provides  a  double  frequency  term  as  well  as  a  constant  average,   so 
that  the  frequency  used  must  be  well  above  the  clock  frequency,   and 
above  the  mechanical  break  point  of  the  system. 

Moving -coil,   or  d'Arsonval,   instruments  also  provide  large  torques 
for  small  size.     This  type  of  torquer  usually  uses  a  permanent  -magnet 
field.     The  difficulty  with  permanent  magnets  is  that  measurement 
techniques  are  not  yet  accurate  enough  to  determine  the  stability  of  the 
field. 

A  problem  that  all  three  kinds  of  torquers  share  is  the  stability 
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of  their  sensitivities.     Good  operation  of  the  instruments  requires  that 
the  same  angular  impulse  be  applied  to  the  float  for  each  electrical 
impulse  applied  to  the  torquer.     If  this  is  not  the  case,  the  instrument 
scale  factor  will  vary  with  time,   and  errors  will  appear  from  this 
source. 

A  related  problem  is  the  variation  in  torquer  time  constants  with 
direction  of  torque.     If  the  time  constants  are  different,   more  impulse 
will  be  applied  in  the  direction  of  the  smaller  time  constant.     Eventually 
a  change  in  the  output  will  be  indicated  although  the  input  has  been  zero. 

To  achieve  the  required  accuracy  in  torquing .  a  very    , 
stable  current  source  must  be  used.     If  a  three-level  relay  is  used, 
either  the  current  must  be  turned  off,   or  it  must  be  fed  into  a  dummy 
load.     Turning  off  the  current,   while  achieving  great  savings  in  power, 
causes  stability  problems  in  the  current  source  above  those  already 
present.     Ac  current  sources  can  be  made  more  stable,  but  ac  torquers 
are  less  sensitive  than  dc  torquers. 
3.  D.  3    Signal  Generator 

With  ac  signal  generators,  the  output  is  a  supressed- carrier  wave. 
Thus,    careful  attention  must  be  paid  to  the  effects  of  circuitry  which 
operates  on  the  signal  generator  signal  so  that  the  gain  for  each 
sideband  frequency  is  the  same,   and  the  phase  shifts  at  each  sideband 
are  the  negative  of  each  other.     If  this  is  not  the  case,   a  quadrature 
component, which  will  cause  a  larger  uncertainty  in  the  position  of  the 
float  with  respect  to  the  case,  will  be  introduced  into  the  signal.     The 
end  result  will  be  an  increase  in  the  minimum  angle  of  the  float  which 
can  be  detected.     For  this  reason,  the  low-pass  to  band-pass  trans- 
formations commonly  used  in  ac  compensation  are  not  sufficiently 
accurate  to  be  used  in  this  application. 

As  noted  above  (Sec.   3.B.4)  when  ac  signal  generators  are  used, 
the  carrier  and  clock  pulse  phases  are  adjusted  so  that  the  clock  pulses 
occur  at  the  peaks  of  the  carrier  sinusoid.     This  adjustment  is  usually 
made  by  shifting  the  ac  signal  at  either  the  primary  or  the  secondary 
winding  of  the  signal  generator.     For  the  reasons  noted  above  this 
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phase  shifting  should  be  done  on  the  primary  to  avoid  adding  a 
quadrature  component  to  the  signal  generator  and  increasing  the 
minimum  discernable  angle. 

Flapper-type  (Ref.   24)  and  d'Arsonval  signal  generators  provide 
attractive  alternatives  to  microsyn  signal  generators,  because  of 
their  small  size. 
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CHAPTER  IV 


KINEMATICS  OF  ROTATING  COORDINATE  FRAMES 


4.  A    Introduction 

In  a  conventional  inertial  guidance  system,   the  inertial  space  reference' 
is  preserved  physically  by  the  gimballed  stable  platform.     The  gimbals 
isolate  the  platform  from  the  motion  of  the  vehicle  through  servo  drives, 
operating  to  null  the  output  of  gyros  mounted  on  the  stable  platform.     In 
a  gimballess  system,  the  inertial  space  reference  may  be  preserved 
analytically  in  a  transformation  computer.     This  computer  operates  on 
a  vector  measured  in  body-fixed  coordinates  (the  A.V  vector)  to  give  the 
components  of  the  vector  in  an  inertially  fixed  frame  of  reference.     The 
computer  uses  angular  information  obtained  from  gyros  mounted  on  the 
vehicle  to  generate  the  rotational  transformation. 

In  this  chapter,  the  properties  of  the  transformation  and  the  relations 
between  the  angular  velocity  of  the  body-fixed  coordinate  system  with 
respect  to  inertial  space  and  the  rates  of  change  of  the  transformation 
quantities  will  be  adduced. 

There  are  essentially  three  methods  of  representing  rotational  trans- 
formations.    Probably  the  most  familiar  is  the  matrix  of  direction  cosines. 
This  matrix,  or  alternatively,  this  second-order  tensor,   may  have  its 
elements  written  in  terms  of  the  nine  direction  cosines  explicitly,  in 
terms  of  a  set  of  Euler  angles,  or  in  terms  of  the  axis  and  magnitude  of 
the  rotation.     Quaternions,   invented  by  Sir  William  Hamilton  in  1853  (Ref.  31), 
use  four  parameters  to  describe  a  rotation.     The  vector  to  be  transformed 
is  cast  in  quaternion  form,   pre -multiplied  by  a  quaternion  representing 
a  rotation,   and  post -multiplied  by  the  conjugate  quaternion.     The  result 


is  a  quaternion  representing  the  rotated  vector.     A  variation  on  this,  the 
Cayley -Klein  parameters,  uses  the.  same  four  parameters  in  a  slightly- 
different  form.     The  third  transformation  uses  a  vector  to  represent  the 
rotation,   and  performs  the  transformation  by  means  of  purely  vector  and 
algebraic  operations. 

These  three  types  of  transformation  will  be  compared  to  determine 
which  is  most  useful  for  machine  computation. 
4.  B    Direction  Cosines:    Matrices  and  Tensors 
4.  B.l    Direction  Cosines 


The  most  familiar  of  the  rotation  transformations  is  the  matrix 
transformation  (Ref.    30;  p.   95ff) 

Vj    =    CV^  (4-1) 

where  V-,  is  a  column  vector  with  components  measured  in  an 

inertial  coordinate  frame. 

V,    is  the  same  vector  with  components  measured  in  a 
body -fixed  coordinate  frame,   and 

C  is  the  square  matrix  of  direction  cosines  of  the  inertial 
axes  relative  to  the  body  axes. 
The  first  component  of  Eq.   4-1  in  vector  notation  is 

Vn  =Tn  *  \  _  _  J 

=  (1n  *   H)i)vbi  +  (1n  '   1b2)vb2  +  (1n  '   Wvb3  i 

(4-2) 
where  L.,  is  the  unit  vector  along  the  inertial  X  axis,   etc.     Since  the 
vectors  involved  are  unit  vectors, 

lTL    '     H)l    =    COS  (1I1    »     W    =    Cll  (4-3) 

and  in  general,  the  elements  of  the  matrix  C  are 

cij =  cos  (ii  -  y  =  hi  -  Vj  (4-4) 

The  rate  of  change  of  the  elements  of  C  may  be  found  by  differentiating 
Eq.   (4-4)  directly: 

dt       ij  dtl    Ii  bj)        dtl   III  j  bj  Ii         dt  I  bj  /j, 
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=  1T.     •      [  (?.   L  .).     +    oo    X    L  .  ] 

Ii  L    dt     bj  b  bj  J 

=  (TIi    '  "WV«    ~    ^Ii    "     V^j+l  (4-5) 

or  using  a  dot  (• )  to  denote  time  differentiation, 

Cij    =    Cij  +  lV2    "    Cij  +  2Wj+l  (4-6) 

Note  that  oo.  is  the  j       component  of  angular  velocity  of  the  body -fixed 

J 

coordinate  frame  measured  in  body  coordinates. 

Eq.   (4-6)  may  also  be  derived  by  differentiating  Eq.    (4-1)  directly: 

=  c(vb  +  c'lb\)  (4.7) 

-1  *  — 
This  is  the  Law  of  Coriolis  in  matrix  form.     If  the  term  C     CV,    be 

b 

identified  with  oo  X  V,,   and  the  angular  velocity  matrix  ft  be  defined  by 

0 
ft    = 


(4-8) 

then 

5'XV,i-flV,    =  c~lcvu.  ,.   n, 

b  b  b  (4-9) 

Thus,   the  differential  equation  for  C  is 

C  =  -Cft.  (4-10) 

which  expands  to  give  Eq.   (4-6). 

There  are  nine  direction  cosines.     Since  a  rotation  involves  only 
three  degrees  of  freedom,   there  must  be  six  constraints  on  these  nine 
quantities.     These  constraints  arise  from  the  fact  that  during  a  rotation, 
the  magnitude  of  a  vector  remains  unchanged.     This  condition  requires 

that 

T  -1 

C1     =    C  (4-11) 

T 
where  C     is  the  transpose  of  C 

and       C      is  the  inverse  of  C. 

Another,   more  useful,   form  of  the  orthogonality  condition  is 
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CTC    =    CCT    =    I  (4-12) 

where  I  is  the  identity  matrix.     Thus  it  is  possible  to  write  the  nine 
direction  cosines  as  functions  of  any  three  of  them,   using  six  independent 
orthogonality  relations  from  Eq.   (4-12).     Another  solution  involves  the 
use  of  Euler  angles. 
4.  B.  2    Euler  Angles 

There  are  many  ways  of  defining  Euler  angles;  the  definition 
shown  in  Figure  4-1  follows  the  most  common  usage  (Ref.    30  p.   104ff.  , 
and  Ref.   39,  p.   23,   for  example).     The  first  rotation  is  about  the  Z  axis 
of  magnitude  A„.     The  second  rotation  is  about  the  displaced  x  axis,   which 
becomes  the  line  of  nodes,  by  an  amount  AT  .     The  final  rotation  is  about 
the  z  axis  by  an  amount  A   .     In  example  being  used,   the  capital  letter 
axes  correspond  to  the  inertially  fixed  frame,   while  the  lower  case  axes 
are  the  body-fixed  set. 


LN 


Fig.   4-1    Euler  Angles 
In  terms  of  these  Euler  angles,   the  direction  cosine  matrix  is 
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8inAL 

sin  Ax 

«lnAL 

cos  A, 

co.  AL 

(4-13) 


The  rates  of  change  of  the  Euler  angles  in  terms  of  the 

angular  velocity  of  the  body-fixed  set  measured  in  the  body-fixed 

coordinates  may  be  determined  by  resolution,   and  are  (Ref.   39; 

p.   24) 

sin  A  cos  A 

.  z    , 

A,~    =  u>^ — : t —  +  <*> 


Z         1  sin  AT  2   sin  AT 

AT    =  co    cos  A     +  w0  sin  A 
LI  z         2  z 

A     =  «,  cot  AT    sin  A     -  w     cot  AT    cos  A     +  w0 
zl  L,  z         A  L  z         3 

(4-14) 
This  computation  suffers  from  two  disadvantages;  trigono- 
metric functions  must  be  determined,   and  a  singularity  exists 
at  sin  AT    =0.     The  singularity  corresponds  to  gimbal  lock  in 
a  three-gimbal  system  (see  Ref.    39,   Figure  8-37,   p.    25)  and 
can  be  eliminated  by  including  a  fourth  parameter  in  the 
description  of  the  rotation,  just  as  gimbal  lock  may  be  eliminated 
by  the  addition  of  a  fourth  gimbal. 
4.B.3    Tensors 

A  third  expression  for  the  direction  cosines  is  most 
conveniently  cast  in  tensor  form.     It  should  be  remarked  here 
that  since  the  transformation  involves  only  Cartesian  coordinates, 
no  distinction  between  covariant,    contravariant,   and  physical 
components  of  vectors  and  tensors  need  by  made  (Ref.   35,   p.   170), 
The  direction  cosine  tensor  is 

Cij  =  C°S  A  6ij  +  (1"C°S  A)  rirj  "  Sin  A  Gijkrk  (4-15) 
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where  the  rotation  which  carries  the  body -fixed  axes  into  the  inertial 
axes  is  about  an  axis  with  direction  cosines  r.  (the  same  with  respect 
to  either  coordinate  system)  with  magnitude  A.     Here  6..  is  the 

J 

Kronecker  delta 

(4-16) 


e..,    is  the  cyclic  tensor 


[+1  if  ijk  positive  cyclic  (e.  g.   123) 


e..,    =  <  -1  if  ijk  negative  cyclic 

0  if  i  =  j,     j  =  k,     or  i  =  k  (4-17) 

and  the  convention  of  summing  on  repeated  indices  is  observed  (Ref.    34, 
p.   96  for  example).     Using  the  tensor  formulation,   Eq.   (4-1)  is 

VT.    =    CV,  .  IA   1Dv 

Ii  ij    bj  (4-18) 


Cij  Cikekjmwm  (4-19) 


and  Eq.   (4-10)  becomes 

4.  C    Four  Parameter  Representation 

The  direction  cosine  computations  are  six -fold  redundant;  the  Euler 
angle  computations  carry  no  redundancy  but  have  a  singularity.     The  use 
of  four  parameters  to  describe  a  rotation  eliminates  the  singularity  and 
diminishes  the  redundancy.     There  are  two  practically  identical  ways  of 
describing  a  rotation  with  four  parameters.     These  are  Hamilton's 
quaternions  (Refs.     31-33,    36,   38,   40)  and  Cayley-Klein  parameters 
(Refs.    30,   40). 
4.C.1    Quaternions 

Hamilton  developed  quaternions  as  an  extension  of  complex  numbers. 
In  representing  a  vector  in  two  dimensions  as  a  complex  number,   the 
unit  real  number  1  represents  one  direction,   and  the  unit  complex  number 
i  =  n/-1,  the  other.     In  analogy  with  this,  Hamilton  let  the  unit  space 
vectors  be  i,  j,  k,   with  multiplication  rules  similiar  to  those  for  *s/  — 1: 

•;     -i2=  j2  =k2  =  -1      ij  =  -ji  =  k,   etc.  (4-20) 

Thus,    a  space  vector  may  be  represented  as  V  =  iV,  +  jV~  +  kV„,   and 
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a  general  quaternion,   with  the  unit  real  1     the  fourth  direction,  as 

q  =  a+  ib  +  jc  +  kd  (4-21) 

Conjugation  is  defined  by 

q"C  =  a  -  ib  -  jc  -  kd  (4-22) 

and  length  is  given  by 

9.  *  9  9  9  9 

r  =  q  q  =  *    +  h'  +  c*  +  d*  (4-23) 

A  simpler  notation  follows  from  representing  a  quaternion  as  a  scalar 
and  a  vector: 

q  =  (X,   p)  =    X+  ipl  +  jp2  +  kp3  (4-24) 

Thus  i 

q*  =  (X,    -p)  (4-25) 

and  if  q^  =  (X^  p^)  and  q2  =  (X2,  p"2),   then 

=  (xxx2  -  prp2,  xlP'2  +  x2px  +Pxxp2)  (4_26) 

Note  that  in  general  multiplication  is  not  commutative.  .     Two  special 
quaternions  may  be  defined,   the  unit  quaternion  q  ,   and  the  zero  quaternion 


*u 


V 


qu  =  (l,0)  qo  =  (0,6)  (4.27) 

Obviously 

qqu  ■  quq  ?  q       qoq  =  qqo  =  qo  (4_28) 

Further,   if  q   f  q  ,  there  exists  a  unique  inverse  or  reciprocal  q      such 
that  qq       =  q     q  =  qu= 

-1         1  * 

q     =  — o  °i 

Iqr  (4-29) 

Two  relations  follow  from  the  definitions  of  conjugate  and  inverse: 

(qxq2)  *  =  q2*<li*  ^2]~    =  V    ^2  (4-30) 


If  the  quaternions  qj  and  %  ^q  defined  by 
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*i =  (ar  vi)}       %  ~~  (V  V,  (4_31) 


and  related  through  the  quaternion  q  having  unit  magnitude. 

q  = 
by  the  transformation 


q  =  (X,p )        X2  +   p2  =  1  (4-32) 


qv  =  qqbq*  (4-33) 

then 

(arVI)  =  [ab'    Vb+2X^xTb  +  2^X(^X^)]  (4-34) 

The  scalars  of  qT  and  q,    are  equal,   and  the  vectors  are  related  by 

Vj=Tb  +  2XpXTb  +  2pX(pX-b)..  (4.35) 

Writing  out  Eq.    (4-18)  in  vector  notation,   and  using  r  as  the  unit  vector 
along  the  axis  of  rotation  gives 

Vj  =  V~h  cos  A  +  (r-  V^)  f  (1-cos  A)  +  rX?bsinA 

■'JT   +  r  X  Vu  sin  A  +  [  (r-  V,  )r  -  (r«  r)VJ  (1-cos  A) 

O  D  D  D 

=  Vb  +  r  X  V^  sin  A  +  r  X  (r  X  Vb)(l-cos  A) 

=  V"b  +  2  sin  ^-  cos  "^    r  X  V"b  +  2  sin2  ^rX(rXVj 

(4-36) 

A 
If  X  in  Eq.    (4-35)  is  identified  with  cos -s- in  Eq.   (4-36),   and  p  with 

-  A 

r  sin-rj-,   then  the  transformation  Eq.   (4-32)  corresponds  to  the  rotation 

of  V,    about  an  axis  r  through  an  angle  A.     Writing  q  out  gives 

q  =  (cos-^,   r  sin-^)  =  cos-^  +  sin-^  (ir^  +  jr2  +  krg) 

Extending  Euler's  formula  for  exponentials  of  complex  quantities 
suggests  the  notation 

4(ir1  +  jr2  +  kr3)  * 

e  =  cos-s-+  sin-^  (ir,  +  jr9  +  krQ) 

2  l         I  2.  6  (4_38) 

2       2  2 

with  r,  +r2     +  r„     =1,   so  that  the  transformation  Eq.   (4-33)  may  be 

written 
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i(irl  +  Jr2  +  kr3>        'T  (irl  +  Jr2  +  kr3) 
qx  =  e  qbe 

The  quaternions  representing  successive  rotations  may  be  combined  as 
follows: 

^2  =  VqA,qlqi,  ]   q'V  (qA9V  ql  (qA9V  ' 

1        l  i  *  2      i  2      l  (4-40) 

The  non-commutativity  of  finite  rotations  is  evident  from  the  non- 

commutativity  of  quaternion  multiplication. 

The  rates  of  change  of  the  rotation  quaternion  parameters  must  now 

be  related  to  the  angular  velocity  of  the  body -fixed  coordinate  system 

with  respect  to  inertial  space.     Let  q  be  the  quaternion  representing  the 

rotation  that  carries  the  body -fixed  space  into  the  inertial  frame  at  time 

t.     The  quaternion  representing  the  opposite  rotation  is  then  q'.     Let 

the  quaternion  representing  the  rotation  carrying  the  body-fixed  frame  into 

inertial  space  at  time|t  +  A  t|be  q'.     The  quaternion  representing  the 

rotation  of  the  body-fixed  axes  during  the  time   At  is  therefore  q'q*: 

AA  AA   - 

q'q*  =  (cos  —£-,   sin  -—^-  r^)  =  (X  +  AX,     p+  Ap)(X,_p) 

=  [(X,p)  +  (AX,Ap)](X,  -p) 
=  (1,0)  +  (AX,  Ap)(X,-p) 


(4-41) 


or 


(cos  ^-1,   sin^A_r^t)  =  (AX,Ap)(X, -p)  (4-42) 

Dividing  this  equation  by  At  and  allowing  At  to  approach  zero  gives 


lim       1 
A  t  -*0      2 


AA      ,  .     AA 

cos      2     -1   AA  sin 


A  A/2  At'/  At       AA/2 


v-  Sm     2       AaI      ,AXApw.         -, 

>rAt  ~aa72  '"StJ,-  (3WX'  "p> 


=  2   <°>  t^   r  At  "ST*  =  (*'P  HX'  "P)  (4_43) 

Since    *^    r  M  ^     A  w,    Eq.    (4-43)  becomes 
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2 


<0,w)  =  (X,    p)(X,-p) 


(4-44) 


or,   post -multiplying  by  (X,p) 


(Xj)  =2  (°^)(*>P) 


whence 


and 


(4-45) 


1    -    - 
X  =  -  2   wp 


p    =   *  (Xw  +  wXp)  (4-46) 

4.  C.2    Cayley-Klein  Parameters 

The  foregoing  discussion  provides  a  method  for  describing  a 
rotation  with  four  parameters.     One  is  redundant;  the  redundancy  is 
removed  by  the  constraint  that  the  quaternion  q  =  (X,p)  have  unit 
magnitude  [  Eq.   (4-32)]  .     This  formulation  suffers  from  the  fact  that 
quaternions  are  not  at  all  widely  known.     To  obviate  this  difficulty, 
Cayley-Klein  parameters  may  be  introduced;  they  retain  the  simplicity 
of  quaternions  while  eliminating  Hamilton's  i,  j,  k.     Defining  the  Cayley- 
Klein  parameters  by 

a    =  X  -  ipQ 


0  =  -P2  -  iP3 
Y  =  +P2  -  iP3 
6     =  X  +  ip3 

and  from  them,  forming  the  complex  matrix  Q, 

fa  j3\  ( q      j3 

kY*J 


(4-47) 


Q  = 


(4-48) 


(4-49) 


gives  a  matrix  which  is  unitary  and  has  determinent  +1: 
*  A  -1 

Q       =  QT  =  Q 

I  Ql    =  aa*  +  j8£*  =  X2  +  P3    +  P%  +  P*   =  X2  +  p2  =  +1 

Here,   Q     is  the  complex  conjugate  of  Q,   Q     is  the  transpose  of  Q,and 
Q     is  the  adjoint  of  Q.     When  using  Cayley-Klein  parameters,   vectors 
are  represented  as  Hermition,   or  self -adjoint,   matrices: 
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V~~V  =        V3  iv2|  (4-50) 


iv1  +  iv2  .3 

Note  that  the  trace  of  the  matrix  V  is  zero,   and  the  determinant  of  V  is 

the  negative  of  the  square  of  the  magnitude  of  the  vector  V: 

r2   ,   Tr2   ,   Tr2x  _     —2 

(4-51) 


VI    =  -(V|  +  V^  +  vp  =  -V' 


Using  the  matrix  Q,   a  similarity  transformation  may  be  performed 
on  the  matrix  V  ,  yielding  another  matrix  V-,: 

Vj  =  QVbQA  (4-52) 

Under  similarity  transformations,  the  Hermition  property  of  a  matrix  is 
preserved,  as  are  the  trace  and  determinant.  Therefore,  V,  must  be  of 
the  same  form  as  V,    [  Eq.   (4-50)]  ,   and  further 

~Vb  =  '  Vb'    =  IVl'    =  ~V\  (4"53) 

Eq.   (4-53)  is  the  orthogonality  condition.     Thus  the  transformation 
Eq.    (4-52)  is  equivalent  to  the  transformation  Eq.   (4-33). 

The  differential  equations  for  the  Cayley-Klein  parameters  follow 
directly  from  the  defining  equation  [  Eq.   (4-47)]    and  the  rates  of  change 
of  the  quaternion  parameters  [Eq.   (4-46)J  . 
4.D    Vector  Representation 

Another  useful  representation  of  a  rotation  is  the  vector  o 

cj  =  r  tan    .  (4-54) 

where  r  and  A  are  defined  above.     In  terms  of  the  scalar  and  vector  of 
the  rotation  quaternion 

_       _    sin  A        _    sinA  p  p 

a  =  r  r~  =  r  — 


cos^     2cos  -7    1  +  cos  ~2  (4-55) 

The  inverse  relations  are 

=  cosT=  2  cos  —  -1  = 2T-I  % 2~A 

sec  -j-       sec  -7- 
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2  -(l  +  tan2^|)      j  _a2 
1  +  tan2  A  1  +  a2 

and  (4-56) 

p  =  (1  +  X)  a  =     2     .       B 

Substituting  these  expressions  for  X  and  p  into  Eq.   (4-35)  gives 

2  -  - 

V"T  =  V:+    2    *  "ff        (     2g  9)XTh  +  2    (— 22-s)  X[  (-22L-)  X  V.] 

1         b  1  +  a2        1  +  a2  b  1  +  a2  1  +  a2  b 

or 

2 

VT  =  V,    +     4(1  "  °>  t  aXf    + §-_  ?X  (aX  V,  )  (4-57) 

1         D  (1  +  aV  D        (1  +  aV  D 

The  vector  representation  of  a  rotation  requires  only  three  numbers; 
there  is  no  redundancy.     There  is  a  singularity,  however,   at|AJ  =  2n , 
corresponding  to  a  full  rotation.     This  singularity  may  be  eliminated  by 
defining  another  vector  a': 
a 


(4-58) 


— ~rA    =  (  -    r)  cot  -£•=  (-r)  tan  (y  -  -g) 

tanT 


a'  =  (-r)  tan    2lt  ~K  (4-59) 

Since  a  rotation  of  A  in  one  direction  is  equivalent  to  a  rotation  of 
27r    -A  in  the  opposite  direction,   a  and  a'  describe  the  same  rotation,   and 
Eq.   (4-57)  holds  with  a'  substituted  for  a. 

The  differential  equation  relating  the  rate  of  change  of  a  to  the 
angular  velocity  of  the  body -fixed  coordinate  frame  may  be  obtained  by 
differentiation  of  Eq.   (4-55): 

a  =    -i iJ— 

1  +  x      <l  +  a;)? 
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(4-60) 


_  1  f  Xoj  +  cjjXp  (cj  »p)p      1 

1      l-o2 
=  s-  [  — o —   co  +  co  X  a  +  (co«  a)a] 

=  ~  [  ^y~  o~>  +  w  X  a  +  (coXa)  X  a  +  a2u] 

J  *     ,      ,2 

d  =  j[  -^y-  u  +  uXa+(uXo)Xa] 

Clearly,   the  same  differential  equation  holds  for  a': 

2 
a1  =  ~[-^~-£+  wX  a1  +  (wXct')  X  a']  (4-61) 

4.  E    Comparison 

The  preceding  three  sections  of  this  chapter  have  presented  several 
methods  for  transforming  a  vector  from  one  Cartesian  coordinate  frame 
to  another,   and  derived  the  differential  equations  for  the  transformation 
parameters  in  terms  of  the  angular  velocity  of  one  frame  with  respect 
to  another.     These  transformations  and  their  differential  equations  are 
collected  below  for  comparison. 
4.  B    Direction  Cosines 


Transformation 


Vj    =    CVb  (4-1) 


Direction  Cosine  Differential  Equation 

C    =    -C«  (4-10) 

Six -fold  Redundant,   with  Six  Constraints 

CTC     =    I  (4-12) 

Euler  Angle  Differential  Equation 

sin  A  cos  A 

z  z 

AZ  =  wl  sin  AL  +  W2     sinAL 

AT     =  u>,  cos  A    +  co     sin  A  (4-14) 

LI  z         2  z 

A      =    -co.  cot  AT    sin  A     -w     cot  AT    cos  A    +  oo 
zl  L  z2  L>  z         o 
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No  redundancy;  singularity  at  sin  AT    =  0 

4,C    Four  Parameters 

Quaternion  Transformation 

V~  =  Vb  +  2Xp  X  Vb  +  2p  X  (p  X  V,) 


A 
2 


X  =  cos   it       P  =  r  sin  -77 


2 


(4-35) 
(4-37) 


Cayley -Klein  Transformation 


X  -   ip3         -p2  -  lp, 


p2  -    lPj  X  +ip3 


b3 


X  +  lp3      p2  ♦   lp1( 

7a  *  V\  x  -  ip. 


(4-52) 


Differential  Equation 


1 


X  =  -  j  w.  p 


1 


p   ~2  (^w  +cjXp  ) 
Once  redundant  with  one  constraint 


X2+P2=l 


4.D  Vectors 


Transformation 


(4-46) 


(4-32X4-49) 


VI  =  Vb  +  4      1_22°  XVb  +  — ~  « 
1  b  (1  +  o  V  b       (1  +  a") 


2-2°X(aXVb) 


(4-57) 


a  =  r  tan 


(4-54) 


or 


vi  =  vb + 4  ,,l:a\2  ° x  vb + -  ^  ,2,2  s' x  <* x  vb) 


(1  +  a'")' 


(1+  o'T 


a'=     -A 


(4-57a) 
(4-58) 
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Differential  Equation 

1       1  +     2 
o    =  y[  g      w  +  w  X  -a  +   (<4Xa)  X  a]  (4-60) 

°f  i  2 

a1  =  j[-^i-^-(L  +  wXa'  +  (wXa')X  a1  ]  (4-dl) 

These  transformations  will  now  be  examined  with  a  view  to  determining 
which  is  most  suitable  for  machine  computation.     The  basic  criterion  to 
be  used  in  this  determination  is  the  size  of  computer  required.     Further, 
since  the  information  being  dealt  with  is  incremental,   incremental 
computation  will  be  the  basis  of  decision  unless  general  propose  computation 
seems  particularly  appropriate. 

The  most  straightforward  computation  employs  the  direct  direction 
cosine  formulation.     The  transformation  requires  nine  incremental 
integrators,    since  each  component  of  the  vector  to  be  transformed  has 
magnitude  AV.     The  generation  of  the  direction  cosines  requires  eighteen 
incremental  integrators.     If  Euler  angles  are  generated,   then  the  sine 
and  cosine  of  each  angle  is  also  required;  further,   the  sines  and  cosines 
must  be  combined  to  give  the  direction  cosines  as  well  as  the  factors  in 
the  differential  equations.     The  use  of  incremental  devices  requires 
forty-six  integrators.     This  does  not  handle  any  problems  raised  by 
the  singularity  at  sin  AT    =  0 

Twelve  incremental  integrators  are  necessary  to  generate  the  four 
parameters  used  in  the  quaternion  and  Cayley-Klein  parameter  representation. 
If  the  direction  cosines    are  formed  from  these  parameters,   an     additional 
twelve  integrators  are  required.     Nine  more  complete  the  transformation. 
Thus,   thirty -three  incremental  integrators  are  required  in  all,   as  opposed 
to  twenty -seven  for  the  direction  cosine  computation.     If  the  formal 
Cayley-Klein  similarity  transformation  is  used,   the  first  multiplication 
can  be  performed  by  twelve  incremental  integrators.     An  arithemetic 
computation  seems  best  for  the  second  matrix  multiplication;  it  requires 
twelve  multiplications  and  nine  additions. 

Twenty-one  incremental  integrators  are  required  to  generate  the 
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three  components  of  the  vector  used  to  represent  a  rotation;  the  vector 
transformation  requires  twenty-five  more.     In  addition,   an  arithmetic 
section  is  required  to  perform  the  branching  operation  used  to  remove 
the  singularity. 

The  foregoing  discussion  indicates  that  the  direction  cosine 
computation   using  a  total  of  twenty -seven  integrators  is  most  compact. 
For  this  reason,  this  formulation  will  be  the  only  one  considered  in 
designing  the  transformation  computer. 
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CHAPTER  V 


TRANSFORMATION  COMPUTER 


5.  A     Introduction 

In  Chapter  IV,   the  various  ways  of  representing  coordinate  ro- 
tations were  derived  and  compared  to  determine  their  suitability  for 
machine  computation.     As  a  result  of  this  comparison,   the  matrix  trans- 
formation was  chosen  to  be  instrumented.     The  rotation  matrix  is  to  be 
generated  by  direct  computation  of  the  direction  cosines.     The  matrix 
differential  equation  to  be  instrumented  is 

C  =  -CO  (5-1) 

where  C  is  the  matrix  of  direction  cosines,   a  dot  (•  )  over  a  quantity 
denotes  time  differentiation,  and  Q  is  the  angular  velocity  matrix, 

o 
"wi 

The  angular  velocity  is  resolved  into  body  axes. 

Eq.   (5-1)  represents  nine  separate  equations;    there  are  three    ■ 
sets  of  the  form 

Cil.=  Ci2w3  "Ci3u2 

^i2=   Ci3wl  -Cilu3  (5_3) 


Ci3  =   Cilt02  "  Ci2wl 


Multiplication  of  Eq. _<5-3)  by  dt  and  integration  of  the  result  yields 

Cil  ■  >  Ci2d83  "  '  Ci3d92 

Ci2  "  J  Ci3d8l   "  >  Cild83  <5"4> 

Ci3  =  /  CudOa  -  /  Ci2d91 

5.B    Algorithms  Available 

Equations  (5-4)  show  that  each  direction  cosine  is  the  difference 
of  two  integrals,  each  the  integral  of  a  direction  cosine  with  respect 
to  rotation  about  a  body  axis.    Since  the  gyro  outputs  are  equal  incre- 
ments of  rotation  about  body  axes,   the  integrals  may  be  approximated 
by  quadrature  rules  using  equal  AG  increments. 

5.B.1    Quadrature  Rules 

Almost  any  standard  text  on  numerical  analysis  (Refs.   41, 
46-48,   52,   54,  for  example)  gives  the  more  common  quadrature 
rules.     Since  only  the  integrand  will  be  available,  no  rules  employing 
derivatives  of  the  integrand,  such  as  Taylor'  s  series,  will  be  con- 
sidered.    The  most  common  quadrature  rules  using  present  and  past 
values  of  the  integrand  are  the  rectangular  rule, 


9(n) 

C(e)d0  =  A6C[e(n-l)] 
e(n-l)  (5-5) 


1 


the  trapezoidal  rule, 

e(n) 


I 


C(8)d9=^{:c|!e(n-1)]  +  C[9(n)]  }  (5-6) 


0(n-l) 
and  Simpson'  s  one-third  rule, 


76 


G(n)    : 
fc(9)d9=^  {C[9(n-2)j+4.C[9(n-l)]  +  C[9(n)]  }  (5-7) 

9(n-2) 

The  first  is  an  open  rule;   that  is,    the  value  of  the  integrand  at 
the  present  time  does  not  appear.     The  remaining  two  are  closed 
rules.     Each  integration  rule  fits  a  (p-1)  degree  polynomial  through 
the  p  points  of  the  integrand  appearing  in  the  formula;   for  this  reason, 
these  rules  are  called  (p-1)  order  polynomial  rules.     The  rectangular 
rule  is  called  a  zero-order  rule;    the  trapezoidal  rule  and  Simpson'  s 
one-third  rule  are  first-  and  second-order  rules  respectively. 

At  each  step,    the  approximation  to  the  integral  is  in  error  by 
some  amount.     This  error,   arising  from  not  considering  higher 
derivatives  or  differences,   is  called  truncation  error.     For  the  three 
rules  above,   the  truncation  error  made  at  each  step  is  (Ref.   47,  p.  7  3)  ; 

Rectangular      E     =  0(A92) 

Trapezoidal     E     =  0(A93) 

Parabolic  E     =  0(A95) 

n 

This  notation  means  that  for  the  rectangular  rule,   the  truncation  error 

2 
at  each  step  is  on  the  order  of  A0    .     The  term  "on  the  order  of" 

2  2 

means  that  the  limit  of  E    /A9     remains  finite  as  A9     approaches  zero. 

5.  B.  2    Influence  of  Size  of  Step 

The  rather  fast  pulse  rate  required  in  the  gyro  loops  suggests 
that  using  higher  order  integration  rules  with  larger  increments 
would  reduce  the  speed  requirement  on  the  computer  as  well  as  de- 
creasing the  truncation  error.     For  example,   if  rectangular  integra- 
tion were  used  with  a  step  length  A9,   the  truncation  error  at  each  step 

2 
would  be  0(A9   ).     If  trapezoidal  integration  were  used  with  a  step 

3 
length  of  5  A9,   the  error  at  each  step  would  be  0(125  A0   )  which  is 

significantly  less  than  the  rectangular  error.     This  increase  in  step 
size  permits  reduction  of  computer  speed  by  about  a  factor  of  two. 
However,   there  is  an  irreducible  uncertainty  in  the  gyro  information 
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2 
of  O(A0   )  in  each  pulse,   so  that  the  apparent  increase  in  accuracy  may 

not,  in  fact,  be  attained.     As  is  shown  in  Chapter  III,  Section  3.B, 

rnTs 
each  output  pulse  of  a  gyro  is  \co-jdt      ,   the  integral  of  one  component 

H[n-1)TS 
of  angular  velocity  resolved  into  body  coordinates.     Since  addition  of 

rotations  is  not  commutative, 


nT„  nTQ  nT  nT 

r*  r*  r*  /■»     s 

\ wdt  £T     \wjdt+j      \u2dt  +  k      \wgdt  (5-8) 

c(n-l)Tg      (n-l)Ts       (n-l)Tg         (n-l)Tg 

in  general.     In  one  sampling  period,   the  angular  velocity  could  be 

equal  about  all  three  axes,  or  there  could  be  successive  rotations 

about  all  axes  in  any  order;  each  input  gives  the  same  output. 

Eq.   (5-8)  adds  these  rotations  as  though  they  were  vectors.     The 

2 
error  thus  made  is  O(A0   )  (Ref.   34,   p.    97).     The  conclusion  to  be 

drawn,   therefore,   is  that  rectangular  integration  of  Eq.  (5-1)  is  the 

most  suitable  algorithm  because  the  angular  information  is  available 

in  rotational  increments.     Note  that  this  comment  applies  to  any 

digital  computation  where  the  rotations  are  handled  in  this  way;  it 

is  not  restricted  to  delta-modulated  gyro  information.     In  special 

cases,   however,   the  use  of  higher  order  integration  rules  to  reduce 

truncation  error  may  be  of  some  advantage. 

The  equations  to  be  instrumented  to  generate  the  direction 

cosines  are  three  sets  of  the  form 

AC,   =  C-0A9Q  -  C,A99 
ll  \l       3  i3       2. 

ACi2  =  Gi3A9l   "CilA93  <5"9) 

ACi3  =  CnAe2  -  Ci2A91 

5.C    Digital  Differential  Analyzers 


The  computer  chosen  to  generate  and  perform  the  coordinate 
transformation  must  be  fast  and  accurate.  The  accuracy  require- 
ment,  with  present  equipment,   eliminates  analog  computation  from 
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Register  Arrangement 


AY 
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AZ 


(b) 

Schematic 
Fig.    5-1    DDA    Integrators 
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consideration;    the  speed  requirement,  based  on  the  necessity  for 
processing  each  pulse  as  it  comes  from  the  gyro,  calls  for  a  com- 
puter whose  cycle  time  is  no  greater  than  the  sampling  period  of  the 
gyro  loop.     The  sampling  period,   as  determined  in  Chapter  III, 
Sec.   3.B.4,  is  the  quotient  of  maximum  angular  velocity  and  angular 
increment. 

A  very  fast  and  compact  type  of  digital  computer,  and  one 
particularly  well  suited  to  the  solution  of  Eq.  (5-9),  is  the  digital 
differential  analyzer  (DDA)  (Refs.   43-45,   53). 

The  DDA  has  been  described  as  an  analog  computer  in  digital 
form;   it  has  the  advantage  found  in  the  Bush  differential  analyzer 
(Ref.   42)  of  not  requiring  time  to  be  the  independent  variable  of 
integration.     The  basic  element  of  DDA  computers  is  the  DDA 
integrator,  usually  shown  diagrammatically  as  in  Fig.    5-l(a).    The 
blocks  labeled  R  and  Y  are  storage  registers.     The  line  labelled 
AY  carries  pulses  representing  a  change  (positive  or  negative 
according  to  the  sign  associated  with  the  pulse)  in  the  least  significant 
bit  in  Y.     The  AX  line  carries  pulses  representing  equal  increments 
in  X;   when  a  AX  pulse  occurs,  the  contents  of  the  Y  register  are 
added  into  the  R  register  or  subtracted  from  the  R  register,  depending 
on  the  sign  associated  with  the  AX  pulse.     The  contents  of  the  Y  reg- 
ister are  left  unchanged  by  this  operation.     The  AZ  line  carries  over- 
flow pulses  from  the  R  register  representing  a  change  in  the  least 
significant  bit  in  Z.     The  AZ  pulses  may  be  used  as  AX  or  AY  pulses 
in  another  integrator.     The  name  "integrator"  is  applied  to  this  unit 
because  it  computes 

AZ  =  YAX  (5-10) 

This  is  an  increment  of  rectangular  integration  of  Y  with  respect  to 
X  [See  Eq.   (5-5)]. 

Each  of  the  variables  X,  Y,  and  Z  represents  a  physical 
quantity  to  some  scale.     If  the  Y  register  is  n  bits  long,   so  is  the  R 
register.     The  least  significant  bit  in  the  R  register,  AR,   is  related 
to  AY  and  AX  by 
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Fig.    5-2    Arrangement  for  DDA  Solution  of  Equation  (5-9) 


81 


AR  =  AXAY  (5-11) 

The  overflow  logic  for  the  R  register  causes  an  overflow  pulse  when 
the  content  of  the  R  register  exceeds  AZ.     The  simplest  kind  of 
overflow  logic  scales  AZ  so  that 

AZ  =  2nAR  (5-12) 

Thus,   the  AZ  overflow  occurs  when  the  R  register  fills  up.     AZ  is 
then  related  to  AX  and  AY  by 

AZ  =  2nAXAY  (5-13) 

The  usual  schematic  representation  for  a  DDA  integrator  is 
shown  in  Fig.    5-1  (b). 

The  solution  of  Eq.   (5-9)  requires  six  incremental  integrators 
for  each  value  of  i  arranged  as  in  Fig.    5-2.     Note  that  no  order  of 
performing  the  additions  is  specified. 

There  is  feedback  in  this  arrangement;   hence  there  are  scaling 
consistency  equations  to  be  satisfied,  just  as  in  the  case  of  analog 
computation.     The  integrators  are  identical.;    thus  consideration  of  the 
loop  with  the  A0g  integrators  leads  to 

A93ACu2n  =  A0g  ACi222n  =  ACi2  (5-14) 

or 

A03  =  2"n  (5-15) 

by  repeated  application  of  Eq.   (5-13).     A  similar  result  holds  for 
AG,  and  AGg.     Thus  for  simplest  overflow  logic,  AG  should  be  chosen 
as  a  negative  integral  power  of  2.     The  value  which  AG  must  not  exceed 
may  be  found  from  system  specifications,   and  n,   the  length  of  the  Y 
and  R  registers,   is  determined.     Requiring  that  the  scaling  of  the 
velocity  increments  be  the  same  in  both  coordinate  frames  yields 
AC  =  2'n. 
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5.  D    Analysis 

With   the  form  and  function  of  the  transformation  computer  de- 
termined, an  analysis  of  its  operation  must  be  made.     Three  points 
are  usually  considered  in  the  analysis  of  a  computer:    dominant  mode, 
truncation  error,   and  round-off  or  quantization  error.     These  will  be 
taken  up  in  turn  after  a  frequency  domain  model  for  the  computer  has 
been  developed. 

5.  D.  1    Frequency  Domain  Model 

Computers  operating  in  real  time  have  been  characterized  as 
sampled-data  systems  and  transfer  functions  derived  for  them 
(Ref.    51).     Since  the  equations  that  the  computer  is  solving  are  dif- 
ference equations,   they  may  be  transformed  in  a  way  similar  to  that 
used  in  transforming  time  sequences  or  samples  of  continuous  func- 
tions (Ref.   61,  Ch.   1). 

The  equations  describing  the  computer  operation  with 

u-,   =  u>2  ~  °>   u3  ^  °  are 

Cn(n)  =  Cn(n-1)  +  A93Ci2(n-l) 

(5-16) 
Ci2(n)  =  Ci2(n-1)-  A93Cil(n) 


AG 

-s 


Let  Zo  =  e      Wo     where  s  is  the  Laplace  complex  frequency 

variable.     Then,  if 


C    (z)  =  ZC:(n)zn  (5-17) 

J  J 


Eq.    (5-16)  may  be  transformed  into 


A9o  z~ 

Cil(z3>=r-^  Ci2(z3> 
L"  Z3 

AG 

Ci2<z3>  =  -  — Cil(z3> 
l-z3 
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(5-18) 


If  u,  =  2u>2,  u>o  =  0,   the  transformed  set  of  equations  is  (see  App.  C, 
Sec.  C.   l.b) 

A9z    l+A9  z 

1  -  z2 
Cn(z)  =  -  i  Ci3(z) 

(1  -  z)(l  +  z)   1  +  AG2  1  +  ,Z_ 

Ci2(z)=^^-  Ci3(z)  (5-19) 

1  -  z 

C.3(z)  =  _*i_  C„  z  -  -*L  C.2(z) 
16  1  -  z2      U  1  -  z      lZ 

A9  0  0  A9  AG 

-s  —  ,2         -s    2  — -         -s 


where  z  =  e       <*>      and  z=e  w,=e      w  o    • 

This  set  of  equations  may  be  well  approximated  by 

C.l(z)  =  --**L  c„(.) 

1  -  Z 

Ci2(z)  =  *?*-  Ci3(z)  (5-20) 

1  -  z 

Ci3<z)=-J^Cil(z).^ci2(Z) 

1-z2  l'z 

AG  AG 

If  the  variables  z,  and  z„  are  defined  by  z,   =  e         u,    and  z2  =  e        u>2  , 

then  Eq.   (5-20)  may  be  written 

C    (.)  -  -  -*?*-  Ci3(z) 
1  "  z2 

Ci2(z)  =_Mz_c.3(z)  (5-21) 

l-zL 

Ci3<z>  "  ~  Cil<z>  "  7^-  Ci2<z> 
1  -  z2  1  -  zx 
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where  z  represents  the  delay  between  the  cycling  of  integrators.     If  the 
components  of  angular  velocity  are  related  by 

kid-,   =jfuo  =  mug  =  w  (5-22) 

AG 
-s 

where  k,  i,  and  m  are  integers,   and  if  z.  =  e      ui   ,   then  the  z.'  s  are 
related  by  /" 

zfm  =  z2km  =  z3k^=  zk*m  (5-23) 

Here  again  z  represents  the  delay  between  the  cycling  of  the  integrators. 
With  the  above  relations  between  angular  velocity  components  satisfied 
the  set  of  equations  [cf.    Eq.   (5-9)] 

Cn(n)  =  Cn(n-1)  +  A03Ci2(n-l)  -  AQgC.gdi-l) 

C.2(n)      Ci2(n-1)  +  A01Ci3(n-l)  -  A03Cil(n)  (5-24) 

Ci3(n)      Ci3(n-1)  +  A02Cu(n)      -  AQjC.^n) 

may  be  approximately  transformed  into 

1  l-z3      l2  l-z2 

d,fe)  =  —  C.,(z)  — —  C,,(z)  (5-25) 

1-z,  l-z3     u 

Ci3W%^LCil(z)--^C12<,) 

i-z2  l-z1 

With  w.   =  2u2,  u~  =  0,   this  set  is  identical  with  Eq.   (5-21),   which 
is  a  good  approximation  to  Eq.   (5-19).     Similarly  with  y,  =  u2  =  0, 
u„  ^  0,    Eq.   (5-21)  becomes  Eq.   (5-18).     Note  that  if  y^  =  0,    the  terms 

in  z.  are  taken  as  zero. 

1 
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Eq.    (5-25)  may  be  represented  by  a  signal  flow  graph  as  shown 
in  Fig.    5-3.     When  the  angular  velocity  of  the  system  is  specified,   the 
substitutions  indicated  by  Eq.   (5-22)  and  Eq.   (5-23)  may  be  made. 
For  example,   with  u,   =  2u>2  =  3cjo,  k  =  I,  jf  =  2,  m  -  3,   and  the  appro- 
priate signal  flow  graph  is  shown  in  Fig.    5-4. 

The  use  of  Mason'  s  method  for  reducing    signal  flow  graphs 
(Ref.   60,  Ch.   4)  shows  that  the  characteristic  equation  for  the  com- 
puter under  these  conditions  is 

1  +    Aq2z      +    Aq2z       +    Aq2z ^1 =0         (5-26) 

o  n     2  q     2  9  "i 

(1-zr       (l-zJ)         (l-zd)         (l-z^d-z"3) 

and,   in  general,   the  reduction  of  the  signal  flow  graph  in  Fig.    5-3 
leads  to  the  characteristic  equation  for  the  computer 

1  +    Aq2z     +    aq2z        +    A82z       _      A83z(l  -  z) 

(1-z/      (l-z2)2     (l-z3)2     a-Zl)(i-z2)(i-z3)  =  0 

(5-27) 

With  the  node-to-node  transfer  functions  determined,   the  remainder  of 
the  analysis  can  proceed. 

5.  D.  2    Determination  of  Dominant  Mode 

The  various  normal  modes   of  the  transformation  computer  are 
determined  from  the  roots  of  the  characteristic  equation.     Each    root  £ 
leads  to  a  time  sequence  of  the  form  t,       with  coefficients  chosen  to  fit 
initial  conditions.     If  t,  is  complex,   its  conjugate  will  also  be  a  root, 
and  an  oscillatory  term  will  be  present.     If  It,  I  >  1  then  the  mode  will 
decay,  while  if  I  £,  I  <  1  there  will  be  exponential  growth. 

The  initial  conditions  as  the  varicus  nodes  appear  as  step  inputs; 
for  this  reason,   the  roots  near  z  =  1  will  provide  the  dominant  modes, 
while  other  modes,   which  occur  near  the  various  n      roots  of  unity, 
will  have  magnitude  A9.     Thus  the  dominant  mode  of  the  computer 
comes  from  the  roots  near  z  =  +1. 


86 


Fig.    5-3    Signal  Flow  Graph 
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Fig.    5-5    Analog  Computer  Flow  Graph 
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The  magnitude  of  the  roots  and  their  associated  modes  may  be 

2 
found  by  root-locus  techniques.     Since  the  loop  gain,  A0   ,   is  so 

small,  linear  approximations  must  be  used  in  conjunction  with  the 

graphical  technique.  . 

In  Appendix  C,  the  analysis  of  two  situations  is  carried  out 

using  root-locus  techniques.     For  the  case  of  angular  velocity  ratio 

two  to  one,   there  are  extraneous  roots  near  the  square  roots  of  unity 

not  +1.     In  the  case  of  equal  rates  on  all  axes,   there  are  only  roots 

near  +1. 

5.  D.  3    Truncation  Error 

Truncation  error  is  the  error  made  by  not  considering  higher 
derivatives  or  differences  when  computing  the  solution  of  a  differ- 
ential equation.     It  is  the  difference  between  the  exact  solution  and 
the  computed  solution  using  computer  registers  with  infinite  length. 

Assumption  of  the  same  angular  velocity  relationships  as  above 
[Eq.   (5-22)]  with  the  added  restriction  that  the  angular  velocity  be 
constant,  permits  the  solution  of  Eq.   (5-1).     With  Eq.   (5-22)  satisfied, 
Eq.  (5-1)  becomes  three  sets  of  the  form 

r>      -  c       w       n       u 
Hi  "  H2  m  "     i3  J 


u) 


io 


H2  "  S3  k  '  Hi  m 


(5-28) 


S3  "  Si  |  "      i2  k 


Applying  a  Laplace  transformation  to  Eq.   (5-28)  yields 


Cu(s)=I 

s 


'-Ci2<s>-T-Ci3(s>  +  Cil<0)  ] 


m 


Ci2(s)=I 


H  ci3(s)  -  Ji  C.,(s)  +  ci2(o) 

k  m 


(5-29) 


Ci3(s)  -1 

s 


7Cil(s)-fCi2(s)+Ci3(0) 

Z  k 
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where  the  C--(O)  are  initial  conditions.     Converting  this  set  of  equations 

J 

to  a  signal  flow  graph  leads  to  Fig.    5-5.     The  initial  conditions  are  taken 
as  inputs  at  appropriate  nodes.     Flow  graph  reduction  may  be  used  to 
derive  expressions  for  the  signals  at  each  node.     For  example, 


Cn(0) 


Cn(s)  = 


s2+-^ 
.  k2 


2i 


+  u 


C.3(0)      ci2(o)i 


m 


s  +u 


Cj3<°>  |  C12(0) 


km 


k* 


(k     +1     +  m   )u   s 


+  1 


_(k2  +  l2  +  m2)w2 


(5-30) 


These  expressions  can  be  transformed  into  the  time  domain  and  the 
exact  solution  to  the  differential  equations  found. 

The  solution  to  the  difference  equations  instrumented  by  the  com- 
puter may  be  found  by  taking  the  z-transform  of  the  solution,   as  deter- 
mined by  root-locus  methods,  and  writing  down  the  corresponding 
time  sequences.     There  will  be  truncation  error  from  three  causes; 
the  roots  near  unity  may  not  have  the  proper  argument,   they  may  not 
be  on  the  unit  circle,  and  there  may  be  other  roots.     The  error  in 
argument  of  the  root  near  unity  makes  the  oscillations  of  the  direction 
cosines  the  wrong  frequency.     If  these  roots  are  not  quite  on  the  unit 
circle,   the  oscillations  will  not  be  of  constant  amplitude  as  they 
should.     The  other  roots  may  lead  to  exponential  or  oscillatory 
terms;  in  any  case  the  size  of  the  corresponding  mode  may  be  deter- 
mined. 

When  finding  the  z-transform  of  the  computer  output,   the  last 

root-locus  construction   to  find  the  roots  of  the  computer  character- 

2  2 

is  tic  equation,   deals  with  a  loop  that  has  a  gain  of  A0    .     Since  A  9 

is  very  small,   the  construction  cannot  be  carried  out  exactly,  but 
linear  approximations  may  be  made  analytically.     These  approxima- 
tions are  sufficient  to  make  the  truncation  errors  predicted  inexact. 
However,   the  order  of  the  errors  predicted  is  correct  and  the  growth 
rate  is  correct,  but  an  exact  error  analysis  is  not  provided.     A  more 
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nearly  exact  error  analysis  may  be  had  by  finding  the  roots  of  the 
characteristic  equation  exactly,   and  finding  the  magnitudes  of  the 
modes  by  root-locus  techniques,  but  the  present  analysis  provides 
enough  information  to  form  a  basis  for  engineering  decisions. 

5.  D.  4    Quantization  Error 

Quantization,   or  round-off  as  it  is  sometimes  termed  by 
numerical  analysts,   occurs  whenever  a  number  is  assigned  to  a 
physical  variable.     Through  the  quantization  process  the  value  of 
the  variable  is  represented  as  the  nearest  integral  number  of  the 
basic  units  in  which  the  variable  is  being  measured.     If  this  basic 
unit  be  q,   the  quantization  process  may  be  represented  by  the  input- 
output  relationship  shown  in  Fig.   5-6, 

While  this  is  the  basic  quantization  process,   the  overflow  logic 
of  many  computers  causes  hysteresis  in  the  quantization  process. 
This  hysteretic  quantization  has  the  characteristic  shown  in  Fig. 
5-7  (a).     The  hysteresis  may  be  represented  by  feedback  around  the 
basic  quantizer  as  shown  in  Fig.    5-7(b).     Note  that  the  quantizer  has 
grain  size  2q.     The  z  in  the  feedback  loop  represents  a  delay  of  one 
cycle  time  of  the  computer. 

Since  quantization  distorts  the  signal,   although  in  a  deter- 
ministic way,   the  distortion  may  be  thought  of  as  noise.     Widrow 
(Ref.    55)  shows  that  if  the  amplitude  of  the  input  to  the  quantizer 
ranges  over  more  than  three  quantization  boxes,   the  quantization 
noise,   to  a  good  approximation,  may  be  treated  as  an  additive  first 
order  random  process,  q(t).     The  probability  density  distribution 
of  the  quantization  noise  is  shown  in  Fig.    5-8.     The  expressions 
for  the  autocorrelation  function  and  the  power  spectral  density  of 
quantization  noise  are 

2 


0nn(T)  =  S_  6(t) 

(5-31) 


qq  12 


2 

$      (z)  =  q— 
qqV  12 
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Fig.   5-8    Probability  Density  Distribution  of  Quantization 
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Fig.   5-9    Representation  of  a  Quantizer  by  Independent  Additive  Noise 
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Fig.    5-9  shows  the  representation  of  a  quantizer  as  independent  additive 
noise. 

Although  one  attempt  was  made  to  apply  Widrow'  s  theory  to 
DDA1  s  in  real  time  (Refs.   49,   50),   some  difficulty  was  encountered 
in  applying  it  to  interconnected  DDA1  s. 

The  appropriate  signal-flow  graph  for  the  computer  with  non- 
hysteretic  quantization  is  shown  in  Fig.    5-10,   and  Fig.    5-11  shows  the 
quantizers  replaced  by  their  equivalent  gains  and  additive  noise. 

If  hysteretic  quantization  is  used,   a  self-loop  with  gain  +z/2 
must  be  included  at  each  node.     With  integrators  cycling  at  different 
rates,   the  z  used  in  the  self-loop  must  correspond  to  the  fastest  basic 
rate:     the  z  of  Eq.   (5-23)  must  be  used. 

Care  must  be  taken  to  insure  that  a  zero  value  of  the  variable 
actually  corresponds  to  the  zero  of  the  computer.     In  many  cases  a 
bias  of  ^  q  is  introduced  because  of  an  ambiguity  of  initial  condi- 
tions (Ref.    49,   p.    51). 

Using  the  statistical  model  of  quantization  and  the   sampled - 
data  model  of  the  transformation  computer,   an  expression  for  the 
power  spectral  density  of  quantization  noise  may  be  developed.     For 
example,  with  u,   =  2^^  =  3uo,   the  power  spectral  density  of  quan- 
tization in  Cm  is 


*n(Z)  = 


2 
A0.  z 


+ 


A0z 


2    2 

A9z         A0    % 

+ 


+- 


A02z 


(1-z)         l-z°       1-z' 


(l-z)(l-z3)       (l-z)(l-z2) 


1+A92z     +A92z 


+ 


A02z 


A93z 


(1-z)' 


(1-z2)2      (1-z3)2 


(l-z2)(l-z3) 


A9 


12 


(5-32) 


with  similar  expressions  for  power  spectral  densities  at  other  nodes. 
Note  that  for  the  case  of  the  computer,   the  quantization  grain  size 
is  A0. 

The  mean  square  quantization  error  may  be  obtained  by 
forming  the  sum  of  the  squares  of  the  sequence  represented  by  the 
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transfer  function  in  Eq.   (5-32)  and  multiplying  this  sum  by  the  mean 
square  of  the  quantization  noise  (Ref.   61,   Ch.   7).     Since  in  the  case 
of  this  computer  the  roots  of  the  characteristic  equation  have  magni- 
tude close  to  unity,   this  mean  square  value  may  become  very  large. 
Thus,  it  is  appropriate  to  consider  only  a  finite  number  of  terms 
in  the  sum-of-squares  series  corresponding  to  a  finite  number  of 
steps  in  the  computation.     With  the  series  truncated  after  a  finite 
number  of  terms,   an  estimate  of  the  error  may  be  obtained.     An 
example  is  shown  in  Appendix  C,  Section  C.  3.     It  is  found  that 
quantization  error  is  much  smaller  than  the  truncation  error;  the 
number  of  steps  required  to  make  the  rms  quantization  error  equal 
to  AG  is  0(A9"2). 

In  the  same  section  (Section  C.  3),   an  expression  for 
hysteretic  quantization  error  is  derived.     It  is  found  that  the  mean 
square  value  of  hysteretic  quantization  error  is  greater  than  that 
of  non-hysteretic  quantization  by  a  factor  of  four  because  of  the 
larger  quantization  box.     For  this  reason,   non-hysteretic  overflow 
logic  should  be  used  in  the  transformation  computer. 

5.  D.  5    Reversibility 

Up  to  this  point,   the  analysis  has  assumed  that  all  the  AG 
pulses  have  been  of  the  same  sign.     Of  course,   this  will  never  be 
the  case,  and  some  examination  must  be  made  of  the  effects  of 
negative  pulses. 

In  the  sine-cosine  generators  shown  in  Fig.    5-12,   first  one 
variable  is  updated,   and  then  the  other  variable  is  updated.     Thus, 
Fig.    5-12(a)  implies  the  difference  equations 

cos  (nA8)  =  cos  [(n-l)AG]  -  AG  sin  [(n-l)AG] 

(5-33) 
sin   (nAG)  =   sin  [(n-l)AG]  +  AG  cos    (nAG) 

If  a  positive  pulse  is  followed  by  a  negative  pulse,   there  should  result 

cos  [(n+l)AG]  =  cos  [(n-l)AG]  (5-34) 
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Fig.    5-12    Reversibility  of  Sine-Cosine  Generation 

The  use  of  Eq.   (5-33)  for  processing  both  pulses  yields 

cos  [(n+l)A6]  =  cos  [(n-l)A9]  +  AG2  cos  [(n-l)A0]  (5-35) 

If  positive  pulses  are  processed  according  to  Fig.    5-12(a),   however, 

and  negative  pulses  according  to  Fig.    5-12(b),   Eq.   (5-34)  is  satisfied, 

the  computation  is  completely  reversible,  and  the  truncation  and 

quantization  errors  for  the  two  steps  cancel  each  other. 

This  idea  may  be  extended  to  the  direction  cosine  computation. 

In  Fig.    5-3,   all  the  pulses  are  assumed  positive,  and  the  order  of 

processing  the  pulses  updates  C.,  first,   then  C-2,   then  C.~.     Since 

the  computation  is  made  up  of  a  number  of  sine-cosine  generators 

interconnected,   the  process  can  be  reversed  by  moving  the  z  in  the 

numerator  of  each  loop  transfer  function  to  the -Other  branch;    this 

may  be  accomplished  by  updating  the  variables  in  reverse  order 

for  negative  pulses. 

There  are  three  pulses,   each  with  two  possible  signs,   which 

3 
drive  the  computer.     This  gives  2     =8  possible  combinations. 
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There  are  only  three  direction  cosines  in  a  set;    the  number  of  possible 
orders  of  updating  them  is  3!    =6.     Thus  one  pair  of  combinations 
cannot  be  reversed.     Using  the  updating  sequence  1,2,3,  for  all  posi- 
tive pulses,   the    computer   cannoti reverse.  A92  +  ;  A9,-, A9„-; 
and  A92-;  A9,+,  A0o+.     With  another  sequence  representing  all  pulses 
positive,   another  pair  could  not  be  reversed;   however, any  combina- 
tion of  two,   as  well  as  any  single  pulse,   can  be  completely  reversed. 

The  non-reversibility  of  this  computation  causes  a  drift  of 
2 
0(A6   )  for  each  step  reversed.     Unless  the  computation  is  operating 

with  marginal  accuracy,   the  added  equipment  required  to  change  the 

order  of  updating  the  variables  is  probably  too  expensive,   in  terms 

of  weight,   power,   and  cost,   to  pay  for  the  additional  precision 

achieved. 

5.  E    Processing  AV  Pulses 

At  each  sampling  instant  in  the  accelerometer  loop  a  AV 
vector,   resolved  into  body-fixed  components, appears  at  the  accel- 
erometer outputs.     If  a  two-level  relay  is  used  in  the  loops,   there 
will  be  pulses  at  each  loop  output  at  every  sampling  instant;  if  a 
three-level  relay  is  used,  output  pulses  will  occur  only  at  certain 


times.     The  AV  vector  may  be  resolved  into  inertial  coordinates 
by  processing  it  with  the  direction  cosine  matrix: 

AVj  =  CAVb  (5-36) 

Since  the  magnitudes  of  the  AV  pulses  are  the  same,    Eq.   (5-36)  may 
be  realized  by  the  arrangement  shown  in  Fig.    5-13.     The  AV,  pulses 
will  be  accumulated  in  VT  registers.     Errors  also  arise  in  this 
computation,  and  they  may  be  determined  by  the  methods  outlined 
above. 

5.  F    Summary 

In  this  chapter,    the  design  and  analysis  of  the  transformation 
computer  was  completed.     Various  quadrature  rules  were  examined, 
and  the  rectangular  rule  was  chosen  as  the  most  appropriate  for 
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integrating  the  direction  cosines  [Eq.   (5-9)]  .     With  the  angular 
information  available  in  A9  pulses,   digital  differential  analyzers 
seem  the  most  suitable  form  of  computer  for  the  computation 
(Fig.    5-2).     A  frequency  domain  model  for  the  computer  was  de- 
veloped (Fig.    5-3)  and  analyzed.     The  truncation  error  was  estimated 
by  developing  an  exact  solution  for  the  differential  equation,   and 
comparing  it  with  the  computer  solution  as  indicated  by  root-locus 
plots.     Quantization  was  treated  statistically  and  it  was  found  that 

hysteretic  quantization  caused  a  greater  error  than  non-hysteretic. 

2 
The  integrators  were  shown  to  drift  by  O(A0   )  for  each  computation 


step  reversed.     Finally,   a  method  for  transforming  the  AV  vector 
was  established. 
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CHAPTER  VI 
SYSTEM  DESIGN 


6.  A    General  Considerations 

The  previous  chapters  have  described  and  analyzed  the  component 
parts  of  the  system;  the  design  parameters  must  now  be  chosen  so 
that  the  system  operates  in  an  optimum  manner  to  carry  out  the 
mission  assigned.     Two  decisions  may  be  made  without  reference  to 
the  mission  which  the  system  will  perform:  whether  the  instruments 
should  limit  cycle,  and  how  the  clocks  in  the  gyro  and  accelerometer 
loop  should  be  related  to  each  other. 
6.A.1    Limit -Cycling  Instruments 

In  Chapter  V,    Section  JxfJl  5,   it  is  shown  that  for  each  step  reversed 

in  the  direction -cosine  computation,   the  integrators  drift  by  an  amount 

2 
0(A#   ).     For  this  reason,   it  is  desirable  not  to  have  the  gyros  limit 

cycling.     With  a  three -level  relay  in  the  gyro  loop,   the  torquer  current 

may  be  turned  off  during  the  quiescent  periods,   thus  saving  a  great 

deal  of  power.     Turning  the  current  on  and  off,  however,   makes 

stabilization  of  the  current  more  uncertain. 

Since  the  velocity  pulses  are  being  processed  by  DDA's,  there 
may  be  a  drift  of  O(A0AV)  in  accumulated  velocity  for  each  reversal. 
Thus,  the  arguments  above  indicate  the  desirability  of  eliminating 
the  limit  cycles  in  the  accelerometer  loops. 
6.  A.  2    Instrument  Clock  Relationships 

In  general,  the  clock  rate  of  the  accelerometers  will  be  slower 
than  that  of  the  gyros;  thus  the  direction  cosine  may  change  many 
times  between  A  V  pulses.     Processing  of  the  A  V  pulses  must 
occur  between  changes  of  direction  cosines;  this  requirement  may 


be  satisfied  either  by  providing  buffer  storage  for  the  A  V  pulses,   or 

by  choosing  the  clock  frequencies  and  phasing  of  the  gyro  and  accelerometer 

loops  so  that  the  A  V  pulses  will  always  occur  at  the  same  point  in  the 

computer  cycle.     The  latter  choice  is  more  attractive,  since  it 

eliminates  the  need  for  buffer  storage  and  the  consequent  delay  in 

processing  information. 

6.B    Basic  Design  Parameters 

The  quantities  needed  to  design  a  guidance  system  are  the  inputs, 
or  maneuvers  it  will  experience,   the  allowable  error  at  the  target,   and 
the  time  of  flight.     The  input  specifications  must  include  the  maximum 
acceleration,  velocity,   and  angular  velocity  which  the  guidance  system 
is  expected  to  undergo.     From  these  data,  the  basic  instrument  and 
computer  parameters  may  be  determined. 

With  increment  size  and  clock  frequency  specified,  Chapter  III  and 
Chapter  V  indicate  methods  for  determining  the  remaining  system 
parameters.     In  the  two  sections  that  follow,   computations  are  carried 
out  to  determine  clock  frequency  and  increment  size  for  typical  aircraft 
and  missile  missions.     Numbers  representing  current  performance  are 
substituted  into  the  expressions  for  increment  size  and  clock  frequency 
to  provide  a  notion  of  the  order  of  the  magnitudes  involved. 
6.C    Airplane  Flights 

Airplane  guidance  systems  are  position  control  systems;  their 
task  is  to  guide  the  aircraft  to  a  specified  position.     Thus,   errors  in 
position  at  the  target  are  errors  in  the  controlled  quantity.     To  determine 
specifications  for  the  instrument  increment  sizes,   it  is  appropriate  to 
analyze  a  "worst  possible"  case,    and  equate  the  error  made  to  the 
error  specification. 

An  airplane  is  essentially  a  constant  speed  vehicle;  thus,  the  worst 
possible  case  for  the  generation  of  error  is  a  flight  at  maximum  velocity 

with  the  direction  of  the  velocity  in  error  by-^-  A0,.  and  the  magnitude  of 

1 
the  velocity  in  error  by  -s-AV.     If  the  allowable  error  in  position  E     is 

the  same  for  gyro  and  accelerometer  quantization,   the  appropriate 

calculations  are 
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1  E 

2  max  Tf 

and 


*  Xf  (6-1) 

where  T,.  is  the  time  of  flight.     Thus  the  required  increment  sizes  are 

2E 

**-- — rf 

max    f 
and  2E  (6-2) 


AV  = 


Tf 


Note  that  Ad  should  be  chosen  as  a  power  of  2.     The  clock  frequencies 
are 


fg  = 

W 
max 

A0 

f    = 

a 

a 
max 

AV 

and  (6-3) 


Specifications  which  represent  a  typical  high  performance  airplane  are 

a  =  5g        =  161  ft/sec 

max  &  ' 


Substitution  gives 


W  =3  rad/sec 

max  ' 


V  =  600  kt 

max 

E  /T.       =  .1  kt 


(6-4) 


A0  =  333  jurad 
AV=  .  2  kt  =  .333  fps. 
Since  A0  should  be  chosen  as  a  negative  integral  power  of  two,   take 

A0  =  2~12  =  244  jurad 

Then  the  clock  frequencies  are 

f     =  12.3  kc 
g 

f     =  482  cps 
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Since  it  is  desirable  to  have  f     a  subharmonic  of  f     set, 

a  g 

P  12300        Aao 

f     =     — rrz—  =  492  cps 
a  25  r 

The  complete  set  of  design  parameters  is 

A0     =  244  jurad. 

f        =  12.  3  kc 
g 

AV    =  .333  fps 

f        =  492  cps 
a 


(6-5) 


The  period  of  the  gyro  loop  is  81.  3  microseconds.     With  18  additions 
to  be  made  in  the  transformation  computer  in  this  time,   an  add  time  of 
4.  5  microseconds  per  word  is  required.     Since  a  first-rate  arithmetic 
computer  with  central  control  would  be  required  in  toto    to  accomplish 
this  task,   a  DDA  is  clearly  the  best  machine  for  this  job.     It  will  be 
smaller,   faster,   and  easier  to  build. 

For  serial  addition  of  a  word  with  thirteen  bits  (twelve  plus  one  for 
sign)  in  four  and  a  half  microseconds,   a  bit  rate  of  about  three  megacycles 
is  required.     This  performance  is  available  in  present  DDA's. 

With  increment  sizes  and  clock  frequencies  fixed,  the  remaining 
parameters  may  be  quickly  determined. 
6.D    Missile  Flights 

In  contrast  to  an  airplane  guidance  system,   a  missile  guidance 
system  is  essentially  a  velocity  control  system,   and  a  specified  error 
in  point  of  impact  must  be  translated  into  an  error  in  velocity  at  thrust 
termination.     Relations  between  increment  sizes  and  velocity  error  at 
cutoff  are  determined  and  translated  down  range  using  linear  perturbation 
techniques  to  find  impact  error. 

The  relationship  between  trajectory  parameters  and  point  of  impact 
is  called  the  hit  equation  (Ref.  65).     The  perturbation  factors  are  called 
range  derivatives  and  are  obtained  by  differentiating  the  hit  equation. 
The  velocity  range  derivative,   6R/6V,   is  the  one  of  interest  in  the 
following  computation. 

The  worst  possible  case  to  be  analyzed  here  is  where  the  velocity 
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is  in  error  by  -*-  AV  in  the  direction  of  maximum  range  derivative,   and 
where  the  error  due  to  gyro  quantization  is  a  velocity  error  of  -9— V     A9 
in  the  direction  of  maximum  range  derivative.     The  origin  of  the  first 
error  is  rather  obvious;  the  origin  of  the  second  is  not. 

Suppose  that  as  each  AV  pulse  is  processed,  the  direction  cosines 


by  -5-   AC    =   p-A9,  then  the  error  at  each  step  is  5-  A&AV. 
Summing  this  over  the  total  velocity  accumulation  yields 


are  in  error  *.      2    —^         9 


E(V)  =  S  i  A0AV    =  ^-A9  2AV    =  4a0V 

a  a  £  CO 

(6-6) 
The  appropriate  calculations  for  increment  size  are 


and  (6-7) 


or 


Iaov 

6R    _ 

r  Tr         ~          ill 

co    6V            p 

i-AV 

6R    _    „ 

6V    ■    Ep 

2E 

A3    = 

-,.       P 
V      6R 

C°6V 

and 


(6-8) 


AV    = 


2E 
P 


6R 

w 


The  clock  frequencies  are  then 

W 
f  max 


g         A9 

a 

f    _       max 


(6-9) 


"a         AV 

Typical  specifications  for  an  1800  nautical  mile  missile  are 

a  =  10g=322  ft/sec2 

max  &  ' 

W  =    .  2   rad/sec 

max  ' 
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V         =    20,000fps 

E       =    300  ft  (6-10) 

P 

Substitution  gives 

A3  =  33.  3  Mradj 

AV=  .333  fps 
Since  AB  should  be  a  negative  integral  power  of  two,  take 

A9  =  2"15  =  30.5Mrad." 
The  clock  frequencies  are  then 


and 


f       =6.55  kc 
g 


f     =  484  cps 


Since  the  accelerometer  clock  frequency  should  be  a  subharmonic 
of  the  gyro  clock  frequency,  take 

~  6550  r-r,  A 

f     =  — TZ~~=    504  cps 
a  13  r 

The  complete  set  of  parameters  is 

A6     =    30.5  jurad 

f       =6.55  kc 
g 

AV    =     .333  fps  (6-11) 

f       =    504  cps 
As  in  the  case  of  the  airplane  mission,   the  computer  must  be 
quite  fast,   requiring  an  add  time  of  about  eight  and  a  half  microseconds 
per  word.    Again,   the  logical  choice  for  the  computer     is  the  DDA, 
since  an  airborne  computer  must  be  compact.     Further,   the  advantage 
of  flexibility  associated  with  an  arithmetic  machine  would  not  be 
realized  because  of  the  fast  add  time  required.     For  a  serial  DDA  to 
perform  eighteen  additions  of  words  sixteen  bits  long  in  the  computer 
cycle  time,   a  bit  rate  of  about  two  megacycles  is  required.     This 
performance  is  available  in  present  DDA's. 
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CHAPTER  VII 
SUMMARY  AND  CONCLUSIONS;  SUGGESTIONS  FOR 


FURTHER  STUDY  AND  DEVELOPMENT 

7.  A      Summary  and  Conclusions 

This  work  has  been  concerned  with  a  broad  theoretical  analysis  of 
gimballess  inertial  reference  equipment  using  delta-modulated  instruments. 
The  equipment  considered  consists  of  three  body-mounted,   delta-modulated, 
single-degree-of-freedom  integrating  gyros,   three  body-mounted,    delta- 
modulated,    single-axis  pendulum  accelerometers,    and  a  transformation 
computer.     The  delta-modulated  gyros  deliver  output  pulses  representing 
increments  of  rotation  and  the  delta-modulated  accelerometers  deliver 
output  pulses  representing  increments  of  velofciity      This  representation  le  ids 
naturally  to  the  choice  of  an  incremental  computer  10  perform  the  coordinate 
transformation  .    An  examination  of  performance  required  from  the  equip- 
ment concludes  the  investigation.     Following  is  a  brief  summary  of  the 
specific  points  covered,   and  the  specific  conclusions  reached. 

1.  Gimballess  inertial  reference  equipment  is  theoretically  feasible, 
and  the  required  computation  can  be  done  with  DDA's  currently  available. 
The  major  obstacle  to  the  realization  of  practical  gimballess  inertial 
reference  equipment  is  the  lack  of  inertial  instruments  specifically 
designed  for  this  application. 

2.  Delta  modulation  is  a  pulse  modul  ition  system  in  which  each  out- 
put pulse  represents  a  fixed  increment  in  the  input,    or  one  of  its  integrals 
or  derivatives.     The  analysis  of  a  deita-modulated  loop  proceeds  in  a 
straightforward  manner,    using  a  moding  plot  to  determine  the  modes  of 
the  system  and  Bergen's  technique:  for  finding  the  detailed  behavior  of  the 
limit  cycles.       In  a  loop  with  a  three -level  relay,   the  width  of  the  zero 
level  may  be  adjusted  to  eliminate  the  limit  cycle.     The  response  of 

the  loop  to  inputs  may    be    found    by    determining    bounds    on    the  con- 
secutive   numbers    of    pulses    possible    under    certain  conditions. 


The  output  sequences  may  be  verified  by  a  modified  form  of  Bergen's 
method.    For  small  inputs  a  modified  step-by-step  transient  approach 
is  best.    This  technique  relies  on  detailed  knowledge  of  the  system's 
modes. 

3.  Delta -modulated  integrating  gyros  and  accelerometers  are 
useful  inertial  instruments.    They  deliver  digital  outputs  in  synchronism 
with  a  clock,   and  can  be  scaled  as  necessary  to  meet  computational 
requirements.    Delta-modulated  instruments  are  fixed  level  torquing, 
which  is  much  more  accurate  than  analog  torquing;  and  the  measurement 
of  integrated  torque,   which  indicates  the  system  output,    is  merely  a 
pulse -counting  process. 

4.  The  transformation  of  the  A  V  vector  from  a  body-fixed 
coordinate  frame  to  an  inertial  coordinate  frame  is  best  accomplished 
through  a  direction  cosine  matrix  transformation.    The  direct  application 
of  the  matrix  differential  equation  is  most  suitable  for  digital  computation. 
Other  methods  considered  for  transforming  vectors,   while  providing 
valuable  insight  into  the  problem,   are  computationally  cumbersome. 

5.  Since  the  instrument  outputs  are  incremental  in  form,   the 
digital  differential  analyzer  is  the  most  attractive  choice  for  instrumenting 
the  solution  of  the  matrix  differential  equation.    Their  speed  and  compact- 
ness are  attractive  features  for  airborne  computers. 

6.  Integration  rules  with  order  higher  than  zero,   using  a  larger 
increment,    cannot  be  used  to  attain  a  slower  computer  speed  with 
increased  accuracy  because  of  an  irreducible  uncertainty  in  gyro 
information.    This  uncertainty  is  a  consequence  of  the  digital  form  of 
the  angular  information.    Thus,   rectangular  integration,   being  the 
simplest,   is  the  best  for  this  application. 

7.  The  transformation  computer  may  be  analyzed  using  frequency 
domain  techniques.    With  the  assumption  that  the  angular  rate  about 

one  axis  is  a  rational  multiple  of  that  about  another,   the  computer 
may  be  modeled  as  a  sampled-data  system,   and  a  simple  signal  flow 
graph  drawn  for  it.    This  signal  flow  graph  provides  great  insight  into 
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the  computer  operation,   showing  the  dominant  mode,   and  giving  an  esti- 
mate of  truncation  error.     Quantization  may  be  treated  statistically. 

7.  B      Suggestions  for  Further  Study  and  Development 

The  areas  of  further  study  and  development  which  can  contribute 
to  the  realization  of  practical  gimballess  inertial  reference  equipment 
are  listed  below. 

1.  Development  of  gyroscopes  designed  specifically  for  gimballess 
applications  should  be  pursued.     The  severe  angular  velocity  environment 
to  which  these  instruments  will  be  subjected  provides  enough  design  prob- 
lems to  warrant  a  separate  design.     A  torque  generator  providing  high 
torque  very  accurately  is  the  most  important  element  in  such  a  design. 

2.  Moving- coil  and  flapper-type  torque  and  signal  generators  should 
be  developed.  Their  advantages  of  small  size  and  greater  sensitivity  com- 
bine to  make  them  an  attractive  choice  for  use  in  inertial  instruments. 

3.  The  frequency  domain  analysis  of  the  transformation  computer 
should  be  extended  to  more  complicated  systems.     A  method  should  be 
developed  which  generalizes  the  results  for  rates  which  are  not  related  to 
each  other  as  rational  multiples  to  cases  where  the  rates  are  irrational 
multiples  and  to  cases  where  the  rates  have  varying  relationships. 

4.  The  reversibility  of  the  computation  is  an  interesting  point  which 
deserves  further  study.     The  logic  required  in  the  computer,   as  well  as 
the  actual  method  of  computation,    should  be  considered. 

5.  Stability  in  the  large  of  the  limit  cycles  of  delta-modulated  loops 
should  be  investigated. 
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APPENDIX  A 


DETAILED  ANALYSIS  OF  A  DELTA -MODULATED  LOOP 
WITH  A  TWO-LEVEL  RELAY 


The  system  to  be  analyzed  is  shown  in  Fig.   A-l.     This  is  the  loop 
considered  by  many  investigators  with  T     =7=1  second,   D  =  1.     The 
moding  analysis  has  been  carried  out  in  Chapter  II,   Fig.    2-4.     The 
Moding  Plot  indicates  the  existence  of  the  1-1,   2-2,   and  3-3  modes. 
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Fig.  A-l    Delta- Modulated  Loop  to  be  Analyzed 
A.l    Zero  Input  Behavior 

A  detailed  analysis  of  the  3-3  mode  follows: 
The  output  of  the  sampler  is 

C(z)  =    (1  +  z  +  z2j(l-z3) 
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(A-l) 


-sT 
Here  z  is  a  unit  delay,   z  =  e        s. 


The  output  of  the  linear  plant  is 
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S  S(7S+1) 

The  z -transform  of  A,   giving  the  samples  of  A,   is 


A(z)  =  -C(z)  G"(s)  + 


A 


(1+z+z2)  (1-z3)  (1-z) 


1-z 


1-z 


Lo      +     7A 


I 


tA 
o 

5l' 

s(ts+1) 

T  z 

s 

7 

l(l-z)' 


1-z 


°    [~PT  ~   1-o-z] 


l-az 


(A-2) 


(A-3) 


-T 

where  en  =  e    "^     . 

If  the  3-3  mode  is  possible,  A(z)  must  be  expressible  in  the  form 


A(z)  = 


A    +  A,z  +  A0z2  +  A.z3  +  A.z4  +  A.z5 
o         1  2  3  4 5 


1  -  z 


6 


(A-4) 


with  the  constraints 

A     >  0  A,  >  0  A„  >  0    A„<0A/0A,<0         (A-5) 

O  1  £  6  4  0 

so  that  the  postulated  relay  action  is  sustained.     If  Eq.   (A-3)  is  to  be  put 
into  the  form  of  Eq.   (A-4),   the  factor  (1-arz)  must  be  eliminated  from  the 
denominator  of  Eq.   (A-3).     For  this  factor  to  be  eliminated,  the  numerator, 
M     (z),   of  the  terms  having  (1-z   )  (l-az)  as  denominator  must  contain  that 
factor.     The  numerator  is 

M(z)  =  -(1-z3)2t-<1-z3)  tA 


o 

-  7  (l-z3)[l-z3  +  A    (1+z3)] 

o 


(A-u) 


1  -  A 


T  (1+AQ)  (l-z°)(l  - 


1  +  A 


o      3  v 
—  z     ) 


If  M(z)  is  to  contain  (1-arz)  as  a  factor  A     must  be  such  that 


^o        3 

■.■■  •'  =  a 


1+A 
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(A-7) 

1      3 
A     =  ,    =    .905 

o         ,\    3 
1+or 

This  is  one  condition  on  the  state  of  the  linear  plant  for  no  transient  in 

the  limit  cycle.     With  this  condition  satisfied,  the  output  of  the  linear 

plant  is 

n      3v2         /    T  z  „  \     ,  ,     .    2,    3,    4,    5 

A(z)  =  -(±^J-  -*  -  -J-  +  1+z+z  +z  +z  +z       (A  +7A    ) 

1-z6  \(l^z?  1-z/  1-z6  °        ° 

27        (l+<*z+a2z2)(l-z3) 
l+<*3  1-z6 

The  last  term  in  Eq.    (A-8)  corresponds  to  M(z)  after  cancellation. 
Collecting  powers  of  z  in  Eq.   (A-8)  gives 

(1-z6)    A(z)  =  Aq 

+  (A     - -IS"   +  1=SL  r-T+l)z 

°       l+cx3        1+*3  S 

+  (A     _J^+il£l7  _2T     +7)z2 
°        l+a3        l+o3  s 

+  (A    +— ^+     l~a„    7-3T-T)  z3 
°       1+*3         l+o3  s 

+  (A    +2£2-+^4-7-2T  -7)-4 

°     l+a3  1+a3  s 

+  (A  +^!i+  i4-T-T  -7>z5 

°        1+^        l+aJ  s  (A -9) 

The  numbers  above  T     =  7  =  1  sec)  lead  to 

(1-z6)  A(z)  =  Aq  +  (A     +  .204)z  +  (AQ-  .  353)z2  +  (Aq-  1. 195)z3 

+  (A     -  1.394)z4+  (A     -  .837)z5  (A -10) 
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By  applying  the  inequalities  Eq.    (A- 5)  to  Eq.   (A -10)  the  following 

bounds  on  A     are  found: 
o 

.353  <  A    <  .837  (A -11) 

o 

Thus,   with  initial  conditions  as  specified  in  Eq.    (A-7)  and  Eq.    (A-ll),  the 

3-3  mode  will  occur  without  transient. 

Similar  computations  lead  to  no  transient  initial  conditions  for  the 

1-1  and  2-2  modes: 

1-1  mode:    0  <  A    <  .076 
o 

^o=W='463  (A"12) 

2-2  mode:       0  <  A    <  .  475 

o 

i-   2 
A     =-L:£o-  =    .762  (A -13) 

°       l+a 

Further  computations  show  that  no  other  modes  are  possible. 
A.  2    Transient  Response 
A.  2.  a    Impulse  Response 

The  impulse  response  of  the  loop  is  found  by  a  modified  step -by -step 
computation.     After  an  input  impulse,   the  output  pulses  remain  one  sign 
for  a  number  of  pulses  determined  by  the  strength  of  the  impulse,   as 
the  increments  of  the  output  corresponding  to  the  integrated  impulse 
build  up.     As  the  output  impulses  pass  through  the  zero  order  hold, 
they  become  a  step.     Thus,   the  appropriate  calculation  to  find  the  state 
of  the  linear  plant  as  the  error  dies  out  is  a  combination  of  step  and 
impulse  response. 

For  a  numerical  example,   assume  that  the  loop  is  in  the  1-1  mode 
centered  at  the  origin.     Let  the  input  impulse  be  applied  halfway  between 
sampling  instants,    as  the  linear  plant  output  is  crossing  from  negative 
to  positive.     Let  the  imput  impulse  have  strength  S  =  5.5. 

With  t  =  0  at  the  clock  pulse  just  preceding  the  application  of  the 
impulse,   the  following  equations  describe  the  linear  plant  output 

A(t)  =  t  -  .  038  -  1.  463(l-et)    0  I  t  ^  .  5  (A -14a) 

A(t)  =  t  -  .  377  +  4.  612  0re"(t"'  5)  \  5  ^  t  ^  1  (A -14b) 


114 


m 
m 
ii 


CO 

r- 1 


a, 
o 
o 


CD 
OT 

G 
O 

a 

TO 

<D 

CM 
< 

•r-i 
fa 


CO 


<D 


C\J 


115 


A(t)  =  -(t-1)  +  2.403  +  3.8(l-e"(t_1))l^    t  (A -14c) 

Eq.   (A -14c)  changes  sign  between  t  =  7  and  t  =  8.     At  t  =  8,   the  state  of 
the  system  is 

A(8)  =  .  792  and  A(8)  =  -1  (A -15) 

These  initial  conditions  lead  to  the  3-3  mode,   as  may  be  seen  by 
comparison  with  Eq.    (A-7)  and  Eq.    (A-ll). 

Note  that  the  input  integral  is  only  determined  to  within  three  units. 
This  is  a  quantization  uncertainty  due  to  moding.     If  the  system  had 
settled  into  the  1-1  mode,  the  quantization  uncertainty  would  have  been 
only  one  unit. 
A. 2.b    Step  Response 

The  first  task  in  determining  the  step  response  of  the  loop  is  finding 

the  upper  and  lower  bounds  on  the  number  of  consecutive  pulses  of  one 

sign.     To  determine  the  upper  bound  on  the  number  of  pulses  of  one  sign 

the  maximum  velocity  and  displacement  of  the  linear  plant  output  are 

assumed,   and  the  number  of  pulses  required  to  drive  the  plant  across 

null  is  counted.     If  the  magnitude  of  the  input  step  is  R  =  rD,  the 

maximum  rate  that  the  system  can  attain  is  D(l+r).     If  a  sampling  instant 

occurs  just  before  the  system  crosses  null,  the  maximum  possible 

displacement  is  D(l+r)T   .     With  these  initial  conditions  on  the  linear 

s 

plant  and  the  linear  plant  being  driven  negatively  by  a  step  of  magnitude 
D(-l+r),  the  time  response  of  the  system  may  be  written  as 

nig 
A(nT  )  =  D(-l+r)nT    +  D(l+r)T    +  2Dr(l-e       7  )  (A -16) 

s  s  s 

If  n  is  to  be  the  upper  bound  on  consecutive  positive  pulses,   it  must  be 

the  least  integer  such  that 

A(nT  )<0  (A-17) 

s 

For  this  to  be  the  case,   the  relation  between  r,   n,    and  the  loop  parameters 

must  be 

T 

7  -n-§ 

n-1-2  ^(1-e       7)  (A_lg) 

r      n"+T 

Note  that  D  does  not  appear  in  the  above  expression.     For  the  numerical 
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example  above  this  relation  is 

>     n-3  +  2e"n 
r  n  +  1  (A -19) 

Table  A-l  is  the  result  of  substitution  of  various  n  in  Eq.   (A-19). 

While  this  eliminates  some  possible  sequences,   the  number 
remaining  to  be  tested  is  still  quite  large;  thus  it  is  desirable  to  place 
a  lower  bound  on  the  number  of  pulses  of  one  sign  which  can  follow  a 
given  number  of  pulses  of  the  opposite  sign. 

The  computation  to  determine  this  lower  bound  is  similar  to  that 
above.     The  initial  conditions  to  be  used  are  minima.     The  initial 
position  is  taken  as  zero  plus,   and  the  minimum  velocity  after  n, 
negative  pulses,   resulting  from  maximum  velocity  at  the  previous 
switching  time*  is 

With  the  above  initial  conditions,   the  expression  for  the  linear  plant 
output  is 

T  T 

s  s 

~~i~  ni  ~~T  n 

A(nT  )  =  D(-l+r)nT    +  2Dr(l-e  )  (1-e  )  (A-21) 

As  before,   if  n  is  to  be  the  required  lower  bound,   it  must  be  the  least 
n  for  which  Eq.   (A -17)  is  valid.     The  required  relationship  is 

T  T 

s  s 

2D  t£-  (1-e"  T~  nl)  (l-e"~  n)  (A -22) 

r  <  1 — 

n 

The  ranges  of  r  for  various  n  and  n,  are  shown  for  the  numerical 

example  above  in  Table  A- 2. 

The  response  of  the  loop  for  r  =  -jr  will  now  be  computed.     For  r  =  -*-    > 

Table  A-l  shows  that  while  the  feedback  and  input  act  together,   the 
maximum  number  of  pulses  (which  are  negative)  is  two,   and  when  they 
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TABLE  A-l 
UPPER  BOUND  ON  CONSECUTIVE  POSITIVE  PULSES 


n 

range  of  r 

1 

K-.634 

2 

-,634<r<-.233 

3 

-.233<r<    .025 

4 

r.02.5<r<   .207 

5 

.207<r<    .336 

6 

.336<r<    .429 

7 
i 

.429<r<    .5 

TABLE  A-2 

LOWER  BOUND  ON  CONSECUTIVE  POSITIVE  PULSES 
FOLLOWING  A  GIVEN  NUMBER  OF  NEGATIVE  PULSES 


range  of  r 

IVi 

1 

2 

3 

4 

1 

-    1     <r<.2 

-     1     <r«:095 

-     1     <r<-.2 

-     1     <r<-.  24 

2      1 

j    .2      <r<.456 

-.095<r<.251 

-.2      <r<.  .  18 

-.24  <r<.   15 

3     ! 

.  456<r<.6 

.261<r<.  456 

.  18  <r<    .  4 

.  15  <r<.38 

4    ! 

.6      <r<.S88 

.452<r<.  579 

.4     <r<    .533 

.  38  <r<.  518 

5 

1    .388<r<.749 

. 579<r<. 856 

. 533<r<   . 622 

.  5 18<r<..6  1 1      ! 

NOTE:    Bounds  on  negative  pulses  may  be  found  by  changing  the  sign  of  r 
and  entering  the  tables  again. 


are  in  opposition  the  maximum  number  permitted  is  seven.     For  the 
feedback  and  input  acting  together,   one  is  the  lower  bound  on  the 
number  of  negative  pulses  following  any  number  of  positive  pulses, 
according  to  Table  A-2.     However  at  least  three  positive  pulses  must 
follow  one  negative  pulse,  and  at  least  four  positive  pulses  must 
follow  two  negative  pulses. 

For  r  =  -p  the  following  is  a  partial  list  of  sequences  which  have 

the  required  average: 

Period  Sequences 

4T  3,1 

s 

8T  4,1,2,1;    5,1,1,1;    6,2. 

s 

12Tg  4,1,3,1,2,1;    4,1,2,1,3,1;    5,1,3,1,1,1; 

5,1,1,1,3,1;    6,2,3,1;    7,1,2,2;    7,2,2,1. 
According  to  the  bounds  found  above,  the  only  sequences  which  are 
possible  are  3,1  and  6,2.     The  sequence     3,1  will  be  verified  in  a 
manner  similar  to  that  used  in  verifying  modes  above. 
The  output  of  the  loop  is 


C(z)  = 


1+  -u   2      3 
1+z+z    -z 


(A-23) 


The  output  of  the  linear  plant  is 


A(s,z)  =  [R(s)  -  C(z)G2(s)]G{s) 


+ 


The  samples  of  the  linear  plant  are  given  by 


A(z)  =  D[r-C(z)  (1-z)] 


A    +  tA 
o  o 


T  z 

s 

L(i-z): 

tA 


+ 


+ 


tA 


S(7S+1) 


1-z 


l-az  - 


(A -24a) 


(A-24b) 


1-z  l-az 

If  the  postulated  sequence  is  possible,   then  A(z)  must  be  expressible 
as 
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A(z)  = 


A    +  A,z  +  A0z2  +  AQz3 
o 1 2 d 

1-z4 


(A-25) 


with 


A    >0    A,  >  0    A,  >  0  A,  <  0 
O  1  2.  6 


(A-26) 

For  this  to  be  the  case,   the  factor  (1-az)  must  be  eliminated  from  the 

denominator  of  Eq.   (A-24).     The  numerator  of  terms  with  denominator 

4 
(1-z   )  (1-az)  is 

M(z)  =  [.5(l-z4)  -  (l+z+z2-z3)(l-z)]7-  7Aq(1-z4) 
=  T  (l-z)[(.5-A   )  (l+z+z2+z3)  -  (l+z+z2-z3)l 


=  -  7(.5+A  )  (1-z)  (1+z+z2 ' 2-  z3) 

o 

.5+A 


(A-27) 


For  (1-az)  to  be  a  factor  of  M(z),   it  is  necessary  that: 


1.5-A 


2,     3 

—    =    a+ce    +  a 


.5+A 


or 


A     = 
o 


1  +  .  5  (l-a-a   -a   ) 

2 — 3 

1+a+a  +a 


With  A     as  in  Eq.   (A-28),  A(z)  may  be  written  as 


A(z)  = 


.5(l-z4)  -  (l+z+z2-z3)(l-z) 
1-z4 


Tsz 


(1-z)' 


1-z 


+  - 


(A-28) 


A  +tA 
o        o 


1-z 


tTr— j-a-z)  [i+(i+a)z  +  a+a+Q!2)z2] 


l+a+a  +a 
Collecting  powers  of  z  yields 

A  •  2rvT 

(l-z*)A(z)  =  A    +  (A  +tA    +.5T-.5T ^ — q)z 

000  s     l+a+a2+a3 


(A-29) 


+(A  +  tA  +.5T-T 
o        o  s 


2a*T    x   2 
)z 


1+a+a  +a 
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2C\+a+a>    )  i 

+(A  +7A     -1.5t-1.5T    + - !T— V)  z  (A-30) 

Substitution  of  the  numbers  of  the  example  into  Eq.   (A-30)  yields 

(l-z4)A(z)  =  A  +(A  +1.  312)z  =  (A  +3. 113)z2  +  A     -  .  25)z3        (A-31) 
Application  of  the  inequalities  Eq.   (A-26)  shows  that  for  3,   1  to  be  an 
output  sequence,  the  acceptable  range  of  initial  plant  values  is 

0  >  A     >  .25  (A-32) 

Similar  computations  show  that  6,  2  is  also  a  possible  output 
sequence,  but  that  none  of  the  others  is  possible. 

The  above  analysis,   as  noted  in  Chapter  II,   is  useful  only  for 
inputs  larger  than  about  10%  of  the  saturating  signal.     For  example,   an 
input  step  of  one -twentieth  of  the  saturating  signal  would  have  an  output 
sequence  of  some  multiple  of  forty  sampling  periods,   which  implies  a 
formidable  computation. 

For  small  inputs,  the  modified  step-by-step  transient  analysis 
described  in  Chapter  II  is  the  most  useful.     The  key  to  the  analysis  is 
the  fact  that  a  small  input  may  be  viewed  as  a  small  disturbance  of  the 
prevailing  limit  cycle. 

For  a  numerical  example,  assume  the  system  analyzed  above 
is  in  the  1-1  mode  with  zero  average  and  that  an  input  step  of  R  =  .  01D  is 
applied  at  the  sampling  instant  when  the  linear  plant  is  positive. 

The  plant  is  assumed  to  move  under  the  influence  of  the  input 
alone  using  the  initial  value  of  linear  plant  output  and  the  velocity  of 
the  plant  relative  to  the  limit  cycle  as  initial  conditions.     These  initial 

conditions  are 

. 

A     =  .038  A    =  0  (A-33) 

o  o 

and  the  equation  describing  the  motion  is 

t 
A(t)  =  .OlDt  +  .036  -  .01DT(l-e"T)  (A-34) 

The  mode  will  change  when  the  plant  output  sample  which  is  now 

negative  becomes  positive.     When  this  occurs*t  =  (2K=1)T  ,   and 
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A[   0C+1)T  ]>.  076,  that  is,   the  linear  plant  output  sample  gets  outside 
s 

its  permitted  range  for  sustaining  the  1-1  mode  [  see  Eq.   (A-12)]  . 

T 

A[  (2K+DT  ]    =  .01D(2K+1)T    +  . 036  -  .  01DT(l-e"  (iK+1)~ ) 

S  S  (A-35) 

For  K  =  2,  A[  2K+1)T  ]    =  .  07607,   and  the  mode  will  change.     The 

s 

velocity  of  the(linear  plant  at  this  time  is  the  drive  velocity,    .  01,   plus 
the  limit  cycle  velocity,    -.  463,   [  see  Eq.   (A-12)]  ,   and  the  position  of 
the  linear  plant  is  .07607  -  .076,   since  the  plant  output  which  was 
negative  is  now  being  considered.     One  sampling  period  later,   the 
linear  plant  output  is  (shifting  the  time  origin) 

A(T  )  =  -.  99  +  .  00007  +  (-.  453  +  .  99)  (1-e"1) 

A(T  )  =  -.651  (A-36) 

s 

The  plant  velocity  at  this  time  is 

A(T  )  =  -.792  (A-37) 

s 

By  comparing  this  state  of  the  linear  plant  with  the  mode  initial 

conditions,   it  is  seen  that  the  system  will  settle  into  the  3-3  mode. 

Again,   the  plant  is  assumed  to  move  under  the  influence  of  the  input 

alone,   with  initial  conditions  which  are  the  actual  plant  position  and 

the  velocity  of  the  plant  relative  to  the  limit  cycle.     Thus  the  initial 

conditions  to  be  used  are 

A     =-.651  A     =  -.792  -  (-.905)  =  +  .113  (A-38) 

o  o 

The  equation  of  motion  for  the  linear  plant  is  (again  shifting  the  time 

origin) 

t 
A(t)  =  .  OlDt  -  .  651  +  T(.  103)(l-e"T)  (A-39) 

The  mode  will  change  when  the  linear  plant  output  sample  which  is 

least  negative  becomes  positive.     This  occurs  at  t  =  (6K+2)T   .     At  this 

time,   the  linear  plant  output  will  leave  the  region  leading  to  the  3-3 

mode  so  that  A[  (6K+2)  Tg]  >  -  .  353  [  see  Eq.    (A-ll)]  . 

T 

A[6K+2)T  ]    =  .01D(6K+2)T     -  .  651  +  7(.  103)(l-e"(6K+2)~:r) 

s  s  (A-40) 
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For  K  =  3,  A[(6K+2)T  ]  *.'  *  348,   and  the  state  of  the  linear  plant  is 

A     =  .005,   and  A     =  +  .  739,   which  leads  to  the  2-2  mode.     This  mode 
o  o 

prevails  for  80  periods,  when  the  3-3  mode  appears.     This  mode 
persists  for  50  periods,   then  2-2  mode  prevails  for  50  periods.     The 
average  output  has  changed  one  unit  in  a  hundred  periods,  giving  the 
required  average  of  .  01.     A  plot  of  this  is  shown  in  Fig.  A- 3. 
A.  2.C    Effects  of  Imperfect  Integration 

With  imperfect  integration  in  the  linear  plant,  the  integrator 
is  replaced  by  a  first  order  lag.  In  that  case  the  linear  plant  in  the 
forward  path  becomes 

Gl(s)  =    (7lS  +  1)  (t2s  +  1)  (A-41) 

To  show  that  a  small  input  cannot  be  detected,    a  sequence  with 
zero  average  and  a  small  input  will  be  assumed.     The    inequalities 
resulting  from  the  output  sequence  will  show  the  size  of  the  input  that 
cannot  be  detected. 

Assume  an  input  R  =  rD  and  test  the  1-1  mode.     Then 


A(z)  =  r 


(1-z)2 
1-z     / 


1 


-T2/71      _      1-Tj/T. 


1-z  1-c,  z  l-CpZ 


A  +T0A  A  +t,A 

o      2    o  o     1    o 


l-a,z  +       l-anz  (A-42) 


T  T 

s  s 

where  a>,  =  e     T.  and  ct^  =  e      7~ 

If  the  postulated  output  sequence  is  to  be  sustained,  A(z)  must  be  of 

the  form 

A  +A,z 
A(z)=-°      l 


1-z2  (A-43) 


with 
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A     >  0 
o 


Ax  <  0 


(A-44) 

Again  the  (l-az)  factors  must  be  eliminated  from  the  denominator  of 
Eq.   (A -42).     The  constraints  on  the  initial  conditions  which  lead  to 
that  elimination  are 


A    +  T9A 
o         2    o 


=    r  + 


13 

1+or, 


l-a« 

A    +  T,A     =  r  +    .  ,    - 
°         l    °  1+a2  (A-45) 

With  these  conditions  met,   the  output  of  the  linear  plant  may  be  written 

as 


l-o 


1 


L-a, 


-7, 


'11+a,  2  1+a„ 

(1-z^)  A(z)  =  r  +  


Tl-72 


l-a. 


1 


l-a, 


1  l+« 


-T 


+  r  - 


1 


2  l+ac 


Tl  -  T2 


(a-46) 


For  the  postulated  mode  to  be  sustained,   the  limits  on  r  are 


l-a. 


1 


l-a. 


-T, 


1  1+a,       2    1+a 


1 


1 


=  IAo'r  =  0 


(A-47) 


Similar  calculations  show    that  the  magnitude  of  input  which  cannot 
be  detected  is  that  which  would,   in  the  steady  state,   drive  the  sample 
of  the  linear  plant  output  closest  to  zero  just  to  zero.     For  design, 
the  dead  zone  for  each  mode  must  be  found  and  the  largest  used  as  a 
limit.     As  was  seen  in  the  response  of  the  loop  for  a  small  input,   each 
mode  was  present,   and  the  largest  dead  zone  would  have  stopped  the 
dtrtput.     As  the  above  discussion  shows,   elastic  restraint  causes  an 
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error  in  the  output.     The  average  position  of  the  linear  plant  and  the 
elastic  restraint  constant  combine  to  give  an  average  force  on  the 
linear  plant  which  is  subtracted  from  the  input.     Thus,  the  output 
indicates  the  integral  of  the  difference  between  the  input  and  the 
average  elastic  force. 
A.  3    Frequency  Response 

The  describing  function  of  the  loop  shown  in  Fig.  A-l  will 
be  generated,  following  the  method  outlined  in  Chapter  II. 

The  input  to  the  relay  is  assumed  to  be  sinusoidal.     In  complex 
notation 

A(t)  =  Aeja,t  (A-48) 

Consequently 


The  error  signal  is 


C(t)    =  —  eJ    L  (A-49) 

7T 


E/t)  -         A         jut  A  iM-/G(ju)M  (A-50) 

Mt)        G(j«)      e  IG(jw)l    e 


The  feedback  signal  is 


sin  s  coT 

s 


B(t)    -±-  «       .»*— ■ g5-)  (A-5D 

s 
2 

The  input  is  the  sum  of  the  error  and  the  feedback  signal: 

R(T)  =  E(t)  +  B(t) 

coT 
A  S  ojT 


■>-,,.   u — e  "l — - —  t tz e  ^  e* 

G(jw)|  ir  coT 

(A-52) 

The  describing  function  is  the  ratio  of  output  to  input  amplitude 

4 
r    (w,  R)  = 


coT 

s  coT 

sin— ^ .  s 


A(R^L  e-JZG^)  +  A     -2  ej         2  g3 


G(JW)I  7T  COT 
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Im 


(R) 


I.  26 
(THRESHOLD) 


*-     Re 


R       1.26      1.4     1.6      1.8      2.0      2.5      3.0      4.0      6.0 
o       .995    .91      .79      .705.635   .51      .425    .32      .215 


R       -9°     -26°    -38°    -46°  -51°    -60°   -65°  -70°    -75° 


Fig.  A -4    Low  Frequency  Describing  Function  of  Loop 
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This  function  is  frequency-independent  at  low  frequencies  if 
there  is  an  integration  in  G(jw).     A  polar  plot  of  this  function  is  shown 
in  Fig.  A-4. 
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APPENDIX  B 


ANALYSIS  OF  DELTA-MODULATED  LOOP 


WITH  THREE-LEVEL  RELAY 


The  same  system  as  in  Appendix  A  will  be  analyzed,    with  a  three- 
level  relay  as  shown  in  Fig.   B-l. 


T 

s 

IF 

1 

A 

*c 

£ 

1 

I.  M. 

%(9>  E- 

C 

—  I 

i 
B 

s(ts  +  1) 

A 
c 

- 

D 

-sT 
1  -  e        s 

s 

Fig.   B-l    Delta -Modulated  Loop  to  be  Analyzed 
B.l    Choice  of  Zero-Level  Width 

As  explained  in  Chapter  II,   the  zero  level  of  the  relay  should  be 
chosen  to  be  as  narrow  as  possible  while  eliminating  the  limit  cycle. 
There  are  two  aspects  of  limit  cycle  elimination  to  be  considered:  the 
low-energy  region,   handled  oy  a  transient  analysis,    and  the  high-energy 


region,  handled  by  a  transient  analysis,  and  the  high -energy  region, 
handled  by  a  describing  function  analysis.     Each  analysis  gives  a  lower 
limit  on  the  zero  level,   and  the  larger  is  chosen  as  the  size  to  be  used. 

The  low -energy  analysis  assumes  that  the  linear  plant  output  is 
at  the  edge  of  the  zero  level  with  zero  derivative,   and  that  one  clock 
pulse  occurs  to  drive  the  plant  back  to  null.     The  zero  level  is  chosen 
wide  enough  so  that  in  an  infinite  time  the  other  edge  of  the  zero  level 
is  not  reached. 

While  the  clock  pulse  is  driving  the  plant  toward  null,  the  plant 

output  is 

t 
A(t)  =  -Dt  +  Ac  +  TD(l-e"T)  (B-l) 

At  t  =  T 

T 

S 

A(T  )    =    -D(l-e    T)  (B-2) 

s 

After  t  =  T 
s 

T 
s  t 

A(t)  =  -DT    +  A    +  TD(l-e     T)  -  TD(l-e'r)  (B-3) 

s  c 

if  A  (oo)  >  -  A  ,      A    must  be  chosen  such  that 

A    >-^-DT  (B-4) 

c       2         s 

This  will  be  the  result  of  any  low -energy  analysis,  no  matter  what  the 

dynamics,  as  long  as  the  plant  contains  one  integration. 

The  high-energy  analysis  makes  use  of  the  describing  function  for 

the  relay  which  is 


N(A)  (B_5) 

This  has  a  maximum,  and  hence  -^-has  a  minimum  at  A  =  n/2A     where 


N<^Ac>    ■    ^A  (B-6) 

c 

The  size  of  the  zero  level  is  chosen  so  that  there  is  no  intersection 

of  the  -  -^rpLocus  with  the  extended  G(jw)  locus.     The  magnitude  of  the 

G(j<*>)  locus  is  critical  for  the  3-3  mode,  the  possible  mode  with  the 
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largest  amplitude.     At  that  frequency 


CO 


GO  HT>! 


CO 


CO 


J7T   <l  +  3  IT  7> 


1+<J 


3T 


J_ 
T 


1     1/2 


(B-7) 


CO 


S  1 

Equating    I G(j  — jH  land   ,„  .yields 


A    >f  DT 
c       2  s 


1+<3    TT> 

s 


2-1     1/2 


(B-8) 


For  slow  systems  (  —=-  <  1 ),  this  width  for  the  zero  level  is  smaller 
than  the  low-energy  requirement,   and  for  faster  systems  it  is  somewhat 
larger.     However,   the  factor  multiplying  -r-DT    depends  on  the  highest 

Ci  S 

permissible  mode.     For  the  case  being  considered  here  (T     =  t)  Eq. 

s 

(B-8)  becomes 

A     >  4"  DT      (2.08)    =    1.04  (B-9) 

c       2,  s 

Since  this  is  larger  than  the  size  specified  by  Eq.    (B-4),   it  will  be 
taken  as  the  size  of  A    for  this  example. 
B.  2    Transient  Response 
B.2.a    Impulse  Response 

The  impulse  response  is  again  handled  in  a  step-by-step  transient 
manner,   recognizing  that  an  impulse  input  changes  the  derivate 
of  the  input  instantaneously  by  the  strength  of  the  input. 

For  a  numerical  example,  take  the  input  impulse  strength  as  5.  5 
and  assume  zero  initial  conditions.     At  the  end  of  one  sampling  period, 
the  linear  plant  state  is 

T 


A(T  )  =  5.5DT(l-e 
s 


T  )  =  3.48 
T 


A(Ts)    = 


5.  5De" 


T 


=    2.02 


(3-10) 


At  this  time  the  feedback  signal  appears;  the  linear  plant  output  is 
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T  T 

-_ s  -    g 

A(nT  )  =  -DnT    +  5.  5Dr(l-e   rj  +  Dt(5.  Se"^  *-  1) 
s  s 

T 

-n-   s 


(1-e        T)  (B-ll) 


The  feedback  signal  will  cease  when  A(nT  )  <  A     or  n  =  6  in  this 
°  s  c 

example  ..     The  state  of  the  linear  plant  is  then 

A(6T  )    =    .5 
s 

A(6T  )    =    -1  (B-12) 

The  plant  continues  under  the  influence  of  the  remaining  velocity,  but 

no  feedback  signal: 

-   t 
A(t)  =  .  5DTs-fD(l-e  T) 

and  #  (B-13) 

A(«o)  =  .5DT   -tD  =  -.5 

s 

Thus  no  negative  pulse  is  generated. 

A  plot  of  the  loop's  impulse  response  is  shown  in  Fig.    B-2. 
B.2.B    Step  Response 

The  step  response  of  the  loop  with  a  three -level  relay  is  computed 
in  much  the  same  way  as  that  of  the  loop  with  a  two -level  relay.     The 
bounds  to  be  determined  in  this  case  are  the  number  of  consecutive 
output  pulses  of  the  sign  of  the  input  permitted  and  the  number  of 
consecutive  pulses  of  sign  opposite  to  the  input. 

The  first  computation  proceeds  from  the  assumption  that  the 

linear  plant  output  has  travelled  a  full  sampling  period  under  the 

influence  of  the  input  alone  at  maximum  velocity.     This  leads  to  initial 

conditions  as  the  feedback  signal  appears 

A      =    A      +    DrT 
°  C  S  (B-14) 

A      =  Dr 
o 

The  linear  plant  output  is  then 

T 

_    s 

A(nT  )  =  -D(l-r)nT    +A    +  DrT    +TD(l-e  "H  (B-15) 

s  s         c  s 

The  maximum  number  of  output  pulses  of  the  same  sign  as  the  input 
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is  that  n  for  which  A(nT   )<+A  •     If  this  is  the  case,   then 

S  c 

T 

T     /i       -n  s  . 

n-  y    <1-e         ~t~) 

r  >  nTT (B-16) 

Substitution  of  various  n  into  Eq.   (B-16)  yields  Table  B-l,  which 
shows  limits  on  r  for  a  given  maximum  number  of  output  pulses  of 
the  same  sign  as  the  input. 

Table  B-l 

Limits  on  Size  of  Input,  Given 

Number  of  Consecutive  Output  Pulses 

of  Same  Sign  as  Input 


n 

Range  of  r 

1 

1  rl    >     0 

2 

Irl    >  .181 

3 

1  rl    >  .378 

4 

1  rl    >  .501 

5 

1  rl    >  .604 

6 

1  rl    >  .668 

The  possibility  of  pulses  of  opposite  sign  from  the  input  must  be 

considered.     If  the  feedback  and  input  act  in  opposition  for  a  full  sampling 

period  before  switching  occurs,   the  state  of  the  linear  plant  is 

A      =    A      -    D(l-r)T 
o  c  s 

(B-17) 


AQ  =-D(l-r) 

Even  for  the  smallest  input,  A    >  -A     so  there  will  be  no  feedback 

r  o  c 

signal  initially.     Thus  the  linear  plant  output  is 

-n  — § 
A(nT  )  =  DrnT    +A     -  D(l-r)T     -   TD(l-e      f    )  (B-18) 

s  s         c  s 

If  ever  A(nT  )  <  -  A  ,  there  will  be  a  negative  pulse;  thus  the  expression 

T 

T  s      2A 

1  +   y  ^-e_n  "  T)  "  W 
0>r> § ?7T1 1  (B_19) 
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is  of  interest.     Clearly  r  has  been  assumed  positive.    For  the  present 
example,  there  can  be  no  negative  pulses  for  a  positive  input,   and  vice 
versa. 

The  output  for  an  input  step  R=o"  ^  w^  now  ^e  computed. 
According  to  Table  B-l,  there  may  be  one,  two,  or  three  positive  pulses, 
but  not  four.     However,  the  only  possible  output  sequences  are  1, 1, 0,  0, 
and  2,  2,  0,  0.     These  may  be  verified  in  the  same  manner  as  the  sequences 
tested  in  Appendix  A. 

For  an  input  step  of  R  =  .  01D,  the  only  possible  output  sequence  is 

1,  99,0,0.     That  no  more  than  one  consecutive  positive  output  pulse 

is  possible  is  clear  from  Table  B-l;  thus,  the  above  sequence  alone  is 

permitted.    A  short  transient  calculation  shows  that  the  first  output 

pulse  will  occur  at  t  =  104T  ,  because  of  the  width  of  the  zero  level. 

s 

B.  2.  C    Effects  of  Imperfect  Integration 

As  in  the  previous  case,  the  effect  of  imperfect  integration  is  to 
mask  a  small  input.     The  limit  on  detectable  inputs  can  be  determined'* 
by  comparing  the  steady-state  response  of  the  linear  plant  with  the  zero- 
level  width  of  the  relay.     Inputs  which  cause  steady -state  displacement 
of  the  linear  plant  less  than  A    will  not  be  detected.     Further,   as  before, 
an  input  corresponding  to  the  average  value  of  the  linear  plant  output 
will  be  subtracted  from  the  indicated  output. 
B. 3    Frequency  Response 

As  in  Appendix  A,  a  frequency  dependent  describing  function  may 
be  calculated  for  the  loop,  beginning  with  an  assumption  about  the  relay 
input.     Assume  this  is 

A(t)  =  AeJwt  (B-20) 


The  output  of  the  loop  is 


the  error  signal  is 


(B-21) 


*»;-<&''?*  '-fair  J*****        B"22) 


and  the  feedback  signal  is 
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2 
The  input  is  the  sum  of  the  feedback  and  error  signals,  so  that  the 
describing  function  of  the  loop  is 


•t 


^ 


[A     \ 
c 

lA 


2 


f(A,«)=  r J1       «Ts    J  wTs 

A_e-J^ti«>  +  i\A-  -£*-'       sin    2    e       2 


lG(jw)lc  7T  Vx       A  ~~ £T 

s 

2 

(B-24) 

This  may  be  plotted  as  a  function  of  R  for  each  frequency.     Unlike  the 
describing  function  for  the  loop  with  a  two -level  relay,  this  is  not  freqency 
independent  for  low  frequencies. 

Due  to  the  zero  level  of  the  relay,  there  is  a  magnitude  which  R 
must  attain,   depending  on  frequency,  before  any  output  appears.     It 
is  required  that 

A    =    RlG(jw)l>Ao  (B-25) 

or 

A 
c 


c 
A 


R>lG(jco)|  (B-26) 


for  an  output  to  occur. 
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APPENDIX  C 


ANALYSIS  OF  TRANSFORMATION  COMPUTER 


C.l    Roots  of  Characteristic  Equation 

The  roots  of  the  computer'  s  characteristic  equation  for  two  dif- 
ferent inputs  will  be  determined  using  the  root-locus  method. 

C.  1.  a    Equal  Rates  on  All  Axes 

In  this  case  uk   =  w2  =Wn,   so  z,  =  z2  =  z~  =  z.     Eq.   (5-27)  becomes 

1  +  3A02 - -  A63 - =  0  (C-D 

(1  -  z)2  (1  -  z)2 

As  is  shown  in  Fig.   C-l  the  roots  of  this  equation  remain  on  the 
unit  circle.     This  indicates  a  non-decaying  oscillatory  behavior,   as 
expected  of  direction  cosines. 

C.  l.b    Different  Rates  on  Different  Axes 

Suppose  that  the  angular  velocities  are  u,  =  2u2,  Wo  =  0.     The 
difference  equations  are  then 


Cn(2n)  =  C-j^n-l)  -  A9Ci3<2n-l) 


Cil(2n+1)  =  Cn(2n) 


C.2(n)  =  Ci2(n-1)  +  A9Ci3(n-l) 

C.3(2n)  =  Ci3(2n-1)  +  A9Cil(2n)  -  A0Ci2(2n) 

C.3(2n+1)  =  Ci3(2n)  -  A0Ci2(2n+l)  (C-2) 


(3A02-Ag3)Z 
(l-Z)2 


A0'  =  ^A0 


/"  A0    3A02   A03    A04 
Vl+1  4  2  4~ 


,     3A0*    A33 
2  2 


Fig.   C-l    Root  Locus  For  Equal  Rates  on  all  Axes 


140 


To  transform  these  equations  into  the  frequency  domain,  it  is 
oo 

V  2n 

necessary  to  compute  y    C-3(2n)z     ,  which  is  the  transformed  sequence 

n=0 
of  every  other  value  of  C.,.     This  is  the  sequence  which  is  fed  into  the 

C.,  node  (register). 

2Ci3(2n)z2n  =  2C.3(n)zn  -  2Ci3(2n+l)z2n+1 

=  2Ci3(n)zn  -  2C.3(2n)z2n+1  +  A0C.2(2n+l)z2n+1 

=  ZCi3(n)zn  -  zZCi3(2n)z2n  +  A02Ci2(2n)z2n+1  +  A022Ci3(2n)z2n+1 

=  ZC3(n)zn  -  z(l  -  A02)2Ci3(2n)z2n  +  A0z2Ci2(2n)z2n       (0-3) 

ZCi2(2n)z2n  =  ZCi2(n)zn  -  2Ci2(2n+l)z2n+1 

=  ZCi2(n)zn  -  2C.2(2n)z2n+1-  A0'.Z:Ci3(2n)z2rl+1 

2C2(2n)z2n  =  J—    ICi2(n)zn  -  A0zLCi3(2n)z2n  .  J     (C-4) 

1  +  z  \  I 


2n  ,  A0z 


ZC,(2n)z2n  =  IC,(n)zn  -  z  fl  +  AQ^  (1  +  -£—   )1  IC.9(2n)z^n+— IC.9(n)zn 
13  i3  L  1  +  z    J        l3  1  +  z        " 


SCi3(n)zn  +  **2-  2Ci2(n)zn 
2Ci3(2n)z2n  = L+5 (C-5) 


(l  +  z[l+A02(l  +  -^—  )]) 


1  +  z 

In  z- transform  theory,    the  z- transform  of  a  function  is  defined  as 
co  oo 

Cr(z)=;J    Ci.(n)zn  (Ref.    61,   Ch.   1).     Let  C! .(z)  =^    Ci;j(2n)z2n.     Then 

n=0  n=0 

Eq.   (C-5)  may  be  written 
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Ci3<z>-—  Ci2<z> 
C!3(z)  =  (C-6) 


(i  +  zHi-fAe2^-2-!) 

1         '  L       l  +  zJ' 

From  the  difference  equations  are  obtained  the  .relations 


C'    (z)  =  J_  C.,(z)  (C-7) 

11  1  +  z       ll 


C.l(z)  =  -^C'    (z)  (C-8) 

11  1-z      16 


Ci2(z)  =A^_C.3(z)  (C-9) 

1-z 


C.3(z)  =-^-Ci2(z)  -    *L  CUz)  (C-10) 

1-z  1-z 

Combining  Eqs.   (C-6)  through  (C-10)  gives  the  characteristic 
equation  for  the  computer: 

AQ2    [.       A02z2\ 
AG    z    1 


1  +  A9   Z       + l-^ =  0      (C-ll) 

(1  -  z)2       (1  -  z)2(l  +  z)[l  +  z  ll  +  A02fl  +-Z- 

L  l  L         1  +  : 

Fig.   C-2  shows  the  root-locus  diagrams  used  to  find  the 
singularities  of  Eq.   (C-ll).     There  are  four  important  zeros  of  this 
equation,  marked  with  squares  in  Fig.   C-2(b).     The  pair  of  roots  near 

unity  represent  the  dominant  mode  of  the  system.     The  pair  near  -1 

A0 
are  oscillatory  with  magnitude    1  — g-    •     This  gives  rise  to  a  time 

sequence  (1  --5— )    .     Since  the  magnitude  of  the  mode  is  O(A0)  the  pair 

2 

of  roots  gives  rise  to  an  error  growing  as  nA0    .     This  leads  to  an 

increment  of  error  in  about  one  revolution. 
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NOTE:    DISTANCES  OF 

SINGULARITIES  FROM 
+  I  AND -I 
ARE    O(A02). 


1  + 


(b) 
A02Z2 

(l-Z)2    V,+    (l  +  Z)(l+Z[l+A02[l+£z])/ 


ReZ 


Fig.   C-2    Exact  Root  Locus  for  co  =2co„,co   =0 
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C.  1.  c    Approximate  Analysis  of  Multirate  System 


For  u^  =  2u2,   u>3  =  0,   Eq.    (5-22)  shows  k  =  1,  I  =  2,   and  Eq.   (5-23) 


shows  z,     =  zo  : 


Substitution  of  these  relations  into  Eq.   (5-25) 
yields  Eq.    (5-21).     The  characteristic  equation  corresponding  to  this 
set  is 


2  2 

1  +A9   z      +  A9    z 


o  2 

a-zr    (i-z  ) 


=  0 


(C-12) 


Comparison  of  this  equation  with  Eq.    (C-ll)  shows  that  it  is  approxi- 
mately correct.     The  root  locus  diagram  for  this  equation,   shown  in 
Fig.   C-3,  is  practically  the  same  as  for  Eq.  (C-ll).     The  extra 
dipoles  appearing  in  Fig.   C-2  have  negligible  effect. 

C.  2    Truncation  Error 

Following  the  method  outlined  in  Chapter  V,  Sec.    5.  D.  3,   an 
analysis  of  the  truncation  error  made  by  the  computer  will  .be  made, 
first  for  equal  rates  on  all  axes,   then  for  u.  =  2to2,  Uo  =  0.     In  each 
case,   the  initial  condition  C(0)  =  I,  will  be  us~ed. 

C.  2.  a    Equal  Rates  on  All  Axes 

The  exact  solution  to  Eq.   (5-1)  in  this  case  is 

2 


C(s)  = 


1/s 


3W2 


+  1 


(C-13) 


The  time  solution  is  then 
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n  ImZ 


IZI  =  I 


(l-Z)2    \      (l+Z)2 


*-Re(Z) 


(b) 


Fig.   C-3    Approximate  Root  Locus  for  w  =  9oj      fii  =0 
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1  ♦  2  cos  VI  ut 
1  -  cos  VI  <■<  +  VI  tin  ut 
1  -  cos  VI  wt  -  VI  sin  wt 


1  -  co*  VI  wt  -  VI  sin  VI  wt 

1  +  2  cos  VI  wt 
1  -  cos  VI  wt  +  VI  sin  VI  wt 


+  VI»in  VI  wt    \ 


I  -  cos  VI  wt 
1  -  cos  ^wt  -  V3  sin  Vlwt 
1  +  2  cos  VI  wt 


(C-14) 


According  to  the  flow  graph  shown  in  Fig.    5-4,   the  computer 
solution,   without  quantization  is 


C(z) 


A92z 
(l-z)2 


TPz) 


,   ,   (3A92  +  £i63)z 
1  + n 

a-zr 


A0z       &0   z 


l-z        (1-zr 


Aflz       A02z2 
l-z        (l-z)2 


A0        A0'z 

+ 2 

l-z      (1-zr 


1  + 


A02z 
(l-z)2 


A0z       A0   z 


l-z        (l-z)' 


(C-15) 


The  normal  modes  of  this  determinant  have  unit  magnitude  and 
argument 


1  + 


AG        3A92         A93  -    A94 


A9'    -  \/'3A9 


(C-16) 


1 


3A92       A93 


2  4 

Thus  the  sequence  corresponding  to  Eq.    (C-15)  is  approximately 


1  +  2  cos  n  AH1 


1  .  M. 
CM  -  _  I    1  -  cos  n  Afl'  +  V3  sin  n  A9' 


1  -  cos  n  Afl'  -  -J5  9ln  n  Afl'  1 


n4«'  ♦<ffilnn  Afl'  \ 


1  +  2  cos  n  A8' 


,1  -  cos  n  A81  -  \fl"sin  n  A8'  1  -  cos  n  AS'  +  VTsin  n  A81 


1  -  cos  n  Afl'  -  VT  sin  n  A81 


1  +  2  cos  n  A9' 


(C-17) 
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The  truncation  error  arises  in  this  case  from  the  fact  that  the 

argument  of  the  roots  near  +1  is  not  quite  correct.     This  leads  to  an 

2 
error  which  grows  as  nA0   . 

C.  2.b    Multirate  Systems 

With  u,  =  2u>2,   Uo  =  0  and  C(0)  =  I,   the  exact  solution  to 
Eq.   (5-1)  is 


C(s)  =    1/s 


.4 


2 

s 

+  ±- 

4 

5 

2 

u 

4 

.8 

s 
u 

u 


(C-18) 


The  corresponding  time  solution  is 


C(t) 


1     iA  *R         *% 

—  (4  +  cos cut) 


2  n               ^    *\ 
—  (1  -  cos ut) 

5  2 

1      .     sR 


—  sin 
ln/5  2 


cot 


2  /I  ^   « 

—  (1  -  cos  —  ut) 
5  2 

l  /i     a        nHT    ** 

-(1  +  4  cos ut) 

5  2 


2      .    nTF    . 

—  sin ut 

■s/5  2 


1  .    sTS    \ 

-  —  sin cjt1 

^5  2 

2  .     4^    . 
—  sm wt 

\/5  2 

n/51      , 
cos cat 


(C-19) 


Under  the  same  conditions  the  model  of  the  computer  is  as 
shown  in  Fig.   C-4. 
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I-Z 


Fig.   C-4    Computer  Model  for  Wj  =  2w2,a>3F0;  C(0)-I 
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I  *• 


M' 


C(z) 


1_ 
-z 


1  + 


Afl2z .  Afl2z 


(1-z)2  (1-z2)2 


(1-z)' 


A92z 


(l-z)(l-z2) 


A92z 


(l-zKl-z") 


1  + 


A92z 
(1-z   ) 


A9 
(1-z) 


(C-20) 


A9z 
1-z 


Aflz 
1-z 


There  are  essentially  two  parts  to  the  time  solution  of 
Eq.   (C-20).     The  first,  desired  part  of  the  solution  arises  from 
the  roots  near  +1.     The  roots  near  -1  give  rise  to  extraneous,   un- 

and 


stable  terms.     The  roots  near  +1  have  magnitude  1  +  —  AG' 
argument  A9'  =  -^-  A9. 
(C-20)  is  approximately 


argument  A9'  =  —^-  A9.     The  time  sequence  represented  by  Eq. 


C(n)=(l-nMi 


-(4  +  cos  nA9'  ) 
5 

-  (1  -  cos  nA6') 


sinnAO'  ) 


n/5 
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—  (1  -  cos  nA9'  )  -— - -  sin  nA9 
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5  \l~5 
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sinnA9' 


cos  nA9' 


+nA9' 
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2 
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1 
2 

1 

1 

2 

4 

A9' 


A9' 

8 
A9J 

8 

0 


(C-21) 


Again,   an  estimate  of  the  truncation  error  may  be  had  by  com- 
paring Eq.   (C-21)  and  Eq.   (C-14).     The  error  arises  from  the  presence 
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of  the  poles  near  -1,  and  the  fact  that  the  roots  near  +1  are  not  on  the 

2 
unit  circle.     The  truncation  error  grows  as  nA9     in  each  direction 

cosine.     It  should  be  noted  that  the  approximations  made  in  the  analysis 

preclude  :taking  Eq.   (C-21)  as  the  exact  computer  solution. 

C.  3    Quantization  Error 


The  power  spectral  density  of  quantization  may  be  calculated 
following  the  procedure  indicated  in  Chapter  V,  Sec.    5.  D.  4.     However, 
perhaps  a  better  indication  of  the  error  due  to  quantization  is  the  mean- 
squared  signal  appearing  at  any  node  due  to  quantization  noise.     This 
will  be  done  for  the  case  of  w,,   =  2u)2,  Uo  -  0>  using  both  hysteretic  and 
non-hysteretic  quantizers. 

C.  3.  a    Non-Hysteretic  Quantization 

The  model  to  be  analyzed  is  shown  in  Fig.   C-5.     The  expression 
for  power  spectral  density  at  each  node'. is: 
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A9z 

i        2 

1-z 

l      A02z             A92z 
(1-z)2       (1  -  z2)2 
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2  2 
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1  +  Aq2z       +    Aq2z 


(1  -  z)2       (1  -  z2)2 


AG' 


12 


A9^ 
12 


(C-23) 
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Fig.   C-5    Computer  Model  with  Non-Hysteretic  Quantization 


Fig.   C-6    Computer  Model  with  Hysteretic  Quantization 
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Since  the  magnitudes  of  iHe  roots  are  very  close  to  unity,   the 
mean  square  value  will  be  very  large.     Thus,   only  n  terms  in  the  sum 
of  squares  will  be  considered.     Note  that  the  contributions  from  the 
roots  near  +1  are  oscillatory,   so  that  the  mean  square  may  be  taken 
as  half  the  amplitude  squared.     Let  qr*.  be  the  mean  squared  value  of 


..th 


ij 


quantization  noise  at  the  ij      node.     Carrying  through  the  computations 
indicated  above  gives  for  mean  squared  value  after  n  steps 

q^(n)  =  .  096  nA04 


q^n)  =  1.  2  nA04 
q^3(n)  =  .  0188  nA04 


(C-24) 


The  rms  quantization  noise  at  any  node,   then,  is  on  the  order  of 
^A02. 


C.3.b    Hysteretic  Quantization 

The  model  of  the  computer  with  hysteretic  quantization  is  shown 
in  Fig.  C-6.     Note  the  power  spectral  density  of  the  noise  is  four 
times  as  great  as  the  non-hysteretic  noise,  due  to  the  larger  quanti- 
zation box. 

The  self  loops  may  be  absorbed  into  the  other  paths,   giving  a 
factor  1/(1 --g  z)  multiplying  each  path  gain.     The  expressions  cor- 
responding to  Eq.   (C-23)  are: 
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(C-25) 


Performing  a  sum-of-squares.    computation  similar  to  the  above 


gives: 


q2L(n)  =  .  384  nA  94  +  1.  07  A04 


q22(n)  =  4.  8  nA04   +   5.  23  A04 


q23(n)  =  .07  52  nA04  +  .  158  A04 


(C-26) 


The  rms  error  is  a  factor  of  two  higher  than  before,   and  a 
constant  has  been  added. 
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21.  Charles  S.   Draper,   Walter  Wrigley,   and  Lester  R.   Grohe, 

The  Floating  Integrating  Gyro  and  its  Application  to  Geometrical 

Stabilization  Problems  on  Moving  Bases,  A.   Sherman  F.   Fairchild 

Fund  Paper  (FF-13)(Institute  of  the  Aeronautical  Sciences,   New  York), 

1955. 

General  discussion  of  gyros  and  their  application 
to  platform  stabilization.     Gyro  performance  is 
considered  in  some  detail. 
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22.  William  E.  Fellows,  "Vibration  Effects  on  Gyroscopes  and 

Accelerometers,  "  Minneapolis  Honeywell  Regulator  Co.  Aero 

Document  R-Ed    29004-1A,  November,  1958. 

Detailed  discussion  of  many  errors  caused 
by  vibrations. 

23.  Richard  H.  Frazier,  Analysis  of  Four  Pole  Torquer  with  Full 
Normal  to  Gap  Faces,   MIT/IL  Rpt.   E-700,   March,  1958. 

24.  Richard  H.  Frazier,  Four  Pole  Alterating  Current  Signal  Generator 

with  Flapper  Type  Rotor,   MIT/IL  Rpt.   E-703,  April,  1958. 

Each   reportis  a  brief  discussion,   with  numerical 
examples,   of  theory  of  flapper-rotor  instruments 
in  different  applications. 

25.  Albert  P.  Freeman  and  Peter  J.    Palmer,  Abridged  Nominal 
Values  and  Specification  for  2FBG-2  Series  Floated  Gyro  Unit,   MIT/IL 
Doc.   G7-100-U-2FBG-2C,   July,  1959. 

26.  Philip  J.   Gilinson  Jr. ,   William  Denhard,   Richard  H.   Frazier, 
A  Magnetic  Support  for  Floated  Inertial  Instruments,  A  Sherman, 

F.  Fairchild  Fund  Paper  (FF-27)(Institute  of  the  Aeronautical  Sciences* 

New  York), May,   1960. 

Discussion  of  theory  of  magnetic  suspension 
system  used  in  many  inertial  instruments. 
Numerical  examples  provided. 

27.  William  B.  Haff  and  Melvin  Meltzer,   Effects  of  Angular  Vibration 

on  the  Performance  of  a  Single -Degree-of -Freedom  Integrating  Gyroscope, j 

M.S.    Thesis,   M.  I.  T.    (MIT/IL  Rpt.    T-257),May,   1960. 

Theoretical  analysis  of  vibration  effects  on 
gyros.     Coning  and  sculling  discussed.     Series 
solution  for  equation  of  gyro  motion  proposed. 

28.  Robert  K.   Mueller,   Microsyn  Electromagnetic  Components 


(Instrumentation  Laboratory,   M.  I.  T.  ),  December,  1952. 

Fundamental  discussion  of  microsyns  in 
their  various  forms  by  the  inventor. 

29.     Herbert  Weinstock,  A  Study  of  the  Response  of  the  Single -Degree- 
of -Freedom  Integrating  Gyroscope  to  Angular  Vibrations,   MIT/IL 
Rpt.    E-885,  January,   1960. 

Consideration  of  simultaneous  angular 
vibrations  about  several  instrument 
axes  at  once. 
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C.     Coordinate  Transformations 

30.  Herbert  Goldstein,   Classical  Mechanics  (Cambridge,   Mass.,  1951). 

Chapter  IV  has  a  complete  discussion  of  matrix 
transformation,   Euler  angles,   and  Cayley-Klein 
parameters.     Good  annotated  bibliography. 

31.  Sir  William  R.   Hamilton,  Elements  of  Quaternions  (Dublin,  1856). 

This  treatise  introduces  quaternions  as  an 
extension  of  complex  numbers. 

32.  A.  S.  Hardy,   Elements  of  Quaternions  (Boston,   1881). 

Basic  text  on  quaternions  with  geometrical 
applications. 

33.  H.  W.  L.   Hime,    The  Outlines  of  Quaternions  (London,   1894). 

Quaternions  introduced  as  the  quotient  of 
two  vectors.     Various  forms  given. 

34.  Sir  Harold  Jeffreys  and  Bertha  Swirles  (Lady  Jeffreys),   Methods 

of  Mathematical  Physics.   (Cambridge,   England,   1956)  . 

Chapters  3  and  4  cover  matrices  and  tensors  very 
well.     Each  chapter  is  headed  by  an  apt  quotation. 

35.  Myron  Kayt on,   Coordinate  Frames  in  Inertial  Navigation,   Ph.D. 

Thesis,   M.  I.  T.    (MIT/IL  Rpt.    T-260)  August,   1960. 

General  discussion  of  coordinate  transformation 
theory  in  Chapter  4  using  tensors. 

36.  P.  Kelland  and  P.  G.    Tait,   Introduction  to  Quaternions  (London,  1873). 

Written  "for  those  who  desire  to  become 
mathematicians.  " 

37.  J.  Halcombe  Laning,   Jr.,  The  Vector  Analysis  of  Finite  Rotations 

and  Angles,  MIT/IL  Rpt.   6398-S-3,   September,   1949. 

An  interesting  development  of  algebra  required 
for  the  use  of  vectors  to  represent  angles  and 
rotations. 

38.  Philip  M.   Morse  and  Herman  Feshbach,   Methods  of  Theoretical 
Physics,   2  vols.    (New  York,  1953). 

Chapter  1  develops  some  transformation  properties 
and  contains  a  brief  mention  of  quaternions. 

39.  William  R.  Weems,  An  Introduction  to  the  Study  of  Gyroscopic 
Instruments.    (Department  of  Aeronautical  Engineering,   M.I.  T. ), 
January,  19  48. 
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Section  8  contains  a  n,ice  discussion  of  Euler  angles. 

40.  E.  T.    Whittaker,  A  Treatise  on  the  Analytical  Dynamics  of 

Particles  and  Rigid  Bodies,   4      ed.   (Cambridge,   England,   1959). 

Chapter  1  contains  discussion  of  rotational 
transformations  using  quaternions,  Cayley- 
Klein  parameters,  and  Euler  angles. 

D.     Numerical  Analysis  and  Computers 

41.  Anonymous,  "Modern  Computation  Methods,  V  Notes  on  Applied 

Science,  National  Physical  Laboratory  (of  England)(New  York,  1957). 

Complete  discussion  of  numerical  analysis  as 
applied  to  computers.     Excellent  extensive 
bibliography. 

42.  Vannevar  Bush,   "The  Differential  Analyzer,  "  J.  Franklin  Inst. , 

CCXII  (1931),   447-488. 

Detailed  physical  description  of  one  of  the 
first  analog  computers. 

43.  Vincent  Rocco  DeMarco,   The  Effect  of  Increment  Representation 

on  the  Accuracy  of  the  Digital  Differential  Analyzer,   M.  S.    Thesis, 

M.I.  T.    (MIT/IL  Rpt.    T-158),   May,  1958. 

Complete  discussion  of  errors  in  DDA  with 
various  increment  representations.     Analysis 
in  time  domain. 

rd 

44.  George  F.  Forbes,    The  Digital  Differential  Analyzer,   3 

ed.    (Pacoima,  Calif.,   1956). 

Discussion  and  analysis  of  DDA. 

45.  Laurence  J.   Gitten,   Incremental  Computer  Solution  of  a  Coordinate 

Conversion  Problem,   M.S.    Thesis,   M.  I.  T.  ,   1959. 

Study  of  a  coordinate  transformation  using 
a  general  incremental  computer. 

nd 

46.  David  R.   Hartree,   Numerical  Analysis,  2     ed. (Oxford,  1961). 

A  standard  reference,   includes  a  section  on 
organization  for  automatic  computation. 

47.  Francis  B.   Hildebrand,   Introduction  to  Numerical  Analysis 
(New  York,  1956). 

Standard  text  on  numerical  analysis. 
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E.   Miscellaneous 

56.  Lewis  Carrol,    Through  the  Looking  Glass,   and  What  Alice  Found 

There,  illus.   John  Tenniel  (New  York,  1941). 

The  quote  is  in  Chapter  I.     The  explanation, 
by  an  early  astronaut  who  had  trouble  on  the 
launching  pad,  is  in  ChapterVI. 

57.  Harold  Chestnut  and  Robert  W.   Mayer,   Servomechanisms  and 

Regulating  System  Design,   2  vols.    Vol  I,   2nd  ed.   (1959),   Vol.  II  (1955). 

Comprehensive  text  on  control  system  theory 
and  practice. 

58.  Ralph  J.  Kochenburger,  "A  Frequency  Response  Method  for 

Analyzing  and  Synthesizing  Contactor  Servomechanisms,"  Trans  AIEE 

Part  I,    LX1  (1950),   270-284. 

The  describing  function  presented  and  applied 
by  its  originator. 

59.  Albrecht  L.  Kosmala,  Feasibility  Study  of  a  Gimballess  Inertial 


SpEce  Reference,   MIT/IL  Rpt.   R-274,  April,   1960. 

Early  attempts  at  analysis  of  gimballess 
inertial  reference  equipment. 

60.  Samuel  J.   Mason  and  Henry  J.    Zimmerman,   Electronic  Circuits, 

Signals,   and  Systems  (New  York,   1960). 

Chapter  IV  contains  detailed  discussion 
of  flow  graph  analysis. 

61 .  Robert  W.   Sittler  et  al.  ,   Pulsed  Data  Systems  (Course  Notes  for 

6.54,   M.I.  T.  ,   Spring,  1960). 

Course  notes  for  sampled  data  systems. 
Includes  material  on  sampling,   signal 
recovery,   compensation  and  statistical 
theory  of  quantization,    z  =  e  ~sJ-s 

62.  Julius  J.    Tou,   Digital  and  Sampled  Data  Control  Systems 
(New  York,  1959). 

Sampled  data  text.     Good  section  on  compensation, 
especially  with  digital  computers,    z  =  e        s. 

63.  John  G.    Truxal, Automatic  Feedback  Control  System  Synthesis 
(New  York  1955). 

Standard  text  in  the  field  covering  practically 
all  aspects  of  automatic  control. 
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48.  Alton  S.   Householder,   Principles  of  Numerical  Analysis, 

(New  York,  1953). 

Numerical  methods  for  sets  of  linear  equations, 
differentiation  and  integration.     Excellent 
bibliography  and  bibliographical  notes. 

49.  Harold  E.   Maurer,   Error  Analysis  of  a  Digital  Differential 

Analyzer,   E.  E.    Thesis,   M.I.  T.   (MIT/IL  Rpt.    T-129),May,   1957. 

Attempt  to  use  Widrow's  (Ref.    5  5)  statistical 
characterization  of  quantization. 

50.  Harold  E.   Maurer,  An  Approximate  Analysis  of  Error  Propagation 
in  a  Digital  Differential  Analyzer  (MIT/IL  Rpt.   E-698),   March,   1958. 

Extension  of  Ref.   49. 

51.  John  M.   Salzer,   Treatment  of  Digital  Control  Systems  and 

Numerical  Processes  in  the  Frequency  Domain,   Sc.  D.    Thesis, 

M.I.  T,  ,  August,  1951. 

A  clean  way  to  handle  numerical  analysis. 
Presents  a  design  method  for  computers. 

52.  James  B.   Scarborough,   Numerical  Mathematical  Analysis,    4 

ed.   (Baltimore,   1958). 

Standard  work  on  numerical  analysis.     Full 
of  numerical  examples. 

53.  Wallace  E.   Vander  Velde,   "Computers  in  Inertial  Navigation 

Systems",    Section  V  in  Inertial  Guidance  -  Terrestrial  and 

Interplanetary,   (Course  Notes  for  16.  45S,   M.  I.  T. ,   Cambridge,   Mass., 

Summer,  1960). 

A  good  basic  discussion  of  all  three  types 
of  computers,  their  operative  advantages 
and  disadvantages. 

54.  E.  T.   Whittaker  and  G.    Robinson,   The  Calculus  of  Observations, 

4th  ed.    (Glasgow,  1958). 

Covers  most  numerical  analysis  topics  as  well 
as  curve  fitting  and  Fourier  series. 

55.  Bernard  Widrow,  A  Study  of  Rough  Amplitude  Quantization  by 

Means  of  Nyquist  Sampling  Theory,   Sc.  D.    Thesis,   M.  I.  T. ,   June,   1958. 

A  detailed  analysis  of  quantization  through 
a  statistical  model. 
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64.  John  G.    Truxal,   ed.  ,    Control  Engineer's  Handbook  (New  York, 

1958). 

As  the  title  implies,   this  is  written  for 
the  practicing  engineer.     Comprehensive 
coverage. 

65.  Albert  D.   Wheelon,    "Free  Flight  of  A  Ballistic  Missile, 

ARS  J.,   XXIX  (1959),    915  =  926. 

One  of  the  few  places  that  range  derivatives 
are  discussed  in  the  open  literature. 

66.  Henry  F.   Blazek,    "The  Performance  of  Inertial  Components 

on  an  Unstabilized  Base",   Presented  at  ARS  Controllable  Satelites 

Conference,    M.I.T.,    30  April  1959. 

A  brief  listing  of  the  effects  of  unstabilized 
bases  on  instruments. 
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